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ABSTRACT 

Atmospheric atomic sodium has been observed over Urbana using a laser 

radar system tuned to 589 nm. The layer typically has a peak concentration 

9 10-3 

of 10 -10 m at 90 km and is about 10 km thick. The photocount data are 
processed using a digital filter to obtain continuous estimates of the sodium 
concentration versus altitude. Wave-like structures in the sodium layer have 
been observed with typical wavelengths of 3-15 km and doimward phase velo- 
cities of 1-3 111 s”^. On one occasion there was evidence for the presence of 
a standing wave in the layer. During two successive evenings the bottomside 
of the layer was observed to undulate with a period of about hours. In 
agreement with worldwide measurements the layer is observed to broaden through 
the night. 

A meteor ablation-cluster ion theory of sodium is developed £uid found 
to be in good agreement with existing atmospheric observations as well ns 
laboratory measurements of rate constants. The Na^ chemistry is connected 
to the Na chemistry by cluster ion reactions involving Na^'N^, Na^*C0.,, and 
Na^‘H.,0. The chemistry includes switching and dissociative recombinatioi of 
cluster ions but does not include their photodestruction. The chemistry i.s 
then used in conjunction with a meteor ablation source and a downward sweeping 
of ionisation caused by an interaction of the neutral wind profile with the 
geomagnetic lines of force (corkscrew mechanism) to explain many features of 
the sodium layer. Sodium is removed from the mesosphere in the form of 
higher order cluster ions which eventually undergo sedimentation to the lower 
atmosphere. 
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1. INTRODUCTION 

1.1 Post Cbs&n)atidnB and Reautta 

Since the discovery of atnos^eric sodium emissions in the late thirties, 
several investigators have measured the features of the sodium layer and 
theorized about its origin and variation, A review of this activity is pre- 
sented in the publications by Bnotan [1973] and Kvifte [1973] . In spite of 
the varied approaches in both ground-based and rocket -borne experiments to 
learn more about the sodium layer, many of the reasons for its origin and 
behavior remain elusive because important correlative measurements have not 
been made. Yet kr.ow ledge of the ihysics and chemistry of atmospheric sodivun 
continues to be important for the understanding of the aeronomy of the middle 
atmosphere. In particular the chemistry of sodium is intimately connected 
to sudi Important constituents as atomic oxygen and ozone as well as elec- 
tron density. Furthermore, it is felt that sodium may be used as a tracer 
for atmospheric motions in the 80-100 km region. 

Prior to 1970, sodium measurements were largely restricted to measuring 
the total flux of resonantly scattered sunlight, giving only the total column 
density of sodium. Rocket measurements [Hunten and Wallaoe^ 1967; Donahue 
and Meier, 1967] were the first to reveal the existence of sharp upper and 
lower boundaries for the layer. Since the initial lidar measurements of 
sodium [Boman et at, , 1969] , nighttime monitoring has established several 
important features of the layer. These measurements combine the altitude 
resolution of rocket experiments with the continuous coverage ability of 
previous ground-based measurements. 

1.2 Lidar Obeervations and Reeutte 

Remote sensing of the sodium layer using a dye laser tuned to one of 


2 


the sodium resonance lines » usually 589,0 nm, has produced detailed night- 
time profiles. In general it has been found that atomic sodium consistently 
resides in the upper part of what is referred to as the middle atmosphere in 

the altitude range of 80-100 km. Its distribution is rather sharply peaked 

9 10-5 

at about 90 km, with a maximum number density of 10 - 10 m . The scale 
height at the top of tl\e layer is about 3-4 km, much smaller than that of the 
surrounding neutral atmosphere; and the lower portion of the layer falls off 
quite rapidly with an order of magnitude drop in about 4 km. 

However, it is the desire to understand the variations of the layer which 
has stirred so much Interest in the sodium problem. Specifically there is a 
winter enhancement at mid-latitudes of about a factor of two and a lowering 
of the altitude of the layer bottoms ide. Tlvis v;ell-defined seasonal varia- 
bility is illustrated by Figure 1.1 [(Jibson rnitf Sandfordt 1971] and also 
appears in the sodium column abundance shown in Figure 1.2 1975]. 

Although ground-based ffSatfaden and Purdti/j 1970] and rocket -borne [Donahue 
and Mo'iex'j 1967] dayglow measurements have indicated a diurnal variation, 
laser radar measurements have indicated no such variation [ffibson £wd 
Sandfoi'd, 1972] . The nighttime sodium concentration shows little systematic 
variation during the night [(Jtbeott and Sccndfovd, 1971] although Hake et a2. 
[1972] have noticed an Increase in the sodium layer during a period of the 
Geminids meteor shower activity. Detailed structure was observed in the 
sodium layer by Kivahhoff and Clemesha [1973] and they attributed this struc- 
ture to atmospheric tidal oscillations. Negi& [1975] reported on 80 nights 
of observations from October 1973 to July 1975 and speculated tltat dynamics 
of the mesopause region were responsible for the time -dependent stratifica- 
tion of the layer. Using a steei'able sodium lidar, Thoma et at, [1976, 1977] 
claim to have observed lateral movement of sections of the sodium layer. 
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Hence, a wealth of infoTtnation about the sodium layer has been provided by 
ground-based lidars, but as yet there is no theoretical study which can simul- 
taneously explain all of the observed features , 

1,5 Theories of the Sodium Layer and Problem 

Detailed theories of the sodium layer have been presented by Hunten and 
f/aZlaae [1967], Donahue and Meier [1967], Sanson and Donaldson [1967], 

Gadsden [1968], Hunten [1971], Nitra [1973], and Hoaao and Visaonti [1973, 
1974] . A problem common to all of these theories is that none of them can 
cxpliiin all of the observed features of the sodium layer. It is the attempt 
to explain these observations which has prompted the theoretical investiga- 
tion of this work, and this theory is bx'iefly summarized below. 

Gadsden [1970] made the interesting suggestion that sodium atoms ablated 
from meteors are ionized in the E region and the long-lived Na ions move 
do\<m to the base of the E region by the corkscrew mechanism [Chimonas cmd 
Axfordj 1968] . Follottfing their transport to the 80-100 km region, the Na* 
ions are converted to neutral Na tlirough cluster ion reactions. In this way 
the dumping of ions by the coi'kscrew mechanism maintains the steep topside 
gradient of the sodium layer in the presence of vertical eddy and molecular 
diffusion. Within the sodium layer in the 80-100 km region, chemical quasi- 
equilibrium exists with the predominant forms of sodiunil being Na, NaO, Na*, 
and Na*H,G. It is questionable whether other forms of sodium compounds such 
as NaO,, and NaO* play an important part in the chemistry. Nevertheless it 
may be shown that sodium concentration is proportional to both electron den- 
sity and atomic oxygen. Hence, the sharp bottomside falloff of sodium is 
probably a result of the steep electron and o^Q^gen density gradients below 
90 km. Below 85 km sodium is removed from the layer, possibly in the form 
of higher order water clusters which undergo sediment ation to the lower 
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Mesosphere. The details of this neteor-ion theory are presented in Chapter 7, 
1.4 Qb^&ativ&s afui Saopa of Tftis Study 

The scope of this work has included an experimental as well as theoreti» 
cal investigation of the sodium layer. As part of this work a sodium Udar 
system was designed, constructed, and operated at the Aeronomy Laborator)*^ 
Field Station, Altltouglt the published system cliaracteristics of many groups 
are rather sparse, some comparison between lasers is jmssible and these 
characteristics are tabulated in Table 1.1. A similar coi^arison between re- 
ceivers is impossible because of lack of data. 

As indicated in Table 1.1, the Urb.ana lidar is quite adequate for the 
first scientific measurements. In addition, it should be noted that the 
higher energy lasers employ an oscillator-amplifier configuration and an 
amplifier would be jm excellent addition fbr future correlative measurements. 
Tlie theoretical work was developed simultatveously and conclusions were drawn 
about tlve experimental observations. 

A brief review of terminology and principles involved with the study 
of sodium v?ipor is presented in Chapter 2. Chapter 3 provides a detailed 
discussion of the upper atmospheric lidar applicsition which is appropriate 
to this study. 

lit^ortant concepts about dye lasers and a description of the IJrbana 
laser fbllows in Gluqjter 4. The Urhana receiving SN’stem and routine oper- 
ating techniques are discussed in Chapter 5. 

The experimental results are presented in Quapter 6. The feasibility 
study for tlie first correlative measurements with tl»e Urbana meteor radar 
facility is also included in this cliapter. 

Chapter 7 develops the meteor-ion tlieory. This study inwlves ,a mecha- 
nistic approach to the problem. That is, individual dvnamical and chemical 
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TABLE 1.1 


Summary of dye lasers used for sodium lidar measurements 


l^nstitution 

Wavelength 

Energy/Sfaot 

Rtitietition 

I^te 

Lineuldch 

aiiJ Space Research Station, liNGlANO 
e.t 1969] 

S89.6 nn 

4 vl) 

U.2 s"^ 

0,15 iiM 

GIJA Corporation Technology Division, 


600 mJ 

— 

0,i-Q,i im 

Bedford, Massadiusetts 
ISdiwZer et 19711 


200 mJ 


12 pm 

SRC Radio and Space Research Station, BNGhAND 
{Oihacm and Sandfo^*dy ^1012] 

589,0 

200 mJ 

2 s'^ 

3 p» 

Service d^Aerotionic du CMRS , FRAMCK 

$89,0 lun 

125 mJ 

O.U17-O.OS s', 

IS pm 

[Blairmt fct ttZ. > 1972); 

S89.0 nm 

300 stl^ 

0.017'0.0S s 

13 pm 

Ciiropean Space Research fnstituce, lTAl*y 
and ScIm&idei'^Murttau, 19721 

S89,0 DA 

600 nJ* 


10 p« 

Stant^ord Research institute, 
Hcnlo Park, California 
iHdket at 1^72} 

589,0 m 

SOO wJ* 

.25 s'* 

<5 pm 

Instituto de Pesquisns hspaciais, BRAZIL 
[Kirahlioff and Clemstha, 1973] 

589,0 nm 

5 mJ 

0,«S5 s"* 

5 pm 

Tohoku IhiivcTslty , JAPAN 
el «l.j 1974] 


2 mJ 

1 s'* 

\^l\ 

2 pm 

Service d'Aeronwaie du CNUS, FRANCE 

569,0 nm 

400 tk? 

10 pm 

[U^lh and Neffie, 1974} 

589,0 nm 

1.2 J* 

1 s 

to pm 

University of tlllfiois, 

Urbanu, Illinois 

♦Denotes an oscillator-ampiif ier configuration. 

589.0 n» 

iOO mJ 

.5 

10 pm 




*sl 


I 


[GINAI* @ 




nechanisns are investigated by nunericai methods and the results are inter 
preted in the light of experimental findings. 

Chapter 8 sunnarises the eiqierimental and theoretical iiqslications of 
this study. Future areas for research are defined and further correlative 
measurements are suggested. 
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2. ATOMIC SODIUM VAPOR CHARACTERISTICS AND SPECTROSCOPIC NOTATION 
2,1 Basic Tvansitians 

Sodltim is the eleventh element of the periodic table and exists at meso- 
spheric heights in the atomic gaseous phase. Being an alkali atom, its structure 
consists of one valence electron and two inner shells of non-radiating electrons. 
A brief review of spectroscopic notation applicable to the sodium atom is helpful 
in understanding some basic concepts important to this lidar experiment. 

It is necessary to designate a number of spatial quantum numbers to describe 
the complicated spectrum of sodiiun. Among these we include the principal or 
total quantum number, n, and the angular momentum quantum number, I. The 
principal quantum number takes on integral values (1, 2, n, ...) and denotes 

the orbit or shell in which the electron is moving. For a given value of Z 

may take on values 0, 1, 2, ... n-1, and characterizes the orbital angular 
momentum of the electron. The states associated with the Z values are sometimes 
referred to as sub-shells. In spectroscopic notation it is common to designate 
these sub-shells by the letters, e, p, d?, /, (altematively 5, P, D, F) corre- 
sponding to the values Z = 0, 1, 2, 3, respectively. The notation grew out of 
optical spectroscopy where sets of spectral lines were characterized as sharp* 

principal, diffuse, and fine. Using this notation the electronic ground state 

2 2 6 1 

configuration for sodium may be written as Is 2s 2p 3s . The number preceding 
a letter refers to the total quantum number of state, n, while the superscripts 
denote the number of electrons in each subshell. 

Sodium, like all alkali atoms, exhibits a series of doublets rather than 
single lines in its spectrum. To describe this phenomenon, we must introduce 
a new quantum number, This total angular momentum quantum number for 
extranuclear electrons is obtaiued by the expression: 


lid 

z » 0 


10 


( 2 . 1 ) 




1*9 

1/2 


where a is the spin angular momentum number and takes on values a t 1/2. Note 
that the Z » 0 state is not split. Hence, all states except the S state are 
actually double. 

We may then label our states as shown in Table 2.1 where J appears as a 
subscript in the entries and the superscript 2 denotes a doublet state. The 
selection rules for allowed transitions become Aj - ±1, 0 with the transition 
3 - Q 3 =0 not allowed and Aj = ±1, 0 with Z = 0 ^ Z = 0 forbidden. 
Therefore, the two D-line transitions of interest to our problem are: 



X = 589.6 nm 

Cz?i 

line) 

(2.2) 


X ~ 589.0 nm 

02 

line) 

(2.3) 


These lines appear in the energy level diagram shown in Figure 2.1. 

2.2 Spectral Properties of High Resolution Linas 

2.2.1 Hyperfine structure. If one observes the line spectrum with an 

rt 

apparatus of high resolving power, one will find a fine grained or hyperfine 
structure to the lines. The hyperfine structure is qualitatively explained by 
assuming the nucleus possesses an intrinsic angular momentum and a corresponding 
magnetic moment. We assign another quantum number, i, to denote this momentum. 

The quantum number i differs from one type of atom to the next. For sodium 
i ^ 5/2, A new quantity, f, called the hyperfine quantinn number is defined to 
represent the total angular momentum and may take values f- 3 +i, 3 + i- l, 

3 + i - 2, |j - i|. If 3 ^ i, there are 2i + 1 values for f line). 

If i - 0 , there are 23 + 1 values for / (Z?j line). The selection rules that 
apply to f are = +1, 0 but / = 0 = 0 is forbidden. 












y: 


J=i-. ; - 1 . 
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2 2 2 

We observe splitting at the ^1/2’ ^3/2 because the 

same interval rule holds as did for A,/, the hypermultiplet will have a similar 
appearance to an ordinary multiplet. Hence, each Z7->line has two groups of fine 
stinicture lines associated with it and we designate these groups as ^-^a* ^12>' 
^2a* ^2h’ splitting applicable to the line is shown in Figure 2,2. 

2,2.2 Zeemm effect. Because the atom has a magnetic dipole moment, the 
splitting of the hyperfine lines can be observed when the atom is placed in an 
external magnetic field {Zeeman effect) . We may estimate the effect of the 
earth's magnetic field on this splitting by the expression 1934]: 


Ao = m 


g / eg \ 
\4ir mcj 


CSI units) 


(2.4) 


where m. is the magnetic quantum number and denotes the sublevel of hyperfine 
structure and g is the Lande g'-f actor. The quantum number, m., takes values 

tJ 

differing by one from -j’ to +j. Lande g^-values depend on the quantum numbers j 
and Z and if knovm, are usually included in tables of atomic energy states 
[e.g., KuhUi 1969, p. 291]. The values for m. and the product m. • g for our 
case are listed in Figure 2.2. The quantity enclosed by parentheses is called 
the Larmor frequency. Quantities indicated have their usual meaning, namely: 

e = electronic charge 
B - magnetic field 
m = mass of electron 
<2 = Speed of light 

_5 

If we assume the earth has a magnetic field of 6 ^ lO Tesla (.6 Gauss) at our 
latitude, then we calculate: 


Ao = S X 10"^ 


(2.5) 


6 

Av = oAo = 1.5 X 10 s" 


( 2 . 6 ) 











^ s 2. This is more than 2 orders of magnitude smaller than the 


IS 


for • g s 2, This is more than 2 orders of magnitude smaller than the 
separation of the hyperfine lines. Hence, it is appropriate to neglect the 
splitting of the line due to the earth’s magnetic field. 

2 . 3 Line Intensities 

Another principle, important to the sodium spectrum, relates the relative 
Intensities of the lines within a given multiplet. The statement, known as 
the sura rule of Omstein, Burgers, and Dorgelo, is as follows: The sum of the 

intensities of alt the lines of a imltiplet whiah belong to the same initial 
ov final state is proportional to the statistical weight 2 q + 1 of the initial 
Or final state , respeotively , Hence, if we let I represent the intensity of 
the line, then: 


D 


2^2 (3/2) +1 2 

■ 2 ( 1 / 2 ) + 1 " 1 


(2.7) 


Hence, it is preferable to tune the system to the line for a twofold 
increase in returns . 

Furthermore, this rule applies to the hyperfine structure if we replace 

2 

3 by /. For example, the ratio of the sum of intensities with ^ 

are to those with = 2 as (2 x 3 + l)/(2 x 2 + 1) = 7/5. Referring 

to Figure 2.3 we note that the only line from the / = 3 state 
2 2 

( ^3/2^ = 3 ^if 2 ^ ^ ^ relative strength of 14. There are two lines 

coming from the f = 2 state (^^ 3 / 2 ^ = 2 = 1 and = 2 -*• = 2) 

both having a relative strength of 5. The ratio of the intensities is therefore 
14/(5 + 5) = 7/5 in agreement with the sura rule. 

All of the relative intensities may be obtained by this rule. The line 
spectrum as shown in Figure 2.3 would be that of sodium at 0®K. However, 

Doppler broadening is dominant for higher temperatures, and this tends to merge 
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SODIUM RESONANCE LINE SHAPES AT DIFFERENT TEMPERATURES 


Figure 2.3 Line profiles for and Z>2 emission at various 
temperatures [Chani>erlain, 1961], 
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the lines into a continuous absorption profile as also shown in Figure 2.3. 

The profile for 200" K is suitable for mesospheric heights corresponding to 
the sodium layer, 

2.4 Absorption Coefficient 

2.4.1 liaturaZ and Doppler broadening. There are several processes that 
contribute to the formation of an absorption line of a gas. For mesospheric 
heights we need only consider the processes of natural broadening and Doppler 
broadening. 

Natural broadening results from the finite lifetime of the excited state. 
For sodium the lifetime of the "*^ 3^7 state is x ® 16 ns [Mitchell and Zmansky, 
1971]. Hence, we calculate: 

Av = ^ = 10^ Hz (2.8) 

n 2tTT 

We will see that this is negligible compared to the effect of Doppler 
broadening . 

Neglecting natural damping, the expression for Doppler broadening of a 
line is given by: 

ln 2 . fT 
AVp = — 5 ''0 V W 

where 

R - universal gas constant 
c - speed of light 
Vp = center line frequency (v^ - 
T - absolute temperature of the gas 
M = molecular weight 

Taking = 589 nm and T = 200"K, we calculate for sodium: 

= 1.07 X 10 ® Hz 


(2.9) 


( 2 . 10 ) 
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Clearly we have Av^ » AVj^^. 


2.4.2 Mcaeimm dbsoi^tion ooeffiaient^ kQ» We will follow closely the 
excellent work by Mitohell and Zemanaky [1971] in developing the expression 


for the absorption coefficient. If we consider only Doppler broadening, the 


absorption coefficient of a gas is given by: 


\ = % '"p - 


2/lii2 (v - Vq) 


( 2 . 11 ) 


where represents the maximum absorption coefficient. Integrating equation 


(2.11) we obtain: 


lo In 2 


( 2 . 12 ) 


Two alternate forms for this expression exist due to Ladenburg: 


r.. 

Jn 


S', p 


8tt T 


(2.13a) 




(SI units) 


(2.13b) 


where 


= statistical weights of the normal and excited states 


respectively 


p = niiraber density of sodium 


T =5 lifetime of the atom in the excited state 


f = oscillator strength or spectral line /-value 


e = electronic charge 


e^= permittivity of free space 




■W 
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m = mass of an electron 


a B speed of light 


If we define the scattering cross section, at the peak of the line by 
the expression = cr^p, then by combining equations (2.12) and (2.13) we may 


solve for o^ and obtain: 


1 2 0 ^ 2 ^ ^ 
n 8 tt T 


TT 4eQnw 


f (SI units) 


Taking f - ,67 [Mitahell and Zemanskyj 1971] we calculate: 
cTq * 1.5S X lO’^^ m^ 


(2.15) 


2.4.3 sti'*cngths and positions of and lines. Since the 

Z?, line on a hyperfine structure level consists of 6 separate lines, we may 
write a more exact expression for k^. Assuming all hyperfine levels of a 
resonance level have tht? same lifetime, we express k^ in sura form as: 


6 2 /In 2 (v - v^) 

k = y exp - 7 — ^ ^ 

^=1 


(2.16) 


where the siun is over the 6 hyperfine lines. 

However, for tliis study, it was found sufficient to use just two lines 
designated as the and lines. The positions of tliese lines are taken 
as the weighted averages of the appropr-iate hyperfine lines. Referring to 
Figure 2.5 we obtain a wave number spacing for the two lines of: 


14(-2.07) . 5(-:.27) . H-2. S9 X , .j.„ „-l 


(®2n 'I"®) 
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3. SODIUM LIDAR PRINCIPLES 

3.1 Inti’orfucjtion 

Lidar (light tfetection and ranging) systems operate in much the same 
manner as conventional radar systems. An alternative name for lidar, laser 
radar, is used to emphasise the similarities. Essentially only the 
operating frequency has been raised. 

Uowe\fer to understand the operation of a lidar system it is necessary 
to cover several topics pertinent to optical atmospheric transmission as 
well as to laser characteristics. In Section 3,2 we will discuss the 
different tjqies of optical scattering which occur in the atmosphere. 

Sections 3.3 and 5.4 will develop the lidar equation for application to the 
sodium layer. Section 3,5 will present the technique of system calibration, 
A design calculation for system analysis will be presented in Section 3,6. 
And finally a check on the possibility of line saturation will be given in 
Section 3*7. 

5.2 Pvoaessi^s 

A laser beam passing through a homogeneous atmosphere of thickness s 
is attenuated and Bouguer*s law applies as follows: 

I = Ip ejq) - [6s] (5,1) 

where Ip is the incident intensity and P is the extinction coefficient 
having units m'^. In general the e.xpression is more complicated since the 
atmosphere is not homogeneous and hence p = p(s) and equation (3.1) becomes 


I = Ip e.xp 



P(e)Ja 




(5,2) 


! 

I 


Tlie quantity $ is usually tlie sum of several compOTients. Lower 
altitude lidar retvirns cleai'ly indicate that the main contribution conies 




wm 


from Rayleigh C^olecular) scattering. Occasionally layers of dust or high 


altitude clouds may cause deviations from molecular scattering. Mie 


scattering, or scattering from aerosols, is responsible for this effect. 


Both Rayleigh and Mie scattering are elastic processes for which excitation 


can take place at any frequency. 


On the other hand, several inelastic scattering processes which are 


wavelength specific may be used to distinguish constituents. Raman scattering 


occurs when the incident light excites vibrational or rotational transitions 


of molecules. Re-emission takes place on either side of the incident 


frequency aiid the shift characterizes the molecule involved. The efficiency 


of the process may be enhanced by as much as a factor of if the incident 


radiation is close to an electronic transition. The process is sometimes 


referred to as resonant Raman scattering. In either event the process plays 


a minor role in the value of the extinction coefficient for this experiment. 


Fluorescence involves the absorption of the incident radiation in one of the 


absorption bands of an atom or molecule and the re-emission in a longer 


wavelength band of the spectrum. The current literature draws little distinc- 


tion between the fluorescence and resonance processes. Resonance scattering 


results from the excitation of electronic transitions . In this case the 


re-emitted light is the same frequency as the incident light and occurs after 


a period of time corresponding to the lifetime of the excited state (1-20 ns). 


Since the process is very efficient, it may be necessary to calculate 


the attenuation of the beam as it passes through the sodium layer 


to compute accurately the densities at the top of the layer. Hence resonance 


scattering may contribute to the extinction coefficient. The final 


contribution to the extinction coefficient comes from absorption. In this 


case the incident energy is not scattered but rather converted to a 




non-optical for«i of onergy, \ culculutlou is wmle later to estimate absorp- 
tion by osono. 

The extinction coofficiont is proportional to the density, p, of the 
scattering medium; and often we express the extinction coefficient for a 
particular process In terms of a scatterliig cross section, o, by the relation 


[1 «* po 


(3.3) 


The scattering cross seetloti is ossentlally the cross sectional area that 
the particle presents to the incident beam for scattoting. As a guide the 
scattering cross sections for those various processes are compared In Table 
o.l. It is apparent from Table 5.1 that resonance scattering may be used 
to detect very small densities of a coi\stltuent since its scattering cross 
section may he 10 orders of magnitude greater titan the other processes 
involved, 

3.5 tiidfW b'ipkttion 

5.3,1 l'\ymuiaiiont The University of Illinois lidar Is a pulsed system 
operating it» the monostatic mode, Ttu' laser is tuned to a resonance line of 
sodium and the received backscatter signal is processed to yield Information 
about the sodium density profile. This translation of the received signal 
into s density profile requires a knowledge of the Interaction of the laser 
pulse with the atmosphere. The system configuration is Illustrated in 
Figure 5.1. The laser transmitter and receiver are separated slightly 
t'dS m) but the separation is small compared to the altitude of the sodium 
layer ('HO km) and can ttierefore be neglected. TIve solid angle of the 
receiver, is detetmined by the field stop iris and the diffraction limit 
of the Fresnel objective lens. The i'resnel lens limits the minimum half 
angle field of view to about I mrad which is comparable to the laser diver- 
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Table 3,1 Typical Values 

for Scattering Cross Sections 

Process 

Cross section Cm^) 

Mie 

10”^'^ - 10"^^ 

Raman 

10"^^ - 10“^"^ 

Rayleigh 

10-52 

Res onance* Fluores cence 

10"^^ - i0"23 


OKIGmAL 

Of* POOR QUAUW 
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gence. Consequently^ the field stop iris is opened slightly so that the 
entire scattering volume is Illuminated by the laser ((1^ > . 

The laser radar equation Indicates the system and atmospheric parameters 
which contribute to the measured signal. Since the received signal is 
small, it is necessary to operate in a photon counting mode. A few remarks 
concerning the derivation of the lidar equation for this high altitude 
application will be helpful. Excellent references for the following are 
Novthend et al, [1966], [1972], and ilfeg’ie [1976]. 

We assume that the laser emits an expected value of photons per 
pulse with a frequency distribution p^(v) . The quantity giy) represents the 
fraction of photons emitted between v and v + cfv and by definition we have 
! 0'Cv) = 1. For simplicity we do not consider the temporal variation of 

Jo 

the pulse at this time but rather only assume that the duration of the pulse 
is less than 2hz/o where da is the range bin size. 

The energy per photon is hv, where h is Planck’s constant and v is the 
frequency of the emitted photon. If P and i represent the power and 
duration respectively of the transmitted pulse, then we may write 


,, Pt 


(3.4) 


This expression is useful for approximate calculations; however, the frequency 
dependence of the pulse, ^(v), is ultimately needed and hence a spectral 
analysis of the laser pulse must be performed. 

Consider the photon flux incident on the bottom of the layer which 
occurs at altitude Zq (Figure 3.2) . If we lump all of the atmospheric 
attenuation which results from elastic scattering (Rayleigh and Mie scattering) 
and absorption into a parameter T, then we may express this flux as 




^ V " ' .r 



INCIDENT LASER BEAM 


Figure 3.2 Range cell convention. 
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where Is the surface illusdnated by the laser beam at altitude a. 

This area is given by the expression 

A^Cs) « flj a (3.6) 

where is the angular halfWidth or divergence angle of the laser beam. 

For a 1 mrad the beam has a radius of about 80-100 m at the altitude of 

2 

the sodium layer and so A^(3 q) = 20000 m , 

As the beam passes through the sodium layer, it is also attenuated by 
sodium 27 line resonance scattering. Hence at altitude s > the flux 1$ 
given by 


21? ffCv)ir 


i iX-) * - 


(3,7) 

where is the optical depth of sodium given by 





T^(^) = 

a^p(e)d8 

(3.8) 


In expression (3.8) is the absorption cross section for sodium and pCs) is 
the density of sodium at altitude s. (Note that cross section 

is not constant with altitude as it depends on temperature. 

We may approximate the scattering atmospheric volume by a cylinder 
since the divergence angle for the laser is small. Hence the number of 
photons re-emitted per unit area from an infinitesimally small thickness, dz, 
is given by 

d2l?(v,3) = 0^p(3)4»^(3)d?3 


(3,9) 
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We are assuming in equation (3 <9) that natural damping is more important than 
deactivation by collisions with atomic or molecular constituents. This is 
reasonable since the collisional frequency at 80 km is about 10^ s ^ and the 
natural lifetime of the excited state is 16 ns. In addition we may assume 
that the re- emitted photon has the same frequency as the incident photon. 

If the area of the receiver is then the receiver solid angle as 
seen from the atmospheric cell is given by 




1’ 



C3.10) 


Assuming that the scattering from the sodium layer is isotropic, the number 
of received photons from the partition of the sodium layer ds thick is 


fv,a) = T exp - [t fail (3.11) 

where the factor of 1/4 tt converts the isotropic scatterer to a steradian 
basis. 

Combining equations (3.7), (3.9) and (5.11) we may write the lidar 
equation for our case. If we define 

s. = a„ + I'Aa i - integer (5.12) 

t u 

then the expected value of the number of photoelectrons, ff.(v), counted 

T* 

during tl\e interval of time corresponding to the returns from the Cth range 
bin is given by 




3 




eM. 


o.. exp - [2T^^(e)]ds 

s 


(3.13) 


where n is the receiver efficiency. 
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3.3.2 Approxinjationa, A number of approximations may be used to sim> 

2 

plify the integration in (3.13). The factor of 1/e may be approximated by 

2 

the constant 1/ (a . + As/2) . This is equivalent to assuming that l/a changes 

% 

slowly over distances on the order of the range bin size, and the approxima- 
tion introduces a maximum range error of Lz/2z. For As = 1.5 km this means 
an error 1% at 80 km and 0.75% at 100 km. 

The exponential factor in (3.13) may also be approximated by the first 
two terms of its Taylor*s series e^qiansion. If we take the column density 
of sodium to be 10^^ m"'^ and let the scattering cross section be 9 x lO”^ m^ 
which corresponds to the peak of the line, then 

s m"^)(9 X 10’^® m^) = 0.09 (3.14 

Hence the error incurred by replacing exp - [2t^Cs)] by 1 - 2tv(s) is less 
than 2%. 

The absorption cross section o^, varies with altitude since it is tem- 
perature dependent. This is because the expression for is inversely pro- 
portional to A^p and hence inversely proportional to y/¥ (see equations 2.9 
and 2,11). The relative error that results from neglecting temperature varia- 
tions depends on the frequency v of the incident radiation. For example at 
V = Vq, the dependence of the exponential factor in on Avp is eliminated. 
The temperature profile at this altitude can be highly variable and exhibit 
a wavelike structure at times. Ten5>eratures can vary between 150-250*K. For 
T - 200 ± 50*K, 0 ^ varies by at most t 15%. However when the laser linewidth 
is greater than the width of the absorption line, the effective scattering 
cross section is independent of temperature (see Section 3.4,2). For this 
reason cr^ will be assumed to be constant with altitude. 

It is appropriate to discuss the consequences of neglecting the time 
response of the laser. If represents the photon time response of the 
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laser pulse, then the total intensity received from the passage of the laser 
pulse through the atmosphere is given by 


P(s)W^(* - ds 


(3.15) 


where F is a fraction which depends on lidar parameters such as number density 
of scattering medium, cross section of scatterer, range, geometry, and 
efficiency factors. 

To obtain the number of counts in a bin, the signal processing equipment 
counts the number of photoelectrons emitted from the cathode of the photo- 
multiplier tube in an interval of duration T = 2Ls/a, Hence for a particular 
bin, the number of counts recorded is 
ft 


d. 

•' t-T ^ 


FCaJff^Ct - -=^) da 


(3.16) 


where the lower limit on the a integrand has been extended to by using 
the convention that Fis) = 0 for a negative. It is convenient to rewrite 
(3.16) in terras of a receiver gating function G, 

0 £ t 5 T 
otherwise 

Then expression (3.16) becomes 


°{o 


(3.17) 


F(a)(i, 


=f: 




dx 


(5.18) 


where = 2s/a is the delay between the transmitted and received pulse. 
Recalling that F(s) is proportional to the density of the scattering mediumi, 
p( 3 ), we find that (3.18) is a relationship involving the two-fold convolution 
of the density, transmitter pulse, and receiver gating function. 

If the receiver gate is larger ti\an the pulse width* the system resolution 
is approximately As = oT/2. The characteristic resolution is the value which 
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results from the convolving of the laser pulse with the gating function and is 
tHe sum of the range bin size and spatial length of pulse. We may estimate 
the error incurred if we neglect this convolution. An approximation to 
based on photographs of the time response of the laser pulse (Figure 3.3) is 




2c? t exp - {at}? t 2 0 


(3.19) 


otherwise 


where ff.(t) is normalized, a « 1/3 and t is in us (Figure 3.4). The full 


width at half maximum is about 3.5 us which agrees with experimental values. 

If we neglect the translation t^, then we may consider the convolution 
to be of the form 


V' 


-w.r 


u(t)w(t - i)dT 


(3.20) 


where 


v(*) - 


2 2 
2a t exp - [at] 


m(*) = 


0 

1 

0 


(3.21) 

(3.22) 


t > 0 
otherwise 
0 < t < T 
otherwise 

or convenience we subdivide the problem into three cases. 

Case 1: r < 0. In this ease the two functions do not overlap and their 

product is zero. Hence 

D*b} = 0 t < 0 (3.23) 

Case 2: 0 i t ^ T. The functions partially overlap but only over the range 

0 ^ T S t. Hence 
ft 


v*w 


7 2 

= 2a‘^T exp - [ar] dt 

Jo 


(3.24a) 
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Figure 3.3 


(a) Five superimposed laser pulses as seen by a 
sodium cell monitor. Horizontal scale is 2 ps/div 
b) Five superimposed laser pulses shown as tuning' 
dritts to. the edge of the line. Horizontal scfle 
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This integral nay be solved by making the substitution z ~ (axy and we 


obtain 


v*w « 1 - exp - [at\ 


0 < * < S’ 


(3.24b) 


Case 3: t i T, In this situation the functions fully overlap and thus 


= f 2a^T exp - [at]^ dx 
it-T 


= exp - [< 2 ^(t - 20^] “ exp - [at]^ 


(3.25) 


t : T 


This convolution is plotted in Figure 3.5. Note that when we consider the 


convolution of the laser pulse and the gating function, there is some overlap 


of signals from adjacent range bins. 


A estimate may be made on the effect of neglecting this overlap. The 


sodium layer can fall by over a factor o£ 10 in 4 km on the bottomside and 


has a scale height of 3-4 km at the topside. This means that adjacent 1,5 km 


range bins at the bottomside may differ by as much as a factor of 4 and by a 


factor of 2 at the topside. Referring to Figure 3.6, we see that the total 


counts received during the sampling corresponding to bin i consists of a tail 


from bin •£ - 1 and the counts received from bin i while the laser pulse is 


entirely within the range bin. This return .s to be compared to the total 


number of counts received form the r-^bin. Referring to Figure 3.6b, if we 


assume that the density of scatters in range bin i - 1 is a factor of 4 


smaller than that of bln i, then we form 


# counts received duting interval of duration T 


+ 4^2 


counts emitted from range bin % 


4(^2 + dj) 


= 0.80 


(3.26) 


where we are integrating the expressions derived in equations (3.24b) and 


(3.25) to determine A^ and A^ respectively. Obviously A^ - Ay Alternatively 


1 ^ 
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(a) Space-time diagram of propagating laser pulse; 

(b) overlap of scattering from adjacent range bins 
caused by convolution of laser pulse with receiving 
gating function. 
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£or the case when the bln is a factor of 2 less dense than bln i - 1, 



the ratio becomes 


2Aj ♦ dg 


(3.27) 


Wo see then that under the most unfavorable conditions where steep gradients 
occur in the sodium profile, returns could produce as much as a 30% error in 
some range bins. However when adjacent range bins are of comparable number 
density, the error should be minimal. In summary. If the density is 
increasing with altitude, then the measured gradient will be less than the 
true gradient. Alternatively if the sodium density is decreasing with 
altitude, the measured gradient will be greater than the true gradient. 

These effects would be enhanced if the range bin size were reduced without 
reducing the laser pulse width. Then deconvolution procedures would need 
to be introduced to solve (3.18). 

For completeness wo note that the minimxnu resolution of the system is 
determined by the lifetime, t, of the excited state. If we allow 4t seconds 
for the excited atoms to return to the ground state, which corresponds to 
92% of the excited atoms, then the minimum range resolution for a monostatic 


Ildar configuration is 


= 9,6m 


(3.28) 


3.3.3 ?itnpHfi&d ^sepvesasion fov r&timm fi*om a vanga Hn. Applying 
the approximations discussed in Section 3.3.2 to equation (3.13) we obtain 


IV, (V) 


4it[c^ + Aa/2)“ 


P(a) exp - [2T^(e)3da 


(3.2.9) 







Recauso the system receives data in the form of range bins having thickness 
As, it is natural to replace the continuous density profile, p(s), by a 
discrete profile, p^^., which represents the average density in the range 
bln, bet us compute il?, for the case £ « 0. If we expand t fa) as defined 
In (3.8) by assvming no variation of cr^ with altitude and consider the 
discrete profile mentioned above, we obtain 


<Vq(v) - 




4v[a,, + Aa/2]' 


pQ exp - (2pQO^(e - SQ))cfa 


0 (3.30) 

Wo may use the Taylor’s series expansion for the exponential factor as 
discussed in Section 3.3,2 and pei'forra the integration as indicated below 


f®l - - 

%J f'o ~ ^ ^0%^ 

(5.31) 

To eliminate the fi'equency dependence of A*q(v) we nmst Integrate over 
all V, Substituting the approximation (3,31) into (3.30) .and performing the 
intogi'ation I'esults In 


A ;r„(i 

» . ,v _J a ^ 

" ^ 4n(a . Aa/23*- 




(3.32) 


0 


This quantity represents the expected number of covints registered in the 
i *» 0 bin. 

To calcxilate the number of photons received from the t‘~ bit», we must 
account for the absorption of the laser line as it passes through bins, 

J, J < i, bet us assume that the laser linewidth is a Gaussian given by 





















where 


^(v) = C7 q exp - [w/a]' 


" ■ 2 if t'' ■ V 


C3.33) 



Av, 



ln2 

IT 


a = Av^/Av^ 


AVp, Av^ = absorption line Doppler width and laser linewidth 
respectively . 

2 

If we express the absorption coefficient as exp [-W ] , then 

the absorption incurred by the beam passing through the ith range bin is 
given by 


A = 


{1 - exp[-?:QZ expC-w^)}} exp[-0«>/a)^3<^w 


exp[-(u/a)^]dw 


(3.34) 


A series expansion exists for the quotient [Mitahell and Zemanakyj, 1971] 
and is given by 


A = 


V 




Vi + Ct^ 2 m + 2a ^ 


^ , 




nr 


,vr: 


+ . . . 


(3.35) 


no 


We may eliminate the higher order terms if the value of k^t is very 

small. To estimate its value we take t = 1,5 km which corresponds to the 

thickness of a range bin. The value of is given by a^p as indicated in 

equation (3.3). If we take the maximum cross section for the D- line of 
-‘16 2 

^Oa " * 10 m and a typical maximum density for sodium of 

9 -3 . -3 

p = 5 X 10 m we obtain kQl = 7,1 x 10 , Hence A may be approximated by 

the first term of the series. 
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If the laser linewidth is * 10 pm, then Av^. •* ^ AX^ ■ 8.6 « 10® 

Hz for X B 589 nm. The quantity Av^ was calculated in equation (2.10). Hence 
as (8.6 X 10® Hz)/ (1,07 x 10® Hz) ■ 8.06. Therefore 

kl 

A « (AX^ » 10 pm,ir = 200®X) (3.36) 

•th i 

To calculate the flux impinging on the i£^bln, we account for the 
absorption in each range bln ^ , d < and write 





(3.37) 




. , CTrtP • As 

■ fio ' 8.06 ^ 


(3.38) 


Often an upper bound on the estimate of the transmission through the 
layer is given by 


i-1 _ 

An estimate on the effect of this transmission on the sodium density calcu- 
lated is givei^ in Section 3.5.3. The final expression for the expected num- 
ber of counts received from the bin is given by 




y- 

4tr[s. + As/2) 


Iq % 9'(v) iv 


(3.40) 


where 2*^^ is given by either expression (3.38) or (3.39), 

3.4 Effective Scattering Cross Section, %ff 

3.4.1 Definition. The expression for the e^qpected number of photons re- 


ceived from the i~ range bin (3.40) contains an integral of the product of the 
laser line and the cross section of sodium. We may rewrite this integral as 


/; o„ dv . 


C3.41) 
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where n represents the total number of photons emitted in one laser pulse 
and is the effective scattering cross section. Hence is defined 
by the equation 


/q (v) dv 




(5.42) 


3.4.2 Appvoximate values for targe and small laser lineaidths. As a 
first approximation, let us assume that the laser pulse is square in the fre- 
quency domain with width Av^ « Av^. Employing the simple formula for the 
integrated cross section (see equation (2.13b)) 


^0 "v 


r-sL / 


(SI units) 


(3.43) 


and taking g-(v) to be constant over the width of the D, line we obtain 


* jf 

j = g f 
eff Av^ 


Av^ » AVpj Av < (3.44) 


The value of this expression is appropriate not only when the conditions 
specified in equation (3.44) are satisfied but also when the laser spectrum 
may cover different sections of the ZJj line on a shot to shot basis because 
of temperature variations of the laser tuning elements. In the latter case 
we assume that g(v) is effectively constant over the line when averaged 
over many shots . 


Notice that as the laser linewidth is narrowed so that it becomes much 
smaller than the resonance line, increases to a maximum 


eff 4e^yn^c Av^, 


AVp » Av^ 


(3.45) 


where again g(v) is uniform over its frequency distribution. 
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3.4,3 ® Gau88icn laser frequenay profile. To compute 

more realistically, we will assume that is a Gaussian profile given by 


= ilry? ('■ - VJ' 


(3.46) 


where = center frequency of the laser line 

= laser line halfwidth 

Since is taken to be the sum of two Gaussians , one centered about the 
and the other about the 1 ) 2 ^ line as described in Section 2.4.3, an analytical 
solution to (3.41) exists and is derived below. 


We may express as 


eff 


{' 

>CD 

{■ 


Oa 


exp - 




2a 


D 


+ Ogj, exp - 




2 


2a 


B 


} 


exp - 


2ito. 


(V - Vq) 


L 2a, 


2 n 


} 


dv 


(3,47) 


where 


Av, 


'D 


8 ln2 
2 


Av. 


= variance of Doppler broadened profile 


= variance of laser line 


£ 8 ln2 

Vp^, = center frequency of D 2 ^ and D 2 ^ line, respectively, 

Tne integral of (3,47) may be divided into two similar integrals and we indi- 
cate the solution to one of these integrals below. The portion of (3.47) 
which involves the ^® written as 
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(3.48) 


Changing variables by letting u ■ v - and defining - Vq, (5.48) 


may be rewritten as 


exp - 


-« L 


^ L J J 


du (3.49) 


By completing the square of the quantity enclosed by square brackets we ob- 
tain a Gaussian with variance 

2 2 2 

— y — ‘ - "y and mean - — » s- . Hence (3.49) may be 

V * V V ’ \ 

written as 


12 2 2 
a„ / 2vo_, 0- V 

Qq / D L a 


(3.50) 


A similar expression may, be obtained for the line contribution. Sub- 

stituting the appropriate values for a^, cr^ and we obtain the final ex- 
pression 


W /-“2“ 2 1 0^ 

Vav/ + Av/ t 


exp - 


4 ln2(v^ - v^)‘ 


* ®*P - 


4 macvp^ - v„)‘ 


(3. SI) 


Values of for various laser linewidths are plotted in Figure 3.7. Note 
that if the laser is tuned to the line rather than the center of the z ?2 
line, an increase in is obtained for small AX^ . 


j- j . . ..til-:- 5 
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We may now compare the value of obtained by the approximation (3.44) 
with the more exact expression (3. SI) for » Av^. Substituting the ex- 
pression for the scattering cross section of the line (2.14b) and using 
the relative weights of the and lines as given in (2.18), equation 
(3.51) becomes 

.2 


2 /in2 


. f r 20 


“I - V 


’'0 V 


Av. 


L 

12 P V‘ ll 




■ST 


(3.52) 


9 _-l 


The Doppler width of the line is Av^^ » 1,07 x lo s as given in equation 


(2.10), The linewidth of the laser is calculated by Av^ » oAAj^Xq* 


For AX, 


10 pm, Av^ = 8.65 X 10^ s 


If we assume the laser emits at line 


center, namely o • 1697340 ra”^, and take and to be positioned rela- 




tive to Vq as calculated in (2,17) then we obtain from equation (3,52) 


=1.82 X 10“^^ m^ 


(AX^ =10 pm) 


(3.53) 


This value may be compared to that by applying (3.44), In this case em- 
ploying equation (3.44) 


(j „ - s 2,05 X 10*^^ m^ 
eff 


(AXx 


(3.54) 


- 10 pm) 

Hence we find that the approximate formula (3.44) gives a value about 12% 
higher than the expression obtained by more realistically modeling the laser 
frequency distribution. 

3.4,4 Sensitivity of tinesvidth and stability. It is worthwhile 

to study the changes in for values of Av^ comparable to Av^ since this 
is the regime in which the laser may be operated in the future. In Figure 
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3.7 we showed the increase in 9^^^ as the spectral width of the laser was 


narrowed and noted an enhancement in ^or small AX^. if Vq • v^. For 


a more realistic presentation., one must realize that tuning elements have at 
least a 30% insertion loss. If we assume an etalon narrows the laser line- 
width by a factor of 10. then we should experience a 30% decrease in pulse 
energy for every factor of 10 decrease in linewidth. Hence in reality al- 
though 9^^^ would increase by roughly two orders of magnitude with the use 
of two etalons. the laser energy would drop simultaneously by about one order 
of magnitude resulting in only an order of magnitude increase in returns. 

From expression (3.51) we also note that is sensitive to the center 
line of the laser. This effect is displayed in Figure 3.8. The series of 
plots show that although 9^^^. may increase as the laser linewidth is nar- 
rowed. laser instability could cancel all the benefits of a narrow linewidth. 
By using both Figures 3.7 and 3,8. it was decided that a linewidth of 10 pm 
was optimum for attaining the highest value of given the frequency sta- 
bility of the laser. Certainly when the stability of the laser is increased, 
then a more narrow linewidth may be considered. 

3.5 System Calibration 

3,5.1 Normalization by Rayleigh scattering rebums. Referring to equa- 
tion (3.40') one may observe that it is not necessary to measure d , n. 


and N. individually to obtain a value for p .. Rather, the sodium returns 


may be normalized by returns from the 30-40 bm region, assuming returns in 
this region are consistent with Rayleigh scattering. This may be established 
if there is little evidence of aerosol • layers or high altitude clouds giving 


rise to Mie scattering. The signal strength of returns from this altitude 


will determine the produce N^A^^'t\. 
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Since Raylei^ scattering is frequency non-specific, we may dispense 
with the integral in (3.40) and write the expression for the number of counts 
received from a range bin in the 30-40 km region as 




(3.55) 


where 7^ is the atmosphere transmission to the altitude and is the 
average density of the atmosphere in the range bin of thickness Aa centered 
aroimd Sp. The quantity Op(v) is the classical Rayleigh backscatter cross 


section. 


Taking the ratio of equation (3.40) and (3. 55) and solving for p. we 


obtain 


— T~r‘ — ^ ~~ 7 ~ 

%ff ^ ^ 


(3.56) 


The ratio 7^/7 is discussed in Section 3.5.3. The absorption due to sodium 
which appears as 7^^ is calculated in Section 3.5.4. the number density, 
of the atmosphere range cell used for normalization is taken to be 

n 

= 3.69 X 10*^^ m ^ at Sjj = 30 km [ciBA^ 1972]. 

3.5.2 Bayleigh saatteving cross section. The laser is only weakly 
polarized and hence we may assume that the incident light on the atmospheric 
range bins is unpolarized (natural) light. The angular Rayleigh scattering 
cross section is then given by 


» o • f • i|»(e) 


(3.57) 


where 


a = classical Rayleigh cross section for isotropic molecules 
f = factor which accounts for the anisotropy of the polariza- 
bility of air molecules 
i(f(e) = angular scattering functifm 






so 


The classical Rayleigh cross section is given by 


where 


32.3 CVIV 


* refractive index of air at the altitude of interest 
X a wavelength of incident light 
a atmospheric number density 

Often the quantity enclosed by parentheses will be approximated by 


yielding the expression 


- 1 small) 


(3.59) 


8.3 


(3.60) 


However the approximation (3.59) under the assumptions of Rayleigh scattering 
is not more correct [Van De Hulst^ 1957] and we will use expression (3,58) 
for the classical Rayleigh cross section rather than (3,60). 

The factor f is given for unpolarized light by 


6 + 3A 
6 - 7A 


(3.61) 


where A is the depolarization factor (A ~ .035) [Penndorft 1957]. The scat- 
tering is not isotropic and assuming isotropic air molecules the angular 
scattering fvmction is given for unpolarized light as ItTohnaonj 1953] 


3 cos 0+1 _ p(cos6) 




where p(cos0) is the Rayleigh phase function and ^^(6) converts cr to a 


(3.62) 


steradian basis. In this ease the phase function is normalized to unity 
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and represents the conservative case of perfect scattering [Ch<xndx>aBekheaPt 
I960]; absorption is not treated in this approach, A more exact value for 
p(cos6) is given by {Femdovfa 1957] 

p{cose) = 0.7629 (1 + 0.932 cos^e) C3.63) 

and will be used for this work. 

3 

Therefore for backscatter 6 = w, and ^^Cir) = ^ [0.983). Note that 
/ = 1.06 and hence we may make the approximation 


= Ir (1.04) = 


[3. 64) 


Taking the above into consideration, the Rayleigh backscatter cross section 
then becomes 

.2 

( 3 . 65 ) 


,, 4.2 C«s-1)' 

O-C”) • — j 5^ 


X pr 


-1 


The units of o are cm sr . To obtain OoCir)* which is used in equation 

I* a 

(3.56), Oj.Cn’) must be multiplied by 4ir. The Rayleigh backscatter cross sec- 
tion then becomes 


- -4 — 

X 


(3.66) 


Hie index of refraction shows a variation with- density and temperature given 
by \Permdorfa 1957] 


vdiere 


(n„ - 1) = [n„ - 1) 


1 + at, 


1 + ot 


R 


% 


n 


( 3 . 67 ) 
+8 


= index of refraction for standard air where - 1) x 10 
is calculated to be 27714 for X = 589 nm 
a = 0.00367 

^R* ~ density of air where Rayleigh scattering occurs and 

of standard air, respectively 







V *s » temperature of Rayleigh scattering medium measured in *C 
and of standard air (15**), respectively, 

25 *3 

The temperature at 30 km is about -40“C and assuming Pg *« 2.5 x lo m we 
calculated with the aid of equations (3.65) and (5,66) that 

OgCTT) = 5.2 X 10“^^ m^ (3.68) 

Notice that Cd ( it) depends on selected definitions and ’’standard" adopted 
values for several parameters. As a result the value of Op(ii) adopted will 
differ between authors. 

3.5,3 Atmoaphevio tranamCeaion in the 30-00 km i*egion. Equation (3.56) 
may be simplified by noting = 3*, We determine this by computing the ef 
feet of Rayleigh scattering and absorption on the laser beam between the alti- 
tudes and 3^ + As/2. For computation purposes we set s^ = 30 km and 
+ As/2 = 90 km. 

Tlie extinction coefficient due to Rayleigh scattering is given by 
^90 km ^ 

P = o„p„(a)c?s = 2 X 10'^ (3.69) 

^ J 30 km ^ ^ 

where we assume that is constant with altitude and of magnitude = 

“32 

7 X 10 m‘‘. The values for p^(s) are those of the CIRA 1972 atmosphere. 

The Chappuis absorption bands of ozone span the sodium line. Ozone 

^25 2 

has a maximum ab'^orption cross section at 590 nm on the order of 5 x 10 m*" 

2”^ -2 

[Banka and Kookax'ts, 1973) . If we assume a column density of 2 x 10 m 
[ShtmasaH and 0^at<Ja, 1974], then the extinction coefficient resulting from 
ozone absorption is 

= (5 X 10’^^ m^)(2 X 10^^ m~^) ^ 0,01 (3.70) 

With the above information we may now relate f and by writing 

T = .Tg exp - [8^ ^ = 0.99 (3.71) 


.Ja.... 


ir 
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Hence it is reasonable to equate T and 7^. If we neglect the attenuation 
in the laser beam due to sodium absorption, then we may write (3.56) in its 
simplest form as 


- 

P . =» 77-^ 


[ 3 . + As/23‘ 






" 1 ? 


(3.72) 


3.5,4 Aaoumay of the method. The calibration method depends on de- 
manding little deviation in returns from those expected by Rayleigh scattering 
in the 30-40 km region. Hence, if enhancements occur in this region due to 
aerosol scattering, calibration may be difficult if not impossible. In addi- 
tion, the method assumes that we know the number density of the atmosphere 
in the 30-40 km region. Hence the calibration is only as accurate as the 
model atmosphere used; this accuracy is 10% typically. 

Furthermore we assume a temperature for the 30-40 km region to calculate 
the Rayleigh scattering cross section. A variation in this temperature of 
± 10"C results in a + 4% change in returns which is not serious. More impor- 
tantly, we find a serious variation in Rayleigh scattering cross section 
values in the literature. Because of the various definitions, values of 
differ by as much as a factor of 2 and one must be aware of this, prior to 
comparing data between groups. 

Finally we determine the error if absorption in the sodium layer is 
neglected. In equation (3.14) we estimated the optical depth for a one way 
path through the layer of » 0.09. Hence a two way path will have a 
transmission of 

exp - 2(0.09) = 0.84 (3.73) 

and the top of the layer should be corrected by up to 16%. However, because 
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the linewidth of the laser can be measured to only a factor of 2 or 3 and 
with the systematic errors discussed above, errors incurred by neglecting 
this attenuation are comparatively small. For a qualitative presentation 
therefore, equation (3.72) applies. 

3.6 Signal-to~Noi86 Caloulat-ion 

3.6,1 ffoiee oontrihutions , In Section 3.3 we obtained an expression 
for the expected number of signal photons received from the scattering volume. 
In this section we will investigate the effect of noise on returns and cal- 
culate the signal-to-noise ratio for the system. The noise considered con- 
sists of background noise, detector noise, logic circuit generated noise, 
signal-induced noise, and self noise of the signal. 

Background, or sky noise, is caused by airglow or heavenly bodies which 
may enter the field of view of the telescope. Scattered light from nearby 
cities would also fall into this category. Its effect is strong enough that 
sometimes it prevents accurate measurements from being taken. Detector noise, 
which includes dark noise, is present when photoelectrons are emitted from 
the photocathode of the photomultiplier tube and as a consequence may be 
minimized by cooling the photomultiplier tube. Data channel failures result 
in logic circuit generated noise. This noise is easily observable and soft- 
ware routines are used to identify the problem^ It is not a prime source of 
noise. In addition the laser and receiver have been separated by a large 
distance to minimize interference problems. Signal -induced noise occurs when 
noise is generated after the irradiation of the photomultiplier tube. Jones 
and Kent [1974) suggest that the electrons are excited into trap levels which 
subsequently require less thermal energy for their removal. It has been 
shown that this effect occurs for low light levels and may be present for up 
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to 200 vs after the signal is quenched \FettifeVt 1975] . Since the scattered 
laser intensity is of relatively long duration, the effect is cumulative. 

That is the i— bin will have signal-induced noise contributions from the 
■£ - 1, i - 2, etc., bins. Our returns from the sodium layer occur approxi- 
mately 600 us after the firing of the laser and even with the cumulative ef^ 
feet, it seems doubtful that these returns are affected. However, since we 
normalize our returns on the basis of scattering from heights around 30 km, 
it is possible that these returns could be affected. We will assume that 
signal -induced noise is not a serious problem for times greater than 200 ys 
after the firing of the laser pulse. To verify this assumption one would 
need to run tests on the C31034A photomultiplier tube performance as different 
photomultiplier tubes vary by orders of magnitude in their amount of signal- 
induced noise [Youngt 1976], Finally self noise, or shot noise, results 
from statistical variations in the received si^al; and to estimate the effect 
we will assume that the received photons obey Poisson statistics. Its effect 
will become evident below. 

Under the above assumptions, the contribution from noise to the number 
of counts registered is proportional to the time the receiver gate is open 
(2Aa/e) and is given by an expression quite similar to that for signal pho- 
tons, If we let N represent the number of noise coimts registered in a range 
bin, then 

TA» fn AX A 

L " 

where 

B. - background brightness at the receiver 

A 

n = receiver efficiency 

£2^ = solid angle of the receiver 
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AX ■ receiver bandwidth centered at wavelength X ■ o/v 
« receiver area 

hi » energy of incident photon 
Njj m detector dark count 
and the other symbols have their usual meaning. 

3.6*2 FomutaHon of the S/N ratio. Several expressions for the 
signal-to»noi$e ratio are used in practice and care must be taken to be 
aware of the particular one in use. For our purposes the received signal 
and noise are statistical in nature since we are in a photon counting mode. 
Hence an appropriate measure of the noise is the standard deviation of the 
received counts rather than the number If we assume Poisson statistics, 
this standard deviation is + W where W.. Is determined by (3.40). Hence 

V V 

the signal -to-noise ratio for a range bin on a per shot basis is given by 
the ratio of the mean number of received Signal photons to the standard 
deviation of the mean number of received photons, namely 


S/a 


h 


/a . '+ N 


C3.75) 


% 

Many authors will multiply ff by a factor of 2 to account for the uncertainty 
in the noise count, Ideally if the noise is observed for a long time, N 
can be accurately determined. 

We note that for x shots, the signal-to-noise ratio is given by 


a; ff . i/x 

S/a = L — = _L— (5.76) 

faTTT 


Hence in theory we can increase the signal-to-noise ratio indefinitely by 
integrating over a greater and greater number of shots. In practice one is 
reluctant to do this since over the time scales involved, the sodium profile 
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may change and hence only an average profile rather than an instantaneous 
profile is obtained. We believe that integration times must be kept shorter 
than 1/2 hour to study the time variation of the layer. In applying (3.76) 
one must be aware of the appropriate approximations ^ The signal -to>noise 
ratio may be approximated by S/S s Vx N. for the low noise condition. In 
this case a high repetition rate, low energy per pulse system might be 
suitable. However the system does not always operate under low noise con- 
ditions. For N » iV., S/s s /a? S. / /s and a much longer integration time 
will be required to achieve an appropriate S/S ratio. In summary then (3.76) 
implies that the sodium Ildar system should incorporate not only a high energy 
per pulse laser but also one with high average power. 

5.6.5 System odlouldtion for the Vrbccna tidccr. It is instructive to 
calculate the signal-to-noise ratio for our system. Values for several of 
the parameters used in this calculation are given in Table 3.2. The re- 
ceiver solid angle may be calculated using the telescope divergence half 
angle, , by 


% ~ 


TT(a^s)' 


ira 


E 


(3.77) 


Hence = 7 x 10“^ sr. The detector dark count is taken from the specifica- 
tion sheet for the RCA 31034A photomultiplier tube (see Figure 5.4). The 
dark count is a function Of temperature and for -20“C is Sj^ = 50 s”^. Alter- 
natively = 5 X lO'^/lO ws. We use a value for - lO”^ W m”^ sr”^ pm~^ 
which is an appropriate value for scatter from city lights [Stewart and 

Hop fields 1965]. With this information the background count, S^ may be 

4-1 

calculated using (3.74) and has a value of = 5.9 x 10 s or alterna- 
tively 0.59/10 ps. For this example the system is background noise 


L 
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Table 3.2 Parameters of the Urbana Sodium Laser Radar System 


Wavelength 

Transmitted bandwidth 

Laser beam divergence (half angle) 

Laser pulse energy 

Pulse duration 

Receiver gate time 

Receiver area 

Receiver bandwidth 

Receiver efficiency 

Receiver beam divergence (half angle) 

Separation of receiver and transmitter 

Repetition rate 


589 nm 
10 pm 

tty < 1 mrad 

17 

100 roJ = 3 X 10 photons 

4 us FWHM 

10 us (1.5 km range bin) 
0.114 

5 nm 
0.05 

<x^ a 3 mrad 
25 m 
0.5 s"^ 
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limited. Usually the experiment is done on very clear nights when is 

more than an order of magnitude lower. Under those conditions the system is 

self-noise limited. For further computation we will take N - 0.6. 

Substituting the appropriate values of Table 3,2 into equation (3.40) 

and neglecting the optical depth of the sodium layer we may calculate the 

mean number of photoelectron counts from the range bin corresponding to the 

— 9-3 

peak of the sodium layer taken to be at a = 90 km and p . = 5 x 10 m . If 
2 

we let T =0.2, then N. =5, Hence employing equation (3.75), the signal- 
to-noise ratio becomes S/N =2.1 on a per shot basis. Hence in 400 shots, 
which corresponds to a 15 min interval, we would expect a 20-fold increase 
in the signal-to-noise ratio, 

3.6.4 Confidenoe interval. More exactly we must assign a confidence 

interval to our results [Hogg and Craig t 1970). We consider the experiment 

as one in which we seek to determine the mean number of counts, p, returned 

from a range bin. Each shot produces a value, p., which is close to p and 

2 

is the mean of a distribution with variance a . Furthermore we assume that 
the process obeys a Poisson distribution, namely 


P(y * fc) » (3*78) 

where k is an integral number of counts. 

Under the assumption of low noise (W < if.), the variance for the dis- 
tribution is simply the square of the signal-to-noise ratio or the total 
number of received counts from the particular bin. Hence, if X is to be 
within 10% of p, we write the probability of this event as 


P(.9p < X < l.lp) 


l.lp 

, I 

k=.9u 


k -p 


l.lif. 

I 

k^,9N, 

't 


iN.y 


exp-N^ 

in 


(3.79) 




If we take S/tl • 2 per shot» then in 25 shots the S/N ^ 10. Using (3.79) « 


the probability of the event of being within 10% of 100 photoelectron counts 


is 0.7. Alternatively for 225 shots the S/N * 30. In this case we have 


a probability of 0.99 of being within 10% of the actual mean. 


3.7 Sodium Line Satiopation 


3.7.1 Statement of the problem. For short duration laser pulses and a 


low divergent system, saturation effects may become important. Although it 


will become apparent that saturation is not a problem for our system, we dis< 


cuss the topic as a reminder for future system inq^trovements or studies of 


other constituents such as potassium [Megie et al,» 1978] where saturation 


effects may need to be considered. 


Let and ^2 represent the population number densities of the ground 


state and excited state, respectively. The absorption coefficient, K> is 


proportional to the population difference and is related to the effective 


cross section by 


*=v ' - "2) V/ 


C3.80) 


Hence for equal populations of the two states 07^ 1 = ^^ 2 ) , the material becomes 
transparent to the incident radiation. This condition of no absorption is 


called saturation. Obviously a change in the absorption coefficient will occur 


before saturation, and it is useful to determine if any effects are noted in 


our experiment. 


3.7.2 Solution to the rate equations. We may write the rate equations 


for the excited state population as 


- W2 - 


C3.81) 


where 


IV, a rate of induced emission or absorption 

% 




k ;Vj»' I 'k I ■ ? nJ i -ifA- ’ ' 1 
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T s lifetime of the excited species 


The rate tf, may be written in terms of the incident intensity of the beam« 


-2 


where I has units W*m , Specifically 




(3.82) 


The quantity in the brackets is the incident photon flux in units per fre- 
quency imit per second per square meter. Hence 




Vf 

.2 


S(o^2)‘ 


(3.85) 


where 


- effective scattering cross section of sodium 


Nit) - photon time response of the laser beam 
T = transmission of atmosphere 


= half angle divergence of the laser beam 


z - altitude of range bin 
We may rewrite (3,81) as 






- ~ - (^2 * 




NityT 


vCa^s)' 


(3.84) 


Let us assume that the laser ^ulse is a square pulse of duration At,, 

£i 


given analytically as 


Nit) = 


■{ 


N, 


0 < t < At, 


At, 


(3.85) 


otherwise 

where N^ represents the total number of photons emitted in the laser pulse. 
The solution to (3,84) is given by 








»2(« 


[1 - - C^A)] Alt. 

*|^ - «xp - exp - 


* - ttr 


C3.86) 


t > At. 


where 


(i 9 *t number density in scattering volume 


T m 

8 a 


v(a^s) 


effW 


I=i. 2 

7 ^ ^s 


3.7.3 Saturation criterion* At this point we may obtain an estimate 
for the condition of saturation. The maximum number of atoms which could be 
excited by a laser beam is given by 


At. 


P y a py 

^max ^ T 


(3.87) 


where 


p = number density of sodium 

y = scattering volxime 

= laser pulse duration 
T “ lifetime of the excited state 

If each photon excited an atom« then under the most favorable conditions the 
total number of atoms excited would be given by 


#excited atoms = 


r 


%ff 

Tv(Oj.a)2 


( 3 . 88 ) 


where the quantity enclosed by parenthesis is the ratio of the maximum ef- 
fective area that the medium presents for scattering divided by area of the 
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laser beam at altitude a. Equating C3.87) and (3.88) we obtain a crude 
condition for saturation, namely 



7T St Ij 


To 


eff 


(3.89) 


Equation (3.89) gives us a restraint on the laser parameters. In particular 
we note that the divergence effect is quadratic. Using the appropriate 
values from Table 3.2 and taking t » 16 ns, we obtain the condition for 
a = 90 km that 




sat 


7.0 X 10^® 


(3.90) 


0 -jj 

For our system * 5.0 x lo** hence saturation is obviously not a 

problem. In fact (3.90) would represent the condition of a laser with para- 

-4 

meters 250 J pulse energy and 10 rad half angle divergence. 

Although the above gives an approximation for saturation effects, we 
may explicitly obtain an expression which will indicate the deviation from 
expected returns for high photon flux conditions. Since stimulated emission 
is emitted in the same direction as the laser beam, only spontaneous emission 
contributes to the signal received at the ground. Neglecting any absorption 
of this spontaneous emission or any stimulated emission caused by it, the 
range bin volume re-emits As ^ 9 (*)/t photons per second. Hence the 

received number of photoelectrons from this range bin is given by 

2As 

(3.91) 
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N 


Tr(o^a) As 
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where tz is the range resolution and t " 0 corresponds to the bottom of the 
range bin. This equation is accurate for « 2/isi/a, and in fact we may 
extend the integration to «> since N 2 Ct) falls off rapidly before t > 2bzlo, 
The integral in (3.91) becomes 


■At 
• n 


^ pT [1 - exp - (t/T)]d:fr 


0 


r ^ (1 - exp - (At^/x)] exp - [(t - 

8 


= + (t - t)[1 - exp - (At^/x)]} 

8 

Substituting (3.92) into (3.91), we may rewite as 

R 


N. 


J?Csat) 


1 + H 


1 + 


Lt^ 1 + i? 


1 - exp - 


fA*^(l + i?)1 


JJ 


where 


(3.92) 


(3.93) 


T = lifetime of excited atom 
At^ = laser pulse width 
R =2t/t^ 

N. = number of counts received under the condition of no saturation 

To observe the effect of beam divergence on saturation, we plot the ratio 

^i?(sat)'^^i figure 3.9 as a function of and see that saturation effects 

depend on At^. In general saturation effects will only become appreciable 
-4 

for Op < 10 rad, and At. < 300 ns. 
h b 

Although equation (3.93) was derived for atmospheric -applications, it 


may be used for determining saturation effects in the sodium cell monitor. 
These effects may be important if returns are to be normalized on the basis of 
of sodium cell scattering. 








4,1 "Bvin&vgles of Dye Laaeva 

4.1.1 Laaing tranaitiona. Tlie general properties of org^ic dye lasers 
have been reviewed in the literature [Snavely^ 1969; Schafer, 1973; Magyar, 
1974], This section will present the information necessary for understanding 
the specific problems common to flashlanip-pumped dye lasers . 

There are many ways in which the excited dye molecule may return to the 
ground state. Some of these paths are indicated in the energy level diagram 
of Figure 4,1; the electronic ground state of the molecule is S^, and 
and $2 represent the first and second excited singlet states. Transitions 
between singlet states are spin allowed and hence correspond to strong absorp- 
tion bands. The molar extinction coefficient, e, is a measure of the strength 
of an absorption band and is defined by the relation 

in ilJX) 

^ ~ m, ^ 

where 

N - concentration of dye molecules 
L - dye solution path length 
Iq, I - incident and transmitted light intensity 
The dye used for this experiment, Rhodamine 6G Perchlorate, has a maximum 
molar extinction coefficient of 11.2 x 10 1 mol“^ cm , 

The lasing process begins with the excitation of the dye molecule from 
the lowest levels of the ground state to the higher vibrational -rotatl C'iial 
levels of the state iA h transition) . The molecule decays in a non- 
radiative process to the lowest level of the state from which it returns 
to the higher levels of the state via stimulated emission (fl ->■ a transi- 


SINGLET 

STATES 


triplet 

STATES 



Lgure 


^ 1- tvnical dye molecule with 

Sr.urCstralgM Un^s) 

a^«nLSrs\« by the 

A->b and B+a respectively. 
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4.1.2 Loea prvoeaaea and ti*iptet formation* To facilitate lasing, 
losses from the state cannot be excessive. The term, critical inversion, 
is used to describe the concentration of molecules in the excited state 
required to produce lasing. Some of the loss processes for flashlan^^pumped 
dye lasers are described below. 

Spontaneous emission from the to state (B ->• a transition) is 
known as fluorescence. The time required to decay by this process is on the 
order of 1-5 ns. Also since the molecule decays non-radiatively prior to 
fluorescence, the fluorescence maximum is shifted to the longer wavelengths 
in the Anti-Stokes region. This shift is important because the unexcited 
dye will not absorb the spontaneous emission. In Figure 4.2 the absorption 
and fluorescence bands are drawn as practically continuous over a band of 
frequencies. Because the dye molecule has about 50 atoms and hence about 
150 normal vibrations, collisional and electrostatic perturbations caused 
by the surrounding solvent molecules broaden these lines, resulting in a 
quasicontinuum of lines superimposed on each electronic state. The result 
is a mirror image of the absorption band displaced toward lower wave 
numbers by a reflection at the wave number of the purely electronic 
transition. 

Intersystera c7:‘ossing is another process by which molecules may be 
removed from the lasing action. In this case a molecule in the state 
may relax to a lower lying triplet state, Tj. The transition is forbidden 
but may be collisionally assisted. If it is a radiative process, then it 
is called phosphorescence.. Since the time spent in the triplet state is 
relatively long, molecules which decay to this state are effectively 
trapped and removed from the lasinp process. 







To compound the problem* transitions between triplet states are spin 
allowed. Hence these transitions can form absorption bands which may over- 
lap the spontaneous singlet state fluorescence emission. Absorption between 
singlet states such as and pr<??ent but is not an important loss 

process* since the population 5.n the state is relatively small. To 
minimize the effects of the triplet state formation, a punqp source whose 
intensity rises rapidly with time is required. Plashlamp punning is an 
accepted dye excitation method. 

4.1.3 Solvent and quenahing agent. The processes discussed so far 

have been related to the moleculur structure of the dye molecule. However 

the environment of the molecule* i.e.* solvent* is very Important. For 

exan^le dimerization of the molecule may provide competing quenching 

processes. The dimers usually have a strong absorption band at shorter 

wavelengths and an additional weaker band at the long wavelength side of 

the monomer band. One way to minimize dimerization* is to use a less 

polar solvent such as ethanol. For safety reasons the solvent was chosen 

to be half water and half ethanol [Grant and Hawley ^ 1975] . The equilibrium 

between monomers and dimers shifts toward the latter with increasing dye 

concentration and decreasing temperature; and at room temperatures* a 

-4 

concentration of 5 x lo M Rhodamine 6 G may stop lasing action completely. 

An optimum concentration for the laser used in this study is 5 x 10"^ M. 

Quenching agents may be added to the dye solution to reduce criplet 

formation [filing et al, * 1970] . The quenching agent provides a lower 

lying triplet state than the dye molecule* facilitating a triplet -trip let 

decay which then effectively increases the intersystera crossing rate 

between the and Sq states. The condensed phase meterial* 1, 3, 5* 7- 

-3 

cyclooctatetraene CCGT] in concentrations of 5 x 10 . M has performed 
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quite adequately [BappaZoa^de et 1970], 

4.2 Laser Mainframe 

4.2.1 Layout* the laset and associated moni toning equipment were 
designed and built to monitor the sodium layer as well as to have sufficient 
versatility to investigate transmitter options and modifications. Among the 
factors which determined the physical layout were financial constraints, 
possible inclusion of the transmitter into a mobile lidar system and space 
limitations. Although some design choices were not optimal, the resulting 
system is adequate for the first basic scientific experiments and has 
allowed us to define precisely the areas of development needed to accomplish 
future system improvements. 

The flashlamp'-pumped dye laser is shown in Figure 4.3. The optical 
bench consists of a 2.5 m I-beam (1) supported by two containers of sand 
(2) which serve to isolate acoustically the I-beam from the floor. Two 
plates of 1/4” steel were bolted to the I-beam. The top plate was tapped 
with 1/4” - 20 holes and the bottom plate became a spacing washer having 


clearance holes aligned with the holes of the top plate. 

Effort was made to avoid mechanical vibrations from directly shaking 
the I-beam. As a result the flashlamp capacitor is mounted independently 
of the I-beam with its terminals straddling the I-beam. The flashlamp 
assembly is bolted to a frame which bridges the I-beam but does not touch 
it. The dye cell mount is clamped to the I-beam and does not touch the 
flashlamp assembly; all optics are similarly bolted to the I-beam. Tubing 
for the flashlamp cooling, dye flow, and vacuum systems is attached to the 
side of the I-beam for convenience. 

Sampling of the laser beam is done by a light-activated computer 
trigger pulse generator, spectrometer, relative power indicator and sodium 
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cell monitor. All of these items are placed on various tables surrounding 
the laser and beam splitters are used to direct the light into the appropriate 
directions. A pictorial display of the laser transmitter system is shown 
in Figure 4.4. 

4.2.2 Ett-t-pHoat occoity 

4. 2. 2.1 Geneval aonsideration. Most flashlamp-pumped dye lasers 
employ an elliptical cavity for focusing the excitation light onto the dye 
cell. The elliptical cavity is usually made from aluminum and the machined 
surface is polished mechanically or by hand. The cavity may be made on a 
lathe or milling machine by feeding the material at an angle to the cutting 
tool as indicated in Figure 4.5. 

As a general rule, the eccentricity of the resulting cavity should be as 
small as possible [/?5ss, 1969]. The ray tracing diagrams of Figure 4.6 
indicate the projection distortion that occurs for ellipsoids of high 
eccentricity. To minimize projection distortion, the dye cell diameter should 
be small compared to the semi-axis of the ellipse. Associated with this rule 
is the idea that the homogeneous illumination is limited to a small region 
which is dependent on the size of the source (Figure 4.7a]. Assuming a two- 
dimensional analysis, we see that the boundaries of the homogeneously 
illuminated area. A, are given by 


2 », , 

1 _ a - a _ 1 - e 

r “a + c 1+e 


and 




1 - e' 


1 + e 


(4.2a) 


(4.2b) 
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Figure 4.5 Production of an exactly elliptical cylinder 
on a lathe or milling machine; <(i is the angle 
between the axis of rotation and the feed 
direction. With long cylinders, the mirrors 
are made in two halves because the axis of 
the cutting tool would otherwise interfere 
[/?d5ss, 1969]. 






Projection distortion for an 
elliptical cavity. The major 
axis of each ellipse is the 
same and the source of rays 
is displaced equally. 
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Figure 4.7 (a) Extreme cases of distorted 
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(b) Shadowing in elliptical 
cavity mirrors [R^ss ^ 1969]. 
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eccentricity' of the ellipse 


We obtain the approximation 


A ^ - (1 - 

nr^ irr^ Cl + 


(4.3) 


Hence If e ~ 0.4, then the homogeneous llltanlnated area is about 30% of the 
source cross sectional area. 

Several other properties of elliptical cavities become evident. The 
length of the cavity should be comparable to or larger than the dimension 
of the seml^axis. Hence multiple reflections resulting from end reflections 
are avoided. Direct light becomes more important as the flashlamp is brought 
closer to the dye cell; this occurs for elliptical cavities having small 
eccentricities. However under these conditions the dye cell would shade part 
of the cavity wall, resulting in an uneven excitation of the dye (Figure 4,7b) 
Furthermore when flashlamps are brought as close as 1-2 cm from the dye cell, 
flashlamp shock waves may create gradients in the index of refraction of the 
dye; the resulting spurious deflection of laser light passing through the 
dye cell then degrades tuning performance. 

Several problems are often overcome by the use of multiple elliptical 
cavities. It is important to remember that the energy density produced in 
the dye cell region is dependent on the beam angle of the flashlamp. In any 
case the beam angle of all flashlamps adds up to a constant 2-ir. Hence on 
first principles, the efficiency of multiple elliptical cavities should 
decrease as the number of lamps is Increased; the maximum attainable laser 
output should remain the same. This however is not observed experimentally. 
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IVhat occurs is a larger homogeneous illumination areii for multiple elliptical 
cavities. Furthermore gradient effects in the dye become more symmetrical as 
the number of ellipses increases, resulting in a more stable tuning arrange- 
ment. As a compromise betxveen maintenance cost and improved laser output, most 
groups have chosen a double elliptical cavity. In the extreme case, a cavity 
having many elliptical components approaches the direct light situation and 
a helical or coaxial flashlamp would then provide greater symmetry. 

4. 2. 2. 2 Conjoining of the vacuum^, dye ftow^ and flashlamp oooling 
systems. The elliptical cavity is central to the operation of the laser] 
and as a consequence, it becomes the interaction area for a number of laser 
subsystems. We discuss these subsystems in this section but postpone the 
construction details of the elliptical cavity until Section 4. 2. 3.1. 

The flashlamps are an ablating version and as such must be constantly 
evacuated. The pump and vacuum system is extremely simple in design 
(Figure 4,8a). A Welch 1402B pump (140 1/min) is used to evacuate the 
system and dry nitrogen is leaked into the system through a needle valve. 
Pressure is monitored on the pvimp side of the flashlamps by a manometer; 
typical pressure for the lan^s is in the 4-8 torr region depending on the 
life of the flashlamps. A trap consisting of a flask is placed between the 
manometer and the flashlamps to prevent water from entering the vacuum pump 
or manometer in the unlikely event that mixing occurs between the water 
cooling and vacuum systems where they are coaxial . In addition the trap 
collects much of the quartz dust before it reaches the manometer or pump. 

The dye flow system consists of a magnetically coupled stainless steel 
pump (Micropump Corp. #12-60-303-374), S ym filter (Filterrite Corp. 

# U001A4UCE), dye cell and assorted Imperial 'Eastman Poly Flo tubing 
(Figure 4.8b), The emphasis of the design was to provide a fast dye flow. 



FILTER PUMP RESERVOIR 


FLOW RATE* 

1.5 LITERS/min 


RESERVOIR PUMP 



FLOW RATE: 
3.3 LlTERS/tnin 


Figure 4.8 Laser subsystems (a) vacuum system, (b) dye flow 
system, (c) flashlamp cooling system. 
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to reduce bubbles to a minimum, to build the system from components 
impervious to a wide selection of solvents and to have the option of jet 
stream operation should laser pumping be employed. The pump system has a 
pressure-activated relay which will turn the pump off when a break occurs 
in the dye system. If a fire results from a ruptured dye flow system, this 
feature prevents more of the combustible solvent from being added to the 
fire. Large volumes of dye are needed to minimize the loss in power due to 
photolysis of the dye, and hence the system employs a stainless steel tub 
which increases the system capacity to about 20 liters. The dye cell 
consists of two detachable quartz windows (15 mm diameter x 10 mm thickness) 
and a 10 mm ID X 280 mm pyrex tube. These windows are sealed with 0-rings 
as indicated in Figure 4.10b. Thick quartz windows were chosen since we 
observed etalon effects with thinner windows [E. S. Richter in Edwards, 

1977]. With our present window holder arrangement, approximately 6 cm 
of dye is not excited while 26 cm length is excited. Hence some loss is 
incurred in the dye cell. 

The flashlamp cooling system consists of an 80-liter reservoir of 
deionized water, magnetically-coupled pump to minimize bubble formation 
(Micropump #10-41-316), flashlamp water jacket, and assorted hosing (Figure 
4.8c). Since the conductivity of the cooling water increases with laser 
use, the water is deionized prior to each run by using a Barnestead deionizer 
IDS041 Mixed Resin Cartridge) , Since portions of the flashlamp electrodes are 
immersed in this cooling water, it is imperative that the resistivity of 
the water is kept high to prevent a short circuiting of the lamp discharge 
path by the water. 

4.2.3 Optical pump 

' i 

4. 2. 3.1 Flashlamp mount and assembly. The elliptical cavity assembly 


J 
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is patterned after that of QU>son [1972] and is shown in Figures 4.9 and 4.10, 
The double elliptical cavity was used originally in a ruby laser; and although 
it does not satis^ the criterion discussed in Section 4. 2. 2.1, it is adequate 
for first results. The eccentricity of each ellipse is 0.7 (major axis = 

75 mm, minor axis = 50 mm) and the surface is marred in many spots from flash- 
lamp breakage during operation. The flashlamps are held by lucite stqjport 
pieces and are allowed to move slightly because of the 0-ring seals. This 
freedom of motion is advantageous since the flashlan 5 )s undergo shock after 
each pulse. Lucite was chosen for the pieces which come in contact with water 
since lucite does not absorb water. All other pieces were fabricated from 
bakelite. Initially there was a problem with an arc which originated at the 
high voltage electrode and traveled down the cooling jacket vintil it com- 
pleted its path to the grounded elliptical cavity. This was overcome by 
placing teflon tape around the water jacket at the appropriate location, and 
replacing the high-voltage end of the cavity with a bakelite- teflon assembly. 

The flashlamps are a variation of the design of Ferror [1969], These 
ablating flashlamps generally radiate 15-20% of the electrical input energy 
into the 200-400 nm range and may be considered as 20,000-30, 000®K blackbody 
radiators for this wavelength range. The lamps consist of a 5 mm ID x lo mm 
OD X 280 mm quartz tube to which two brass electrodes are connected as shown 
in Figure 4,10 and 4.11. The electrodes are not symmetric because the walls 
of the elliptical cavity are of different thickness. 

Brass was chosen as a suitable electrode material because brass tubing 
was readily available, easy to machine, and did not sputter excessively. 

The electrodes will clog eventually with glass deposits in spite of the 
nitrogen atmosphere, but the lamps are replaced well before the constriction 
is con^jlete. A piece of metal tape (Scotch #5524) provides a mechanical 
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Figure 4.9 Photograph of the high voltage end of the elliptical 
cavity. 
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Figure 4,11 Flashlamp assembly. 
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seal between the quartz and brass and prevents a blast from exiting at this 
junction. Gum rubber tubing (5/16 IDj wetted with silicone grease (Doiv 
Corning #9 70V) is placed over this metal tape to provide a vacuum seal. 

Lamp assemblies last about 10,000 shots and are entirely replaced at that 
time; laser design insured that flashlaraps replacement did not upset the 
optical alignment. 

4. 2.3.2 TH.ggering and disehcage airauitry. The series combination 
of flashlamps is powered by a 19 kV, 0,5 A power supply. This poiver supply 
charges a 5yF/24kV low- inductance capacitor (Maxwell Laboratories, Series C, 
#33017) to 18.5-19 kV. The capacitor was chosen using the criterion of 
Hayworth [1971]. Once the capacitor is charged to the desired voltage, 
a pulser ionizes the gas in one of the lar^s causing both lamps to discharge 
the capacitor rapidly. 

The pulser incorporates an SCR circuit which drives an automobile spark 
coil (Figure 4.12a). Inputs to the pulser may be either from a timing 
circuit which generates a TTL level pulse every 1-5 sec (.variable) or from a 
manual input. The resulting 15 kV output pulse is passed through a SOOpF/ 
30kV D.C. blocking capacitor to the common junction of the flashlamps (Figure 
4.12b). Generally some adjustment in lamp pressure is required to obtain 
consistent triggering; but once this is accomplished, the laser will trigger 
faithfully throughout a run. 

4.3 Timing 

4.3.1 Mirror- grating conhination. In the early stages of the project 
a grating was used for crude tuning of the dye laser. It was found that the 
grating burned easily and ways were sought to minimize this. The present 
system employs a mirror- grating combination (MGC) [BgorkhoZm et al., 1971] 
and uses a high damage threshold grating (PTR optics #UTF-26, 1800 g/ram) . 
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Figure 4.12 (a) Spark gap initiator; (b) flashlamp electrical 

circuit . 
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The grating is mounted in auto coilimation (Littrow arrangement), and 
for this application the grating equation is expressed by 

n\ = 2d sin 0 C4.4) 

where 

n ~ order of diffraction 
X = wavelength 
d = groove spacing 
9 = angle of diffraction 

The incident beam having a center wavelength X and linewidth AX is diffracted 
into a multiplicity of discrete beams leaving the grating at different 
angles. Each of these beams is a spectral replica of the incident beam. 
Because the grating is blazed for the first order, most of the energy is 
concentrated in this order. The term echellette is sometimes used to describe 
a grating blazed for loitf orders (n = 1 or 2) - 

We chose a grating for coarse tuning because of its simplicity of use, 
availability, and relatively low cost. The grating was mounted such that 
all orders were parallel to the optical bench. This can be easily verified 
with a He-Ne laser by rotating the grating about an axis parallel to the laser 
beam until the laser beam is reflected back on itself for each order. Each 
order may be checked by rotating the grating about an axis perpendicular to 
the bench beam until the condition is met. Three strong orders are observed 
with our grating, corresponding to n = 1, 0, -1; the n = 1 order is easily 
resolved since the power in the « = 1 order is greater than the power in the 
r: = -1 order. The order n = 0 of course corresponds to specular reflection. 

Several problems exist with the use of gratings. Cleaning of gratings 
is virtually impossible; and in the dusty environment of the Aeronomy 
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Laboratory Field Station, if dust is not blotm off the grating with dry 
nitrogen, the small particles will absorb light and the resulting heat will 
destroy the reflective surface. In addition the reflectivity of the 
aluminized surface is about 90% and hence the surface is constantly absorbing 
heat; oxidation of the surface compounds the problem. Eventually the metal 
surface is evaporated and localized heating causes a warping of the plastic 
backing of the replica grating. Thus the groove configuration is eventually 
destroyed and tuning suffers. 

To minimize the grating incident energy, we adopted an MGC as shora 
schematically in Figure 4.13. The MGC may be considered as an etalon which 
has one mirror strongly reflective over a narrow wavelength range. Within 
this range the peak reflectivity of the MGC is given by 


R 


max 


^/R^ 

X * ^ 


C4.S) 


which is the expression for an etalon in antiresonance. As Figure 4.14 
indicates, the MGC can actually reduce cavity losses because the reflectivity 
of the MGC is higher than that of the grating. 

In addition the MGC reduces the incident power on the grating by a factor, 
n , where 


n 


1 - R. 


[1 + 


C4.6) 


For R. = 0.44 and R^ = 0.8, equations (4.S) and (4.6) imply that R =0.96 
and n = 22%. These values are appropriate for the present tuning configuration. 
It was found for i?^ > 0.5 that the laser v?ould lase off the partial 
reflector with the grating having no tuning effect. 


MIRROR GRATING 

OUTPUT COMBINATION 



Figure 4.13 Optical cavity presently used in the dye laser. 
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Figure 4.14 Curves of R^ax ® mirror- grating combination (MGC) 
when considered as an etalon in antiresonance. The 
reflectivity of the grating is denoted by R2 ivhile Rj 
represents the reflectivity of the mirror. 
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4.3.2 Etalon tuning. The use of etalons for tuning high energy dye 
lasers has become an accepted technique. An etalon is a Fabry-Perot inter- 
ferometer whose gap is held fixed by pressing the interferometer plates 
against a rigid spacer. For our application the medium between the plates is 
air; and in the ensuing calculations, the index of refraction for air will be 
assumed to be 1. 

The etalon may be considered as a filter which has maximum transmission 
at wavelengths satisfying 

nX = 2d cos e (4.7) 

n 

where 

n = an integer 

d - spacing between the plates 

B = angle of tilt of the etalon to the incident beam 

The integer n is termed the order of the interference; and because the 

optical wavelength is much smaller than the, physical separation of the plates, 

4 5 

n is usually on the order of 10 to 10 . 

Suppose plane leaves of equal intensity are incident over a range of 
angles, and the transmitted light is collected by a lens as shovm in Figure 
4.15. There will appear on the screen fringes of equal inclination in the 
focal plane of L, following loci of constant 0 when the order of interference 
has integral values, 

Minima will occur when n has half integral values. Considering two 
adjacent fringes for the same wavelength we write 

«Xj = 2d cos 


C4.8) 
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and 

(n = 2d cos 02 


C4.9) 


If we consider a spectral component X 2 close to Xj^ such that this component 
in order n falls on the order n - 1 of Xj^, then for AX = Xj - X 2 

2d cos 02 = - AXO (4.10) 


From equations (4.9) and (4,10) we find 
Xj = n AX 


(4.11) 


Substituting the value for n obtained from (4,8) into (4,11) vie obtain 


AX 



2d cos 0j 


(4.12) 


for small. This is the wavelength interval between fringes of adjacent 
orders for the same wavelength. It is the spacing between successive 
maxima of transmission and is called the free spectral range (F.S.R.) of 
the etalon. 

A measure of the sharpness of the fringes is given by the finesse, F, 
of the etalon. By definition the bandwidth of the fringe is given by 


B.W. 


F.S.R. 

F 


(4.13) 


and the relation between these quantities is shown diagrammatical ly in 
Figure 4.16. The finesse can be thought of as the effective number of 
interfering beams involved in forming the interference fringes. For our 
particular case the major factors which limit the net finesse are mirror 
reflectivity and lack of parallelism and/or planeness of the mirrors. The 
net finesse may be calculated as 



■z. 


(4.14) 


TRANSMISSION 



Figure 4.16 Transmission of an etalon as a function of 
wavelength and the relation between free 
spectral range and bandwidth. 
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where the represent the component finesses, 

The reflective finesse is given by the expression 


'R 1 - i? 


(4.15) 


where R = is the geometric mean of the reflectivity of the two plates. 

If M is the fractional deviation from planeness across the mirror aperture, 
or the measure of the lack of parallelism of the mirrors within the aperture 
being used, the flatness of figure finesse is given by 

Fp = M/2 (4.16) 

for a \/M plate. This plate flatness is normally specified at 500 nm. 

In the above discussion on bandwidth, the divergence of the laser has 
not been taken in account. If we consider only a plane parallel beam, the 
bandwidth of the lines transmitted by the etalon is given by 

B.W. = (4.17) 

However a contribution to the bandwidth also comes from the spread A0 in the 
angle of the laser beam. If we take the derivative of (4.8) and consider 
only magnitudes of the quantities involved, we, obtain 

= tan 0 ie (4.18) 

A 

For 6 small we may therefore write 

AX = X0A0 (4.19) 

To prevent light reflected from the etalon from feeding back through the gain 
medium, we require 0 > A0/2, Hence for a value of A0 = 5 mrad, the 
theoretical minimum value of the bandwidth is 7.4 pm. Realistically since 
the laser is nonlinearly sensitive to losses, the linewidth will vary with 
the pump level and may be much smaller than the computed value. 
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Several practical choices should be made when specifying plates. The 
overall finesse of the etalon should not be excessively high. For flashlamp- 
puraped dye lasers the finesse of each etalon should >s about 5, Standard 
interferometer plates have a flatness of X/20; and if the plates are as large 
as 1.5" in diameter j then over the area of the laser beam the flatness 
specification \flill be better than X/50 which is adequate. Plates should be 
wedged to avoid etalon effects from the outside surfaces, and these outside 
surfaces may be anti-reflection (A.R.) coated. Highly reflective plates will 
cause excessive walk-off losses because the plates are not perfectly flat. 
Reflectivities of 30% and 50% were found to be quite satisfactory, but 85% 
coatings were unsatisfactory unless the plates were exceptionally flat. 
Coatings should be able to withstand the energy of the laser beam; and since 
the etalon forms a resonant cavity, the energy density will be higher within 
the etalon than outside of it. We have found electron beam deposited 
dielectric coatings to be satisfactory but have burned some of the more 
traditional coatings. 

Etalon holders were fabricated from invar and are shoim in Figure 4.17. 
Spacers were made from either thin microscope slides or pieces of aluminum 
foil. For slightly larger gaps \ve understand that shafts of drills or 
bearing balls are quite acceptable. Temperature variation through a night 
required slight adjustments about every hour indicating an oven housing is 
perhaps ultimately needed. Day-to-day variations in temperature and/or 
pressure did cause sufficient drift that careful tuning was required prior 
to each run. 

It is important to remember that although these etalons are crude, they 
do work and cost perhaps 1/3 of a suitable commercial unit. One plate rests 
on a machined ridge while 3 screws are used to press the plates together. A 
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thin invar washer prevents direct contact between the screws and the plate. 

For satisfactory performance it is inqjortant that the coated surface of 
the plates be parallel. By illuminating the etalon with a He-Ne laser and 
looking at the output on a screen, a train of dots will be seen. These 
result from multiple reflections from the misaligned plates. By adjusting 
the three screws, one can collapse these dots to a single spot. 

For final alignment several methods exist. If one can illuminate the 
etalon with a large, collimated monochromatic beam, a couple of straight line 
fringes will be seen on the screen. Fine adjustment of the screws will 
yield even transmission across the aperture of the etalon indicating that the 
etalon is tuned to the wavelength of the incident beam. 

However most of our etalons were aligned using an alternate technique. 

In this case the emission from a sodium lamp was passed through a sheet of 
ground (frosted) glass. If one views the light source through the etalon 
with the eye focused at infinity, a system of rings will be seen with the 
reflected image of the pupil of the eye in the center of the ring pattern. 
^Vhen the eye is moved in any direction, the ring pattern will also move while 
the eye pupil remains at the center of the pattern. If for example the ring 
pattern expands in size as the eye is moved to the right, the plates then are 
further apart at the right than at the left. Suitable adjustments of the 
alignment screws will result in an aligned etalon which is characterized by 
observing that the rings remain the same size as the eye is moved to any 
point in the field of view, 

4.3.3 Monitoring the spectral emission 

4. 3. 3.1 Spectrometer . A 5" divided circle spectrometer is used for 
monitoring the tuning operation, and a used laser grating having a dispersion 
of 0.5 nra/mrad €500 nm is used as the dispersive element. Part of the laser 
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output and the focused output of an Oaram sodium lamp are superimposed onto 
the spectrometer slits. In this way the sodium emission appears as a 
doublet, and the laser may be tuned to the vicinity of the line at which 
point the sodium cell monitor must be used. 

The resolution of this spectrometer is such that the Z? lines appear to 
be spaced by about 2 mm. Hence since we know this spacing corresponds to 
0.6 nm, we judge the linewidth of the laser to be about 0.2 ran when only 
an MGC is used. With an etalon we estimate the linewidth to be 10 pm. 

Because of the resolution of the present equipment, we feel these estimates 
are accurate to a factor of 3. 

4. 3. 3. 2 Sodium vesonctnoe oeVL, The sodium cell is simply designed, 

yet adequate for first measurements. A more elaborate cell is desc ibed 

in Lory-Chanin [1965]. The cell used in this work consists of a glass 

envelope and allows for either absorption or scattering measurements. 

Because the laser linewidth is much greater than the linewidth, the latter 

measurement has been used exclusively. Two sodium dispensers are spot- 

welded to three electrodes and assembly-mounted in the base of the cell. 

The cell’s exhaust hole is then connected to a vacuum system and pumped down 

—6 —5 

for about 3 days resulting in a pressure in the 10~ - 10~ torr range; the 
vacuum is limited by the pump system used. Outgassing of the glass may be 
accomplished by heating the cell while under evacuation. Finally the cell 
is sealed. Sodium is released by passing a current of a few amperes through 
either dispenser. 

The cell is mounted in an oven built from fire bricks CFig^ire 4.18). A 
1500 watt heater coil is placed under each window and the input power to the 
heaters is controlled by a variac. With this configuration the windows become 
the hottest part of the cell and hence sodium is less likely to condense on 



Sodium cell and oven 
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the windows. Temperature is monitored by a mercury thermometer placed near 
the center of the cell and the equilibrixim temperature is kept to 160 ± 1®C. 
This temperature allows for bright scattering without excessive secondary 
scattering. 

It is instructive to calculate the temperature dependence of the 
scattering process to estimate the temperature control required. Since low 
pressures are involved, the line is mostly Doppler -broadened and equation 
(2.14a) applies, namely 



£ 

92 " 


(4.20) 


Employing the perfect gas law and taking 92^9i ~ ^ ot>tains 


kp = 1.57 X 10 


14 


T' 


,3/2 


(4.21) 


-1 

where has units m , P has units of atmospheres, and T has units “K. 

We note then that for 160 ± l^C, varies by less than 1%. 

Equation (4.21) may also be applied to determining the length of the 

cell for 1/e attenuation. In this case k^l = 1. If one knows the expected 

pressure for the evacuated cell, then by referring to the vapor pressure of 

sodium 1956] one may calculate 1. Generally if the cell can be 

—6 “5 

evacuated to 10~ - 10~ torr, significant attenuation will occur for 

T > ISO^C. Practically a sodium lamp source may be used to calibiate the 


cell . 
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5. RECEIVER DESCRIPTION AND SYSTEM OPERATION 
5 . 1 Reosiver Hcopdnare 

The initial receiver incoi^orated the Aeronoray Laboratory Field Station 
PDP-15 computer. The first observations of sodium over Urbana were made 
with this system; and during the early experimental phase, several areas 
for system improvement were defined. This first lidar receiver is described 
in Appendix I. The present receiving system is based on a Hewlett-Packard 
2114A computer and was designed for eventual inclusing in a mobile lidar 
system. In addition it affords more flexibility in a data collection. 

The present receiving system consists of three basic subsections: a 

telescope with photomultiplier tube, a lidar system interface with memory, 
and a HP2114A minicomputer (Figure 5.1). Although a complete description of 
this system may be found in Kinter [1977], a general discussion of the system 
hardware and operation follow below. 

The telescope is a 0.38 m (15 inch) diameter refractor which uses a plas- 
tic f/1.6 Fresnel lens as the objective. A major problem with a Fresnel lens 
is that collimated light is focused to a circle of least confusion which 
for this lens is 1 mm. This circle is larger than the diffraction limited 
spot size of a comparative perfect lens and hence the field of view for the 
telescope is larger than that of a telescope employing a more perfect objec- 
tive lens. Yet, because the lens is so large, the Fresnel lens is judged 
more cost effective. A subassembly which contains a field stop iris, a 
collimating lens, an aperture stop iris, and a 5 nm interference filter may 
be removed from the telescope for precise alignment under laboratory condi- 
tions. The schematic layout of the telescope is shown in Figure 5.2. 
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Figure 5.2 Telescope optical system. 
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In the photograph of Figure 5.3, the portion of the telescope (1) con- 
taining the subassembly can be clearly seen. A RCA C31034A/2A photomulti- 
plier tube (PMT) mounts in an assembly beneath the telescope (2) and is 
cooled to -20° C. This is accomplished by inserting a heating coil into a 
dewar of liquid nitrogen and circulating the resultant vapors around the 
PMT. A variac (4) is used to control the heat input and an ohmmeter C5) moni- 
tors the resistance of a calibrated thermistor mounted inside the PMT housing. 
The PMT views the received light through a double paned window to prevent 
frosting of the window. 

The C31034A/2A PMT has a quantum efficiency of 0.2 at 589 nm as shown 
in Figure 5.4. Unfortunately the tube has a characteristic of losing this 
high quantum efficiency with age. Presumably iirpurities emitted from the 
inner dynodes may eventually coat the gallium arsenide photocathode chip and 
thus reduce the efficiency. From the experience of other groups, this degra- 
dation in quantum efficiency is proportional tc use and can be accelerated 
if high anode currents frequently occur. Low light levels should be main- 
tained under all conditions of operation or storage. Finally the rectangular 
geometry of the chip is not a perfect match to the circular geometry of the 
telescope but problems do not result from this mismatch. 

Negative-going pulses from the photomultiplier tube are routed by a 
matched cable to an amplitude discriminator CUniversity of Illinois Elementary 
Particle Physics Research, Trigger TIOOA) where a threshold level of 50-100 mV 
is appropriate for good signal discrimination. The discriminator eliminates 
a large number of the pulses which originate from the inner dynodes of the 
PMT. The output of the discriminator is a NIM standard = OV into 50 £i, 
"1" = 0.8V into 50 n) and is fed through a matched line to the vicinity of 
the lidar interface unit where the signal is converted to TTL levels using 
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Figure 5.4 RCA C31034A photomultiplier characteristics. 
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a fast ECL circuit. The pulses are then presented to the interface unit. 

Referring to Figure 5.5, a brief summary of what transpires during 
data collection will now be given. The system interface includes a photon 
upcounter, processor, clock and random access memory (RAM) with the processor 
communicating with the computer and coordinating operation within the system 
interface. 

IVhen data are collected the following sequence of events takes place. 

An initialization pulse is given by the processor to prepare for data collec- 
tion. The lidar software allows the operator to load into the computer the 
desired sampling rate and delay, and these values are transferred by the 
processor from the computer to the interface. The sampling rate determines 
how frequently the binary upcounter will be sampled and therefore determines 
the size of the range bin. The sampling period then is programmable and may 
take on values from 1 ps to 128 ys. The delay corresponds to the time dif- 
ference between the firing of the laser and the initial sampling of the up- 
counter. The delay time is a multiple of the sampling period and may be from 
1 to 4095 times the sampling period. 

After this information is loaded, the system iifaits until the laser fires; 
the laser pulse triggers the system into operation. After the desired delay, 
the upcounter is sampled at the desired rate for 256 times. Each sample is 
routed to an addressed portion of the RA!*i. After this sequence the processor 
disables the clock signal and informs the computer that data collection has 
been completed. Between laser pulses the data are transferred from the RAM 
to the conq)uter. Adjacent bins are then subtracted to yield the number of 
counts registered in each bin. Tliis information is integrated with previous 
shots and the sequence continues until tlie desired number of shots is taken. 









m 


Since direct memory access was not available in the HP2114A computer, 
the inclusion of the device’s separate memory wss necessary, l^hile the maxi- 
mum available memory of 8k of core in the HP2114A is a severe limitation, 
the options of variable sampling rate and delay provide a more flexible sys- 
tem over the initial system described in Appendix I. Furthermore, the re- 
ceiving system may be easily included in a mobile lidar system. 

5.2 System Evaluation 

It is instructive to conpare the improvement in system performance over 
the duration of this work. Tivo profiles are superimposed in Figure 5.6 
which were taken about 1% years apart. As noted in Chapter 3, we define the 
maximum signal-to-noise ratio as the number of "signal" counts obtained from 
the peak of the layer divided by the standard deviation of the count from the 
peak of the layer. Explicitly this is written as 

S/N = — - (5.1) 

STTW 


where 

N. = "signal" count from the peak of the layer 

N = background noise count 

Referring to Figure 5.6 the signal-to-noise ratio is calculated to be 4.8 for 
the June data while by November 1977, S/N = 13. More importantly, it should 
be emphasized that the 1976 data required the integration of 1181 shots while 
the 1977 data was obtained with only 246 shots. Taking into account the dif- 
ference in integration time, the small change in range bin sizes, and the sea- 
sonal change in sodium density we find that the signal strength has been in- 
creased by a factor of IQ. 
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Figure 5.6 Con^arison of profiles taken years apart reflecting 
system improvements. Dotted curve was from the evening 
of June 25-26, 1976. Solid curve is taken from the evening 
of Nove,.ber 11-12, 1977. The November data are significant 
in themselves, since a primary peak centered at 90 km 
developed during the evening, indicating the variability of 
the layer. 
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5.3 System Alignment 

Alignment o£ receiver and transmitter is crucial to the proper operation 
of the laser radar system and a rather simple approach to the problem has 
been used. Our telescope, once leveled, is stationary. We therefore align 
the system by steering the laser beam to seek a peak in returns at some alti- 
tude range bin in the 20-30 km region. Generally 20 shots are sufficient to 
give statistically significant returns from these altitudes. Micrometer 
heads are used to steer the laser beam in either an east-west or a north- 
south direction. 

Figure 5.7 depicts the type of alignment returns we observe. In this 
particular case the final micrometer settings were taken relative to the 
arbitrary starting position as 10 mils east and 10 mils north. Alternatively, 
the angular movement from the starting position was approximately 5 rarad east 
and 7 mrad north. Finer adjustments of the micrometer heads would show a 
rounding of the curves near the peak. Improvements in the laser output have 
noii( allowed us to align at much higher altitudes than those indicated in 
Figure 5.7, although the procedure remains unchanged. 

Lengthy alignment is performed only when some major laser optical align- 
ment has occurred between runs. This alignment of transmitter and receive 
beams may take up to an hour to perform. However, once this alignment is 
achieved, the system requires only fine adjustments before future runs. This 
fine adjustment usually takes less than 15 minutes to accomplish. 
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6. EXPERIMENTAL RESULTS 

6. 1 Introduation 

The initial years of this project were spent designing, constructing, and 
optimizing the sodium lidar facility. During this time the techniques for 
successful lidar measurements of the sodium layer were developed, and 
monitoring of the layer concentrated on obtaining one or two profiles for an 
entire evning. By October 1977 the system signal-to-noise ratio had been 
increased substantially and experimental efforts focused on studying the 
short-term variability of the layer. Emphasis was placed on observing relative 
structural changes in the layer rather than measuring the absolute value of 
the sodium concentration versus altitude. In conjunction with this effort, 
signal processing techniques were developed to study the apparent wave motion 
often noticed in the layer. 

The theoretical tifork has indicated the need for several correlative 
measurements. Joint experiments with the Aeronomy Laboratory Field Station 
meteor radar facility and ionosonde are being planned. The signal processing 
techniques are briefly discussed in Section 6.2. Results from the nighttime 
monitoring of the sodium layer are presented in Section 6.3. Section 6.4 
presents the results of a feasibility study for joint meteor radar and laser 
radar measurements. Section 6.5 describes the interesting data obtained from 
monitoring the pre-sunrise sodium layer and simultaneously sounding the 
mesosphere ivith a 2.66 MHz transmitter. And finally a brief conclusion is 
given in Section 6.6. 

6.2 Sodium Lidar Data Collection and Analysis 

6.2.1 Introduction. Our lidar system has been developed to measure 
the sodium concentration profile over Urbana, Illinois. Because low signal 
levels are encountered, the receiving system operates in the photon counting 
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mode. During the past year J. R. Rowlett and C. Sr Gardner have been 
investigating various techniques for processing photocount data to obtain 
estimates of the sodiian profile. The usual procedure is to divide the time 
axis into disjoint intervals and to count the number of photons detected 
during each time interval. Using the lidar equation, the photocount in each 
interval or range bin can be related to the density and is then plotted 
versus time or distance as a histogram. 

The main advai.tage of the histogram approach is the simplicity of the 
photon counting electronics. However, the technique only yields the density 
at discrete points and since the measurement is actually an integrated 
density, the range resolution is limited by the receiver gate time. In 
addition photon counting does not fully utilize all the available information 
in the signal. Useful information contained in the arrival times of the 
detected photons is suppressed by the histogram data. If the arrival times 
are known, it should be possible to obtain a more accurate estimate of the 
sodium profile. Unfortunately, because the arrival time measurement requires 
sophisticated and expensive electronics, it is seldom used in lidar 
applications. 

Initially the signal processing work concentrated on the problem of 
estimating the sodium density from the histogram data using statistical 
techniques such as least-squares and minimum variance estimators. Particular 
attention was given to deconvolving the effects of the receiver gate time and 
laser pulse width. In addition it was necessary to determine the receiver 
gate time effect on the accuracy of the sodium density estimate. 

A gate time of zero corresponds to an arrival time measuring receiver. 
Since histogram data can be generated from knowledge of the photon arrival 
times but not vice versa, it is expected that the density estimates become 
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more accurate as the gate time is decreased. This work has led to the 
development of a promising signal processing technique based upon frequency 
domain filters which can be used to generate continuous density profiles from 
the histogram data. Furthermore it has been shown that under appropriate 
conditions, histogram data can yield density estimates whose accuracy is 
comparable to that obtained with arrival time data. A detailed presentation 
of this itfork is given by Rowlett and Gardner [1978], 

In this chapter data are presented sometimes as a histogram and sometimes 

til 

as filtered data. In the former case the counts registered in the i — range 

bin are plotted at the center of the range bin; unless otherwise noted, the 

range bin size is 1.5 km (10 ps] . Error bars are obtained by equating the 

square root of the received count to the standard deviation of the returns. 

The processed (filtered) data can easily be recognized because the filter 

cutoff frequency always appears with these data. A horizontal scale labeled 

2 

"counts'' implies that the data have not been corrected for the l/s dependence 
mentioned in Chapter 3. In most cases these data have been taken with a 
150 m (1 ys) range bin size; occasionally the data were taken with a 450 m 
(3 ys) range bin size, and these cases are clearly indicated. 

6.2.2 Signal processing. In Section 3.3.2 it was mentioned that the 
lidar equation has the form of a two-fold convolution of the receiver gate 
function, the transmitter pulse, and the unknoum sodium profile. The profile 
is "smeared" by the effects of the laser pulse shape and the receiver gate, 
and io is necessary to deconvolve these smearing effects to obtain an estimate 
of the sodium profile. However, the observed count also contains noise 
counts from background radiation and from "dark counts" generated by the 
photomultiplier tube in the absence of any signal. The presence of noise 
complicates the deconvolution process. 


nr 


Deconvolution is an inherently' ill-conditioned process, because the 
deconvolution filter is essentially an inverse filter to the original 
convolution. From the frequency domain vie\jpoint, the convolution process 
attenuates the high-frequency spectral components of profiles. The 
attenuated components will then be amplified in deconvolution, but the 
measured count will contain primarily noise at these frequencies. Amplifying 
the noise components produces an unstable solution. To reduce this 
instability we incorporate a priori infonnation about the solution. 

We can gain an intuitive idea of the effect of using a priori information 
by examining the deconvolution procedure from the frequency domain vie\tfpoint. 
In the frequency domain, the lidar equation C3.181 becomes 

(7.5) 

t* V 

2 

where D(o)) is the transform of p(r)/r . To recover an estimate p(r) of the 
profile, we should pass the data through a filter 

This filter has poles at w = — , k an integer, k ^ 0 from the zeros of 

(?((!)). However, if we assume that DCm) is bandlimited to oj < we can set 
S(w) = 0 for u) larger than the cutoff of PCoj). The physical meaning of 
assuming D(u)) is bandlimited to some w is that p(r») has no sinusoidal 

O IT 

component whose wavelength is less than ■s; — , where a - speed of light. We 
would choose the cutoff as small as possible to exclude as much noise 
whHe at the same time retaining as much information as possible. Figure 
6.1 shows the deconvolution filter applied to the data collected on October 
19, 1977 with three different filter cutoffs. As w is decreased, the 
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estimates become less noisy, but some infoTmation is also lost. Presumably, 
a middle ground can be reached where most of the information in the signal is 
retained while most of the noise is excluded. The present opinion is that a 
cutoff of C2.7 km)"^ is near optimum for these data but ultimately the exact 
value will vary from one data set to another. 

6.3 The Nighttime Sodium Layer 

6.3.1 The initial monitoring of the layer and system performance tests. 
The initial monitoring of the sodium layer in the spring of 1977 concentrated 
on obtaining one or two profiles for an entire night. The peak of the layer 
was consistently found to occur near 90 km. The layer was typically quite 
narrow and showed on occasion the characteristic bifurcation. Evidence of 
secondary peaks on the topside of the layer also existed. Several of these 
profiles are shown in Figures 6, 2-6.5, 

By October 1977 the system signal -to-noise had increased enough to 
enable us to obtain several profiles in a single evening. Such an example 
is shown in Figure 6.6. In this figure the Rayleigh scattering is also 
indicated to emphasize that lower altitude returns fall off exponentially 
and are not affected by Mie scattering. Furthermore this night of data 
indicated tne variability of the sodium layer. The profiles in the latter 
part of the evening show a slight decrease in the altitude of the peak of the 
layer and a lowering of the bottoms ide. The layer was observed to broaden 
through the night. Each of these profiles represents the integrated return 
from a 520 shot sequence and indicate the improvement in signal-to-noise over 
the Spring 1977 profiles. 

Several runs were hampered by haze or moonlight. Since returns from 
these evenings were usually noisy, time was spent on these night checking 
system performance, A standard test for a resonance scatter experiment is to 
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tune off the 589 nm resonance line and monitor the returns. The absence of 
returns in this case would indicate resonance scattering is the dominant 
process at 90 km. In Figure 6,7 the results of such a test are presented. 

In most of this work the emphasis will be placed on observing relative 
changes in the layer. However it is important to verify the value of the 
peak concentration of the layer. Using the data of February 22, 1978 and 
applying equation (3.72) to the returns, the sodium profile is shown in 
Figure 6.8. The concentration scale is accurate to a factor of 3 which results 
from our inability to measure the linewidth of the laser precisely. 

6.3.2 deployment of the layev struature with time — obaervatione of 
wavelike atumatures and vconationa of the topside soale height. Some of the 
more interesting lidar data we have collected was obtained during the late 
night and early morning hours of November 11 and 12, 1977, and the night of 
November 12. The data were processed using the unbiased estimate filter and 
plotted on a linear scale in Figures 6.9 and 6.10. The amplitude spectra of 
the raw data indicated that the signal usually drops into the noise at a 
spatial frequency somewhere between (1.8 km) ^ and (2.7 km) 

The raw data for each profile were obtained by firing the laser 250 
times over a period of approximately 10 minutes. An additional 5 minutes is 
required by the computer to output the data on punched paper tape. Samples 
of the November 11 and 12 data are replotted in Figures 6,11 and 6.12 (solid 
lines) along with curves representing plus and minus one standard deviation 
about the estimated profiles (dashed lines). These plots are typical of the 
remaining data and illustrate the relative accuracy of the estimated sodium 
profiles. 

An anaraolous feature of these returns is that the peak of the layer 
remained at 85 km for much of the November 11-12 monitoring. On the basis 
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Figure 6.7 Returns showing the effect of laser tuning. Each graph is an 

integrated return from 260 shots on October 21, 1977. Starting 
time for the returns are from left to right 2135, 2157, and 2209 
hrs. The center set of returns was taken when the laser was tuned 
to resonance scattering by sodium atoms . 
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Figure 6.9 Time evolution of the sodium lidar data obtained during the late evening 
of November 11, 1977 and the early morning of November 12, 1977. The 
receiver gate time was 1 ysec. 


NOVEMBER IH2, 1977 FILTER CUTOFF = (2.7km) 
















■ 'I 

r 


Locol Stort Time = 3206 CST 


Loco) Start Time = 2206 CST 



* I * 

*■ I * 

• I * 

% 

i: 


Filter Cutoff =(L0km)‘*’ 
Receiver Gate = l/&sec 


} 2 *.. »\\ 


Filter Cutoff =(2.7 km)"' 
Receiver Gate = lyttsec 


te.o 24.0 
COUNTS 


18,0 24.0 

COUNTS 


Figure 6.H Samples of the November 11, 1977 data of Figure 6.9 (solid line) plotted along with 
curves representing plus and minus one standard deviation about the estimated 
profile (dashed lines). 
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of these data one wonders if the seasonal lowering of the peak of the layer as 
observed by Gibson and Sandfovd [1971] is caused by a systematic lowering or 
rather is the result of averaging several profiles, many of which have 
substantial contributions well below the average peak of the layer. The data 
also indicate a relation between the topside scale height and the altitude of 
the peak of the layer. The profile of 2120 hrs has a topside scale height of 
about 12 km; by 0026 hrs the scale height is closer to 5 Ion and the peak of 
the layer iias changed from 83 km to 90 km. These data are in good agreement 
with the observations of Megie and Blcanont [1977] which show a relation of the 
topside scale height to the altitude of the peak of the layer. 

The November 12, 1977 data (Figure 6.10] provide an equally interesting 
aspect to the sodium layer. On this night the layer bifurcated, forming 
peaks at 83 km and 95 km. The peak at 95 kiu is then observed to move downward 
with time with a phase velocity of about 5 km/2 hr or 70 cm s . This is the 
same order of magnitude for the domward transport of ionization by the wind 
shear mechanism as discussed in Chapter 7. 

Another interesting feature of the time evolution of the sodium lidar 
data shoim in Figures 6.9 and 6.10 is the change in the height of the bottom- 
side throughout the night. In Figure 6.9 the height of the bottomside goes 
through a minimum in altitude at about 2324 hrs and a maximum near 2206 and 
0026 hrs. On the night of November 12, 1977 the same phenomenon occurs but 
now the maximum height of the bottomside is at 2132 hrs while the minimum 
height occurs at 2302. In both cases the period of this oscillation appears 
to be near 2 1/2 hours, ^itt et at. [1974] have suggested that the bottomside 
gradient of the sodium layer may undergo undulatory motions which result from 
dynamical processes and our observations appear to support their conclxision. 

A similar oscillation was noted by Thomas et al. [1976]. 
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The data of January 3, 1978 indicated a highly variable layer (Figure 
6,13). Wave structure seemed evident but a bulk movement as suggested in 
the November 12, 1977 data could not be found. It is possible that several 
waves are interacting in the region causing cancellation at one time and 
reinforcement at another. Indeed the amplitude spectra of these data indicate 
the possibility of many waves being simultaneously present (Figures 6.14 and 
6. IS). 

It is well knoim that the wave structure in this region may be quite 
complicated. Evidence for the reflection of waves over Urbana exists [Eesa^ 
1976] and standing waves are therefore possible. 'Such a standing wave in the 
lidar data would be typified by no vertical movement of the nulls and peaks 
of the wave structure but rather a change in amplitude vdth time at any set 
altitude. Using a notch filter to study only a portion of the spectrum of 
the returns, Rowlett and (kapdnev [1978] obseirved such a structure. The wave 
profile appears to return to its starting profile about every hour. The 
vertical wavelength of the wave is slightly greater than S km. 

6,3.3 Broadening of the layer with time. Evidence for broadening of the 
layer with time has already been shoiim in Figure 6.6. The monitoring of the 
layer during the night of March 16-17, 1978 indicated less structural varia- 
tion than on previous nights (Figure 6.16). However the layer did broaden 
during the night, starting at a full width at quarter maximum of 16.4 km 
(2231 hrs) and becoming 19.2 Ion wide at 0350 hrs. If the distribution topside 
scale height is taken to be the altitude increase from the peak of the layer 
at which the density is equal to 1/e of the peak density, then the distribution 
scale height increases from 8,5 km to 10.7 km in this time range for an 
average rate of change of 400 m hr'^. The peak of the layer which started 
at 90.3 km fell by 1.7 km by the end of the observation period. This 
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Figure 6.14 Averaged spectrum from the first half of January 3, 1978 data. 
Time period covers 2230-0100 hrs . 
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broadening of the layer has also been observed by Megie and Blamont [1977] 
and the change in scale height with time is in good agreement with their 
observat ions. 

6,4 Purpose and Pretitninary Tests for Joint Meteor Radar and Laser Radar 

Measurements 

The Aeronomy Laboratory Field Station meteor radar facility can be used 
to infer the neutral wind velocity in the 90 km region by measuring the radial 
Doppler shift resulting from specular reflection of the 41 MHz transmitter 
pulse from the ionized trails of meteors. Central to the meteor-ion theory 
described i Chapter 7 is the wind shear mechanism. It becomes important, 
therefore, to determine if the wave motion detected in the sodium layer 
corresponds to the motion of the ajnbient neutral atmosphere oi’ results from 
layers of ions caused by the wind shear mechanism. Since the mechanism is 
sensitive to the direction of the wind, and this direction can be sensed by 
the meteor radar experiment, it should be possible to lend experimental 
evidence to the corkscrew mechanism as a source of sodium. If the sodium is 
affected predominately by the ambient neutral atmospheric motion, then sodium 
may be an effective tracer for atmospheric motion at nighttime since the layer 
is under d)"iamical control at this time. The verification of this hypothesis 
could be made with the meteor radar facility. 

The joint experiment is not a trivial procedure because the laser and 
lidar data collection system reside in the meteor radar transmitter room. 

On the evening of October 19-20, 1977 wo attempted a run in \vhich the laser 
radar and meteor radar systems were operated alternately. This metliod was 
chosen to avoid possible meteor radar transmitter^ interference problems in 
the lidar receiver and to guarantee the laser operator’s safety since the 
laser is in the vicinity of some high voltage points in the meteor radar 
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transmitter. Lidar runs of 20 minutes duration ^^ere alteraated with 40 minute 
meteor radar runs; and during this 40 minute period, the ambient temperature 
increased sufficiently to require laser tuning adjustments prior to each 
lidar run. In the summer it is expected that the ambient temperature during 
meteor radar runs could reach as high as 32 “C. 

Unfortunately the meteor radar data were lost as a result of a system 
malfunction. However the attempted run was useful because it was discovered 
that temperature control of the laser tuning optics would be required for 
joint runs. In addition it was observed that if the sodium cell monitor was 
moved to a remote location, the operator might safely monitor the laser 
emission. If tuning became necessai-y, the meteor radar transmitter could then 
be turned off to perfom the adjustment. This arrangement would greatly 
extend the operating time of both systems. 

On January 10, 1978 such a simultaneous run was attempted, but we found 
that the lidar receiver system was severely affected by the operation of the 
meteor radar unit; the interference could be eliminated if the data collection 
system was removed from the transmitter room. In addition the meteor radar 
transmitter prevented the laser from pulsing properly. Much time was spent 
shielding properly the pulsing unit but the effort was not successful. Hence 
it was concluded a simultaneous run between the two systems would require 
that the laser be moved to a room shielded from the fields of the neteor 
radar transmitter. 

6.5 The Ppe-Sunrise Sodium Layer 

As mentioned in Chapter 1, the concentration of sodium is believed to be 
proportional to electron density and atomic oxygen. In particular the steep 
gradient of the bottomside of the sodium layer may coincide with the electron- 
density ledge which is kno\\rn to exist in the 80-85 km region. A simple test 


141 


of this theory would be to monitor simultaneously the bottoraside of the so- 
dium layer and the electron-density ledge. An appropriate experiment for 
performing the latter measurement would involve the ionosonde. Unfortunately 
the Aeronomy Laboratory Field Station ionosonde was not available for this 
study, and plans are now underway to reactivate the unit. In the future, 
incoherent scatter radar measurements of the electron-density profile will be 
conducted thus permitting detailed study of the electron- density ledge. 

Since a suitable e:q)eriment was not available for monitoring the election 
density at night, it was decided to probe the D and B regions with a partial 
reflection experiment operating at 2,66 MHz. Tlie partial reflection tech- 
nique, though, is a daytime experiment, while the lidar operates only at night- 
time. As a compromise it was decided to attempt a near simultaneous run at 
da\vn. Since the D region is not formed at this time, partial reflections 
would not occur but perhaps main reflections from ionized layers in the lower 
E region would be present. Furthermore the pre- sunrise measurements would 
allow an opportunity to observe any significant changes in the sodium layer, 
possibly due to rapid chemical changes in the sunlit atmosphere. 

Gough [1975] made the interesting suggestion that photoionization of 
sodium atoms can be a significant source of ionization in the upper D region 
during twilight. If this is true, then it follows that the sodium layer may 
be responsible for the initial buildup of ionization in the 80-90 km region 
during the sunrise transition period. 

The pre-sunrise sodium layer \tfas monitored first on January 21, 1978 

(Figure 6.17). The data indicate p significant notch in the profiles at 

about 85 km. On February 22, 1978 (Figure 6.18] the pre-sunrise layer was 
again monitored and this notch became prominent about 45 minutes before 

ground sunrise. At 23 minutes before ground sunrise the returns from the 
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Figure 6.18 Sunrise sodium layer. Time of profiles left to right are 042S-0436, 

0439-0510, 0454-0507, 0510-0527, 0523-0534, and 0537-0549 respectively. 
Ground sunrise is 0636 for this day. 
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2.66 MHz radar system indicated a strong echo from 87 1cm and a weaker echo 
from 93 km (Figure 6,19). These returns are related to the gradient in 
electron density and coincide with the negative gradient in the sodium 
returns below the notch and above the peak of the sodium layer in the last 
profile of Figure 6.18. 

If one assumes a screening height of 40 km and neglects the atmospher. 
refraction of sunlight, then the geometry for this sunrise condition is shown 

3 

in Figure 6.20. Taking the radius of the earth to be 6.37 x 10 km, one 
finds that for solar zenith angles of 97.5° and 96.4°, the altitudes of 95 km 
and 80 km respectively are first illuminated. Since the earth rotates at 
4 min deg~^, these heights are illuminated 30 minutes and 26 minutes 
respectively before ground sunrise. Sodium profiles on both occasions could 
not be taken after 45 minutes before sunrise. 

it the layer did not change from the time of the last profile of Figure 
6.18 to the time that the photograph in Figure 6.19 was taken (At ~ 25 min), 
then it is possible that ionization of the soditmi is contributing to the 
increase in the electron density in the upper D region during sunrise. In 
this case sodium would be a significant source of ionization. 

6.6 Conolnsion 

The experuaental results of this chapter indicate the highly variable 
nature of the sodium layer. Wavelike structures have been observed in the 
layer by the use of frequency domain filters which can be used to generate 
continuous density profiles from the histogram data. These waves have 
wavelengths of 3-15 km and may travel downward with a phase velocity of 
50-100 cm s”^. On one occasion there was evidence for the presence 
of a standing wave in the layer. During two successive evenings the bottom- 
side of the layer was observed o oscillate with a period of about 2 1/2 hours. 
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Figure 6.19 A-scope trace of returns from the 2.66 MHz 
transmitter taken on February 22, 1978 at 
0613 hours. The main reflection from 105 km 
is due to the E layer. Ihe reflection from 
87 Km was intermittent approximately 20 minutes 
prior to the taking of this photograph. Ground 
sunrise is at 0636 hours for this day. 
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Many observations of the layer are also in agrooinent with those of 
Megie and Blamont [1977], The layer is observed to broaden through the night 
and the topside scale height has been observed on one occasion to change on 
the average of 400 m hr”^. In addition it appears that as the altitude of the 
peak of the layer lowers, the topside scale height increases. 
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7. A METEOR ABLATION-CLUSTER ION SODIUM THEORY 
7 . 1 Introduction 

The experimental data presented in the previous chapter indicate the 
highly variable nature of the neutral sodium layer. From rocket mass spec- 
trometer measurements and incoherent scatter radar measurements, it is known 
that Na resides in approximately the same altitude region. Furthermore, 
charged metallic species are dynamically controlled by an interaction of the 
neutral wind system and the geomagnetic field. Tliis control evidences itself 
by the fomation of metallic layers and their subsequent downward motion. It 
is therefore possible that variations in Na^ are responsible for at least 
some of the perturbations in the Na layer. 

However the theory demands that the return of Na^ to Na be rapid and 
previous chemical schemes for sodium do not predict this (Appendix II) . But 
an important part of the sodium chemistry has been largely neglected. Clus- 
ter ion reactions have been shown to occur in the laboratory and evidence 
exists for their presence in the upper atmosphere. In this work it is shown 
that the cluster ion scheme may provide a fast return of Na^ to Na and thus 
remove a major obstacle in the theory which relates the dynamics of Na to 
those of Na. 

S<ction 7.2 will cover the meteor ablation source of sodium and the 
do^mward sweeping of ions. Diffusion time scales are discussed in Section 
7.3. A cluster ion scheme is developed in Section 7.4. In Section 7.5 a 
chemical equilibrium model is presented and particular attention is given 
to the seasonal variation of sodium and the lack of diurnal variation. The 
experimental results are reviewed in Section 7.6 in the light of the meteor 
ablation-cluster ion theory. Finally the chapter is concluded with a brief 
summary in Section 7.7. 


149 


7.2 The Soio’oe of Sodium and the Dowmand Sleeping of Sodium Ionization 

7.2.1 TheoHee for the souroe of sodium. The source of metallic ele- 
ments in the upper atmosphere is a subject of much controversy. Theories 
have proposed either a terrestrial origin, such as a particulate matter from 
evaporated sea spray [Aliens 1970], or an extraterrestrial source. In the 
latter case another contested debate concerns the form of this source. 

Donahue and Meier [1967] suggest that it results from a layer of dust con- 
taining aerosols with radii less than 10”^ cm having a residence time of 
several years. Hunten and WaZlaoe [1967] suggested that this layer of dust 
could also act as a sink for sodium. Usually any hypothesis about a common 
source for the atmospheric metallic species is interpreted in the light of 
the numerous studies of the sodium layer. Fioooo and Visconti [1973] discuss 
the seasonal variation of sodium in the light of the dust theory and present 
a numerical model of the sodium layer which involves both photochemistry and 
transport [Fioooo and Visconti^ 1974] . Tliey also indicate that the dust 
layer source may be either extraterrestrial or terrestrial in nature. 

Gadsden [1967, 1968, 1969, 1970] has outlined a theory which employs 
meteor ablation as a source of sodium. His model includes a dotvnward sweeping 
of ionization but depends largely on the shape of the meteoric deposition 
profile to predict the general features of the sodium layer. A sink for 
sodium is assumed below the ablation region in a portion of the atmosphere 
where water vapor exists. The total influx of meteoric material required to 
sustain the sodium layer is calculated and is found to be in agreement with 
similar values derived from different methods. Although it is thought that 
the deposition of sodium from meteor ablation is in the atomic form, the 
truth is not certain. The consequences of significant Na deposition result- 
ing from meteoric ablation are discussed in Section 7.6. Also an apparent 
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weakness of Gadsden’s method is its inability to explain the seasonal varia- 
tion of sodium. The impact of a cluster ion scheme on this process is dis- 
cussed in Section 7.5. 

Megie et at, [1978] have shed some new light on the source of sodium 
by their simultaneous measurements of sodium and potassium using the lidar 
technique. They find that the abundance ratio (Na/K) varies seasonally, 
being a factor of 4 to 5 higher during the winter; and they argue persua- 
sively for a dual source of sodium. A common meteoric source is taken to 
be constant over the year; and a terrestrial source resulting from a high 
latitude, winter vertical transport of salt particles, supplements the Na 
concentration. This conclusion is based on the assimiption that the chemical 
and ionic properties of sodium and potassium are very similar. However, the 
cluster ion scheme developed for this work may make this assumption invalid 
and this point is discussed further in Section 7.6. 

The most convincing arguments for the meteoric source are the clear 
correlations of metallic ion enhancement during meteor showers [NaJt'o-Lsi^ 

1968; Zhinden et at,, 1975] and of neutral sodium [Hake et at., 1972; Megie 
and Btamont, 1977], In addition the relative atmospheric abundances of 
metallic ions are roughly consistent with corresponding relative abundances 
in meteoric material [Bx>oi3n, 1973; Hevimann et at., 1977]. For these reasons 
and the sweeping mechanism described earlier, we will adopt meteor ablation 
as the source of sodium. 

7.2.2 Abtation of meteors. The process of meteor ablation as it applies 
to the sodium problem has been well presented [Gadsden, 1968, 1969]. The 
rate of deposition per unit height of meteoric material is a function of 
initial radius and velocity of the meteor. I'fhile atoms may be released by 
the process higher than 250 km, the main contribution is in the 80-120 km 
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region. In general meteors having velocities greater than 40 km and 
radii of a few microns ablate above 100 km [Gadsden^ 1967]. 

The number of meteors, n(r) , entering the earth's atmosphere per unit 
time per area is related to the radius, r by 

n( 2 »] « r~^dr C7.1) 

where the exponent is subject to some uncertainty. Since the ablation rate 

3 

is approximately proportional to r , it is obvious that small meteors are 
just as likely to contribute to the deposition as large ones. A typical 
entry velocity for many meteors is near 20 km s'^. The number of atoms per 
unit height interval ablated from a meteor having a starting radius of 100 
ym is plotted in Figure 7.1. Gadsden [1968] shows that this deposition is 
only weakly dependent on the velocity and initial direction of the meteor but 
that the height at which ablation occurs is strongly dependent upon the ini- 
tial velocity. In a later model calculation, Gadsden [1970] shows that 
meteor ablation can explain the small scale height on the topside of the 
layer if the gradient of the deposition profile is chosen correctly. Al- 
though Gadsden included the downward sweeping of ions in his calculation, he 
did not consider the layering of ions which can occur. The 95 km region is 
an area where sodium ions are dumped [Bcwe^ 1973; Smith and Medhtly^ 1972], 
hence the wind shear mechanism may also play a role in the control of the 
topside scale height, perhaps perturbing it as layers of ions move downward. 
Furthermore since several parameters can be adjusted in Gadsden's calculation 
to achieve a reasonable sodium layer profile, it is not clear if the deposi- 
tion profiles have been developed to agree solely with sodium. The downward 
sweeping of ions then provides the opportunity to collect Na^ in the 95 km 
region while keeping the process largely independent of the ablation profile. 
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7,2.3 l-ifetime. If sodium enters the atmosphere above 100 km, the 

eddy transport of the neutral species is greater than 5 days. However, the 


lifetime of Na against photoionization, 

Na + hv -> Na^ + J = 1.7 x 10 ^ s ^ (7,2) 

and daytime charge exchange 

Na + O 2 '*' Na'^ + O 2 fe = 6 x 10"^° cm^ s"^ (7.3) 

Na + NO"^ Na'^ + NO k = 8 x 10“^^ cra^ s"^ (7.4) 

is 16 hours and 10 hours respectively where it is assumed that [NO^] = [^2^^ ” 


4 -3 

4 X 10 cm . Hence it is clear that sodium is rapidly ionized during the 

daytime. Direct recombination of Na with electrons has a rate constant 
-12 3 -1 

fe = 10 cm s and is obviously too slow of a deionization process. Fur- 
thermore other reactions which convert Na to a shorter lived ion are slow 
because of the rarity of the atmosphere at this altitude. The only effective 
loss of Na^ at midlatitudes is thus the doivnward transport of ions by the 
corkscrew mechanism. 

7.2.4 Yerttoat movement of long lived ione: the oorksoreut meahanism. 

Dungey [1956] first described an interaction of the wind system and the 
geomagnetic field which resulted in the formation of ionized layers in the 
ionosphere. Since that time, the wind shear mechanism has been developed in 
quite some detail [e.g. lihiteheadt 1967; Axfovd, 1967] and has received ex- 
perimental support [Rowe^ 1973, 1974] . 

In addition to the formation of the layers, there is a net downward 
motion of these layers. This process has been termed the "corkscrew" mecha- 
nism because the effect is analogous to the pulling of a cork from a bottle 
by a corkscrew. Tliis drawing doim of ionization from the ionosphere applies 
to long lived ions, such as Na , and the process may dump these ions into 
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the lower E and D regions. The theory has been developed in Assford [1961, 
1963, 1967]. 

The predominant winds in the E region are those caused by the solar 
tides. In addition, perturbations in this large scale wind are caused by 
gravity waves which originate below the E region. The latter have been sug- 
gested as a source of sporadic-ff layers [Ckirnonas^ 1971], In both cases the 
doivnward directed phase velocity is thought to be responsible for the dom- 
ward motion of the ionization. If we neglect recombination and photoioniza- 
tion, we find that the corkscrew mechanism behaves differently in the E region 
than at lower altitudes. In both cases vertical movement of ions occurs, 
but different wind shears are involved. 

At the highest altitudes the ions are constrained to move along the 
geomagnetic lines of force because the collision frequency of ions with neu- 
trals is much less than their gyrofrequency. Hence the component of the wind 
parallel to the magnetic field drives the ions to form layers near the nulls 
of the wind profile. For midlatitudes it is the (geomagnetic) N-S winds 
which are responsible for this convergence. 

In the lower E region, calculations indicate that the vertical motion 
of the ions is dependent upon the (geomagnetic) E-W wind system [Axford 
and Cumdtda 1966], Since the collision frequency to gyrofrequency ratio for 
ions is large, positive ions move vertically because of the vertically induced 
electric fields which result from the E-W winds and the N-S component of the 
geomagnetic field. The ratio is still small for electrons, and therefore 
they continue to move along the geomagnetic lines of force. Since a vertical 
component to the magnetic field exists at midlatitudes, charge neutrality 


can be maintained. 
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7.2.5 Equation of motion. The motion of an ion is to a great extent 
affected by collisions with neutrals and interactions luith electric and mag- 
netic fields. The equation of motion for an ion is given by 


m. do, 
% 


^ X 5) + vm^(n - v^) - ~ 


(7.5) 


where 


m,, q, V. = mass, charge and velocity of the ion 
% = electric and magnetic field 

V = ion-neutral collision frequency 

= number density and partial pressure of the charged 
particle 

g = gravitational acceleration. 

It is apparent that the collision frequency term is a simplified expres- 

Sion. In this case otv(m - y^) is taken to be the force per unit volume 
exerted on the ions by collisions with the neutral gas. Collisions between 
ions may be neglected because of the large neutral to ion concentration ratio. 

A number of simplifications and approximations may be made. Tlie mag- 
netic field may be taken as being equal to the geomagnetic field and is 
written as 

^ = V ( 7 - 6 ) 

where is the field magnitude and r is a unit vector in the direction of 
the geomagnetic field. This is reasonable since at midlatitudes, ionospheric 
currents produce only minor perturbations in the magnetic declination 
[Cloutiev and Sandel, 1972], 
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If equation (7.5) is divided by then it becomes 


0) * dt 

tr 


E ^ 

^ + U . X f 


4- Cm - 5^) - „ s n 


Vp , + 2.— 


(7.7) 


where uj^ = is the ion gyrofrequency and is the ratio of the ion 

gyrofrequency to collision frequency. Since the gyrofrequency for Na^ is 
large (lo* = 230 s"^) the acceleration and gravity terms may be neglected 
[MacLeod et at, a 1975], 

By comparing the pressure term with v-fR, it is apparent that the former 
term may be neglected [Dungeya 1956] . This term becomes important at higher 
altitudes as temperature grows; or if the shape of the particular layer were 
a concern, then the pressure terra would be important. 

Finally by assuming that the neutral wind is a function of time and 
altitude, the vertical ion drift velocity may be obtained from the steady 
state solution of equation (7.5) 

q(^ + u. X B r) + m.v(M - V.) - 0. (7,8) 

'Z' U V 


To analyze this equation, consider the right handed coordinate system 
of Figure 7.2. The vertical coordinate is taken to be g \idth & and g 
pointing magnetically east and north respectively. Then i = where 

|s| = Bq- The neutral wind is assumed to have only horizc tal components 
and hence u = + y(3,i)5- Equation (7.8) may now be solved for the 

vertical component, w of the ion wind v.. By writing a similar eJipression 

'I' 

for electrons and equating the vertical velocity of ions to that of electrons 
so that charge neutrality is maintained, Ckimonas [1973] obtained the expres- 
sion for w at midlatitudes 


w 


—2 ^ [M(s,t) - R^V[3,t) 

1 + R. 

% 


(7.9) 
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where ff. = The quantity 5. increases with altitude as shown in 

Figure 7.3. 

Generally M(s,t) and have the form of tidal winds. Therefore, u 

(and analogously v) may be expressed as 

M(s,t) = exp[(s - 120/2g} sin k(.B - 120 - v^t) (7.10) 

where H is the local scale height and is the downward directed phase velo- 
city. Typical values for these constants for the midlatitude E region are 

E/q = 50 m s”^ 

H = 10 km 

k - 2ir/X = .2 km"^ 

= -1 m s”^ 

s - 120 = distance in km from 120 km 

yh - 

[ChimonaSi 1973; Axfords 1963]. 

To determine a numerical value for i?., values for and v. are needed. 

t' \J 'I 

For roidlatitudes, = 5.5 x 10*^2’ Gauss [Dudsiak et al.^ 1963]. Unfortunately 
the definition of the collision frequency, v, may take several forms [SUibbe^ 
1968] . One must keep in mind that each definition is a convenient quantity 
having dimension s . From the work of Banks [1966] the approximation for 
is given as 


v . 

1 - 


6.94 


n<7. 


^■n 


% 


[X] 


(7.11) 


where 

q = polarizability of the neutral particle 

= ionic charge 
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y. = m,m /{m, + mj) = reduced mass 
1 n' ^ % n*' 

m = ion and neutral particle mass 

*1^ ft 

[X] = neutral number density. 

Using the definition of Datgccmo [1961] = T<T/m^^ and the expression 

for the diffusion coefficient D. given by Chapman and CatDling [1970], the 

u 

collision frequency becomes 


V . = 6.94 ^ 




vn 


h 


[X] 


(7.12) 


Finally Stuhbe [1968] defines a "real" collision as one in which two parti- 
cles have direct contact. Under these conditions 

[X] (7.13) 


2 . 




%n 


For our purposes we will choose, somewhat arbitrarily, the value of 
V. given by equation (7.12). If the neutral atmosphere is taken to consist 


solely of O 2 and N 2 J then 


= 


1.89 X 10"^^[0,] + 1.9 X 10‘-"^[N„] 


- 12 , 


(7.14) 


-3 

where the number densities have units of cm . Implicit in this calculation 

is the use of the values for the polarizability of O 2 and N 2 given in Banka 

-24 3 -24 3 

[1966], i.e. Tiv, = 1.76 X 10 cm and = 1,59 x 10 cm . The varia- 
1^2 ^2 

tion of with altitude is plotted in Figure 7.3. 

The ion velocity may now be easily calculated as a function of time and 

altitude using equation (7.9) and (7,14). The profile for the ion drift 

velocity for t = 0 is shown in Figure 7.4 and compares favorably with the 

experimental results of Getter et at, [1975] . 
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From Figure 7,4 it is apparent that above 100 km large vertical ion 
drift velocities will cause the formation of ionized layers. This wind pro- 
file then moves downward, drawing the ionization along. Below 100 km the 
mechanism breaks down and the ions are "dumped" in this region. More appro- 
priately the ions are no longer able to continue traveling with the phase 
velocity of the wave but rather move downward with a velocity less than 
However a net doimward movement still exists as each passing wave lowers 
the ions left behind from the previous wave. 

7,3 Eddy and Molecular Diffusion 

* + 

7.3.1 Eispi!>e8s%ons . Since Na and Na are both distributed in narrow 
layers, there is always a tendency for these constituents to try to reduce 
the gradients of their boundaries. Molecular diffusion describes the trans- 
port of gas when the average molecular velocity of the minor constituent 
differs from that of the background gases; the process tends to separate the 
constituents according to theii mass. Eddy diffusion accounts for differences 
in composition and describes the motion of molecules due to turbrlence; this 
tends to mix the constituents and to restore the minor constituent's distri- 
bution height to the atmospheric scale height. Before chemical equilibrium 
can be assumed, the time constants for diffusion must be compared to the 
chemical time constants. 

It is a reasonable approximation to assume that sodium diffuses through 

t h 

a stationary N 2 atmosphere. The flux of tnis i — constituent may be expressed 
as [Colegvove et at*, 1965] 




fidr ^ 1 ■ 


\yr* 


n . 


(7.15) 
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where 


D 


H. 

V 


= molecular diffusion coefficient 
= krlm.g - scale height of the constituent 


rN, n. - mass and concentration of the minor constituent 
1 

g = gravitational acceleration 

k = Boltzmann s constant 

T = absolute temperature 

a = altitude. 

In an analogous manner the flux due to ddy diffusion may be iNfritten as 


^ -K 


dn . 

t ^ 

fl dT . l' 


[r Si"" «J 


1 


( 7 , 16 ) 


where 


K = eddy diffusion coefficient 
H = scale height of the background atmosphere. 

7.5.2 Diffusion coefficient values. If wc treat the molecules as per- 
fect spheres, then the molecular J fir >n coefficient may be calculated by 
using gas kinetic theory (r/wpmm.’ and i Ing^ 1970], 


D = 


5 

8fX]o^2 


* 


(7.17) 


where 

= (oj^ + a-,1/2 and Oj , arc the molecular diameters of the 
minor and major gases respectively 
k = Boltzmann’s constant 

T = absolute temperature 

= masses of the two tyjies of molecules 
[X] = molecular density of the main gas, taken to be N^. 
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Unfortunately the values for many of the molecular diameters are not Icnovm. 

A method for estimating the diameters of simple constituents, which involves 
knowledge of the molecular bond, has been described for a simple sodium 
chemistry 1973]. However it is known that Na and Na have collision 

diameters of 314 pm and 190 pm respectively. The altitude dependence of 
and is plotted in Figure 7.5. Molecules or ions containing sodium will 
have even lower molecular diffusion coefficients than those of Na and Na 
respectively. 

Values for the eddy diffusion coefficient are even more difficult to 
estimate. The problem is that all scales of motion are being parameterized 
by this coefficient. This includes not only local vertical mixinjj, but also 
larger structures such as meridional motions, large scale eddy fluxes, and 
tidal motions. Tlie coefficient also accounts for horizontal as well as verti- 
cal motion. Obviously there is no satisfactory way to parameterize these 
motions into a single coefficient and this problem is the main disadvantage 
of the eddy diffusion coefficient. In fact values are often picked to give 
the best results for a particular measurement. In the 90 km region, values 

vary by about two orders of magnitude depending on the measurement technique, 

7 2 

Analysis of chemical trails lead to K values slightly greater than 10 cm 

s"^ [Zitmeman and Tvowbridge^ 1973]. Studies of inert gases such as helium 

5 6 2 " 1 

or argon lead to much lower values between 10 and 10 cm s" . For this 

6 2 - 1 

work it will be assumed that X = 3 x lO cm s” for the sodium region which 
is the same order of magnitude as that deduced by Hunten [1975] and derived 
by Megie and Btamont [1977] through their analysis of the Na layer dynamics. 

7.3.3 Cornpayison of tijne soales, Tlie turbopause is thought to occur 
slightly above the sodium Ifi/er. Above the turbopause, molecular diffusion 
dominates vdiile below the turbopause eddy diffusion is most important. 
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Referring to Figure 7.5 it can now be seen that most likely eddy diffusion 
will be the more important form of diffusion in the region of the sodium layer. 
If a minor constituent is distributed with a concentration n = exp(-s/ff), 
then a diffusion time scale useful for comparing to chemical time constants 
is given by iP" /K. For = 3 km which is appropriate for the topside scale 
height of the sodium layer, the eddy diffusion time scale has a value of 
about 8 hours. It will be shoxm that most of the cluster ion reactions at 
95 km still have a reasonably fast time constant, on the order of 1 hour. 

The wind shear mechanism described in Section 7.2 will sweep ions dowward 
with an ion velocity of about 30 cm s^ [Gad^den^ 1970]. The time scale for 
this process between 110 km and 90 km will be 20 x 10^ cm/30 cm s “20 
hours. Since the topside gradient lifetime is estimated to be 2% times 
smaller than the time constant for the doivnward motion of Na , and Na can 
be converted to Na rapidly, the topside scale height should vary as layers 
of ions are brought doifli. The topside may show a diurnal variation. More 
importantly, since cluster ion reactions generally are three body reactions, 
their time constants will decrease with decreasing altitude. The lifetime 
of Na against photoionization and charge exchange is less than Jg of a day, 
hence chemical equilibrium is an appropriate assumption for the daytime. 
However at night the lifetime of Na increases substantially and dynamics 
should play a more important role. 

7.4 The Sodium Clustev Ion Chemistvy 

7.4.1 Introduation. The neutral and metal oxide ion chemistry is re- 
viewed in Appendix II. A chemical equilibrium calculation, which includes 
36 reactions, is performed and the calculation indicates that Na should 
dominate at ail altitudes. This result does not agree with experimental 
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evidence and implies that the model is not a realistic chemical scheme be- 
cause the return of Na to Na is too slow. 

In this section a new sodium chemistry based on cluster ion reactions 
is developed. The reactions are temperature sensitive and depend greatly on 
the atmospheric density since many of them are 3 body reactions. Both of 
these properties are investigated in Section 7.S. The scheme provides the 
fast chemical return of Na to Na necessary to relate the dynamics of Na to 
those of Na. 

7.4.2 Sodiwn aluster ion veaotione 

7.4. 2.1 Temperatupe dependence. Several measurements of sodium clus- 
ter ion reactions have been made and the results of these measurements ai*e 
shown in Table 7.1. Three pieces of information are lacking to make the 
table experimentally complete: first, kj, is not kno:m for the Na *^2 cluster 

reaction; second, the relative efficiencies of third bodies is not given; 
and third, tenperature relations for the rate coefficients are not tabulated. 
Referring to Figure 7.6 we see that the rates of Na^ cluster reactions should 
be roughly equal to those of NO^. Since we know kj> - 2 x lo^^ cm^ s^ at 
room temperature for the analogous reaction 

NO"^ 1^2 + N2 N0‘^*N2 + N2 C7.18) 

[Johmen et at., 1975] we will take this to be the value for in the Na^«N2 
reaction. Heimerl [1976] has tabulated the ratios of rate coefficients rela- 
tive to N2 for third bodies X = He, Ar, O2, and NO. In general these ratios 
/fey) are between .9 and 1.5 with the exception of He which exhibits 
rates almost 4 times slower than N2 {Heimerlj 1976] . Reactions involving 
H2O and CO2 are thought to proceed with the same rate as if N2 were the 
third body. 


Table 7.1 


Experimental data available for cluster ion reactions 


O Q 
^ S 

1 1 


Cluster Reaction 


Na + N 2 + N 2 Na 


[kcai niDl*^]* 
9.5 


Rate ** 

= K+1,-0.5) X 10 "^® cin^ 

esio'K 


= 2 X 10 
/c^ = 2 X 10“*^ 


Na"^ + COj + COj * Na**C02 + CO^ 


12.5 


ft = 2 X 10-29*“ 

7 -14^0 5 ®510-K 

= 1 X 10 


Na* +02+02^ Na**02 + 0^ 


5.7 


= 7 (+7,-4) X io‘2“ cm^ 
8310®K 


ftf = 5 X 10"^2(+1.S,-0.S) gjiQO^ 
ft^, . 8 X io-13(+1.5,-0.5) 


Source 

Iffilea et 1972] 
see text 

\Keller> emd Bei/eVj 1971] 

[Mlea et al*^ 1972] 
[Keller and Beyer^ 1971] 


Na* + tl20 + H 2 O ? Na*>H20 + HjO 24.6 

♦Calculated by Speara [1972] 

** Forward rate constant has units cm^ s”^. 


ft^= l(+.S.-.3) X 10-2® 

ft^ = 0 


Reverse rate constant has units cm^ s”*. 


[Johnson et 1971] 
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TOTAL EFFECTIVE POLARIZABILITY 
(IN UNITS OF 10'®® cm^) 
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Figure 7.6 Map of miniraum values for rate constants of the 
initial clustering reactions at about 300 K. 

The symbols denote measured values by the 
itforkers indicated [after iHtes et at., 19721 . 
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The temperature relations for the rate coefficients may be estimated in 
the following manner. As ^ ntioned in Appendix III the foward rate coeffi- 
cient is proportional to (7/ where b is between -1 and -2. We will there- 
fore take b = -1.5 for all the cluster reactions of Na. It is more difficult 
to estimate the temperature variation of the reverse reaction. The activa- 
tion term involves a change in enthalpy and therefore the exponential factor 
in equation CHI«7) may be written as 

expC-e/r) = expC+Aff/UT) (7.19) 

Bond energies for Na *02, Na *N 2 » Na *C02» and Na *H20 have been calculated 
by Spears [1972] and we may equate these values to the enthalpy change. 

Hence the temperature dependence of the reactions may now be established. 

As an example consider the reaction which produces Na^»N 2 (see Table 
7,1). The forward rate coefficient will be taken as 



X 10 


-31 


31 o' 

T 


1.5 


6 

cm 



(7.20) 


The change in enthalpy for this reaction is -AH =9.5 kcal/mole [Spears^ 
1972]. Applying equations (III. 7) and (7.19) the reverse rate coefficient 
may be written as 


, f-4750l ,, ,,, 

where the constant of proportionality is to be determined. The extra factor 

of T results from the change in units of the equilibrium constant, as dis- 

- 1 8 3 

cussed in Appendix III. From Table 7.1, = 10 cm at 310°K. Therefore 

the constant may be determined for 310° and the final expression for 
becomes 


A: = 9 X 10“^ 

r 



3 

cm s 


(7.22) 
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V. 


One may perform several checks on the above calculation. For example 
one sees from equations (III. 9), (III. 10), and (111.14) 

that the constant factors in the rate coefficients are related to the change 
in entropy by the relation 

~E ~ effC3i0"K) (7.23) 

Spears [1972] has calculated -2'AS = 7.2 kcal moi”^. Taking R = 1.99 cal 
(®K)~^ mol~^ and y = 310°K, 

^=-11.7 (7.24) 


3 -1 •"! 

For the right side of (7.23) we must take i? = 82 cm atm(**K)~ mol since 

3 23 -1 

a/c has units of cm . Using Avogadro's number of 6.02 x 10 mol , 


- In C2 X I0'^h(6.02 X IQ^^) _ 
(9 X io-‘7) (82) (310) " 


(7.25) 


Tliis is reasonable agreement; the errors result from the simultaneous use of 
the theoretical calculations of Spears [1972] and the experimental results 
of Niles et at, [1972]. 

The other rate coefficients are obtained in a similar manner and are 
summarized in Table 7,2. 

7. 4. 2. 2 Recombination, Recombination coefficients for the Na cluster 
ions have not been measured. In general they are in the range of lO"^ to 
10 cm s [Leu et al.^ 1973], and it is thought that at temperatures below 
300°K, the recombination rate coefficient displays a dependence; the 

dependence is predicted by theory and established in the laboratory [Bardsley^ 

1970]. Only at higher temperatures (~1000°K) will the dependence take on a 

— 1 5 +■ 

T ' law. Estimates for the recombination rate coefficient of NO »H20 have 

taken on values from 2 x io"“ to 3 x cm"^ s“ at 300°K [Arnold and 


Table 7.2 


Assumed rate coefficients for cluster ion reactions 


Reaction 

reference 

number* 

R, 


Reaction 


Na"^ + N 2 + X ^ Na'^<-N2 + X 


Rate Coefficients** 


K = 2 X 10 


-31 f310]^*^ 6 -1 


I ^ J 


cm s 


"-1 




2.5 


exp 


(^J 


cm 


R. 


Na"^ + CO 2 + X ^ Na'*’*C02 + X 


6 -1 
cm s 


h -2“ 10-^" (™) 


cm 


Na’^ + O 2 + X ^ Na"^-02 + X 


= 5 X 10 


.32 l^j 


1.5 


6 -1 
cm s 




= 8 X 10 - exp 




3 

cm s 


Na* + H 2 O + X ^ Na'^*H20 + X 


- X X X0-- 


1.5 


6 -1 
cm s 




* Reactions lacking a reference number are not used in the final scheme. 

**Positive subscripts denote forward reactions. Negative subscripts denote reverse reactions 
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Krarkomsky 3 1977; Reid, 1977; fhomae^ 1976], Although Keller and Beyer 

-I 3 

[1971] suggested a rate coefficient of 5 x lO"'' (y/SOO)" cm s" for recom- 
bination of Na"'**C 02 and Na ^*02 with electrons a more conservative view sug- 
gests that this be lowered. Therefore the recombination rate coefficient for 
all clusters will be taken to be a = 1 x lo”^ (r/300)~^*^ cm^ s“^. 

Recombination with negative ions can also take place and these rate coef- 

-7 3 -1 

ficients have values on the order of 10 cm s and temperature dependences 
of However since the model considers altitudes above 80 km, where nega- 

tive ion concentrations are thought to be small, dissociative recombination 
of the Na cluster ions with negative ions will not be considered. 

It is interesting to note that the reaction Na*»H20 + a” Na + H 2 O is 

2 

probably sufficiently exothermic to excite Na into its P state [Suddert 
1970; Sipler and Biondia 1978]. Therefore the reaction scheme may provide 
an additional source of sodium nightglow and in fact one may expect a corre- 
lation between the dotmward motion of Na^ layers and sodium emissions. In 
this theory the downward motion of Na* layers may also produce an apparent 
dovmward motion of Na as Na is converted to Na. Hence a correlation between 
do»«iwaxd moving Na layers and sodium emissions is predicted for altitudes 
above 90 km where the corkscrew mechanism is most evident. 

7. 4. 2. 3 Snitching. In previous studies switching reactions for 
the Na cluster ion scheme have not been include i and as a lesult the cluster 
ion scheme was thought to be untenable [Brown, 1973; Fergusonj 1972]. A 
switching reaction is a txtfo body reaction in which the weakly bound molecule 

is replaced by a more tightly bouui molecule. On the basis of laboratory 

“9 3-1 

experience the switching rates have values on the order of 10 cm s and 
do not exhibit a great temperature dependence [Ni-les et al .3 1972]. Switching 
with CO 2 is more effective than with H 2 O because of the higher concentration 
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of CO 2 in the atmosphere. Ultimately the schane should produce Na *H20 which 
is the most stable cluster ion of the four considered (Na *02, Na *N 2 , Na ‘002, 
and Na'^'Il 20 ). 

7.4. 2.4 PhotodsstrucHon of aluster ions> Photodissociation of Na^ 
clusters represents another area where the experimental evidence is sparse; 
this process could be another source of Na . It has been found that 

2 ’ °2^’®2» photodissociation cross sections 

in the 676.4 to 454.5 nra range; although, H’*'*(H20)^ =1-4 
have small or zero cross sections in this range [Beyei> and Vcmderhoffa 1976]. 
The cluster ion N0‘^*H20 was added to the list of ions with a sizeable cross 
section by Vanderkoff and Wiles [1976] and on the basis of this information 

4 

Petevson et at. [1977] suggest sunlight photodestruction rates for O 2 *H20 
and N0‘^*H20 of 0.6 and 0.06 s“^ respectively. However Na"^ is more tightly 
bound than NO*" in the cluster ion. For example the binding energy of N0"^*N2 
is about .2 eV which corresponds to an enthalpy change of 4.9 kcal/mol 
[2^i^er and Coway 3 1975]. However the calculated enthalpy change for Na'*"*N 2 
is 9.5 kcal/mol [Speans^ 1972]. On this basis it is likely that photodes- 
truction of the Na'*’ clusters is negligible, although it must be remembered 
that laboratory measurements have not been made far into the ultraviolet and 
hence the uncertainly remains. 

4 4 

It is probable that the photodestruction rate of Na *02, Na •N 2 , and 
Na^*C 02 will be much greater than that of Na^*H 20 since the latter is a 
more stable cluster. In fact Na^*H20 may not photodissociate at all. Since 
it is likely that the collisional breakup rate of the weaker bound ions is 
comparable to or greater than the photodestruction rate, it will be assumed 
that photodestruction for all of the cluster ion is not an important process. 


7.4.3 !Tke eodiim dhemtoal sink, A likely candidate for the sodium 
sink mentioned in Section 7.2 is the production of higher order cluster ions 
from Na *H20. The cluster ions are known to form a ledge between 80 and 85 
km and therefore conveniently approximate a perfect sink. A test of this 
hypothesis though, must include an estimate of the magnitude of the sink and 
a comparison of this value to the magnitude of the source from meteor ablation. 
If we assume that the cluster ion scheme repeats itself for higher order 
clusters as shown in Figure 7.7, then the controlling reaction will be 

Na'^*H20 + N 2 + X ->■ Na^-H20(N2) + X (Rll) C7.26) 

where the rate must be determined. Dsidic and Kehante [1970] obtained the 
enthalpy and Gibb’s free energy for the reactions 

Na"^ + H 2 O + H 2 O Na"^'H20 + H 2 O [7.27] 


and 


Na'^‘H20 + H 2 O + H 2 O Na'^‘CH20)2 + H 2 O (7.28) 

From the ratios of the equilibriiun constants and activation energy factors 

of the reverse rate coefficients for reactions (7.27) and (7.28), it is 

clear that the higher order clustering foarward rate coefficient is slower. 

Assuming this also holds for the clustering of N 2 to Na'*’*H 20 , k^^ may be 

-1 

conservatively set as 10 times the forward rate coefficient of the reaction 
Na'*' + N 2 + N 2 ^ Na‘‘'*N 2 + N 2 (Rl) (7.29) 

The value for may be derived eventually by aeronomic techniques, and 
this point is discussed in Section 7.5,1. 

Recombination of these higher order cluster will be neglected since [e~\ 
is small in the 80-^85 km region. The breakup of Na^*H20(N2) should not 
affect the estimate severely since the switching reaction R13 is undoubtedly 
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No+CHgOg 


Figure 7.7 Suggested chemical map for a sodium 
sink in the 80-85 km region. 

Na"^' (H 20)2 is then lost to higher 
order water clusters and possibly 
aerosols . 
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faster. Hence the forward rate of Rll should give a reasonable estimate to 
the rate of the loss process. 

Zhinden et at, [1975] measured a Na^*H20 layer centered at 80 km with a 

-3 

width of 6 km. The concentration at 80 km was about 10 cm . Since their 
measurement was taken during the Geminids meteor shower, values of the ions 
are certainly enhanced. Since the maximum Na concentration was about 5 times 

_3 

the peak value normally quoted (100 cm ), it is reasonable to assume that a 

+ -3 

nominal value for [Na ’H20] is 2 cm . This value would also account for the 

sporadic observations of Na *H20. Using the values of [N 2 ] and [X] at 80 km 

from Table II. 2, the column loss rate of Na* is about kj^j[Na**H20] [N 2 ] 

5 3-2-1 

[X]*6 X 10 cm = 4.8 X 10 cm s . If we take the total influx of meteoric 
9 -1 *“1 

material to be 10 gm yr (earth surface) \Zhivuien et at., 1975], then 

3 -2 -1 

this corresponds to a total influx of 1.5 x 10 cm s , in good agreement 
with the column loss. In the steady state the velocity for the downward mo- 

* 4 * * 1 * •!* 

tion of Na is given by the ratio of the influx of Na to the Na density, 
3-2-1 2 -3 -1 

or 1.5 X 10 cm s /I x 10 cm = 15 cm s which compares with Gadsden's 

_1 

[1970] suggestion of 30 cm s . The timescale for eddy diffusion in the 80-90 

km region is given by H^/K = (6 x 10^ cm)^/3 x 10^ cm^ s~^ = 1.2 x lo^ s 

where H is taken to be the scale height of the atmosphere. Hence the velo- 

. 5 -1 

city needed to move an ion from 90 km to 80 km in 1.2 x 10 s is 10 cm s 

7.4.4 The life of a sodiym atom. The picture that emerges is the fol- 
lowing. Ions are layered regularly above 95 km by the semidiurnal tides. 

These relatively high concentration packets of ions are then brought do(«i- 
ward with a velocity of about 1 m s ^ [Chimonas^ 1973] and have the charac- 
teristic that the average flux travels with velocity of 15-30 cm s~^. The 
layers may perturb the Na topside gradient in the process. Once the ions are 
dumped into the 90-95 km region, they diffuse downward to the sink at 80-85 km 
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by the process of eddy diffusion. This sink is most likely total. That iSj 
Na* is rapidly converted to higher order clusters, which may attach them- 
selves to aerosols, and are then removed from the sodium chemistry system by 
a process of sedimentation. Tlie sodium layer takes a small role in this 
process, acting somewhat like a reservoir. It provides Na^ by photoioniza- 
tion and charge exchange processes and receives Na from the dissociative 
recombination of Na^*H^O with electrons but at all times remains in chemical 
equilibri'om with NaO. 

7.4.5 Adopted veoatton scheme. The results of Sections 7.4. 1-7. 4. 4 
may now be combined to form the adopted reaction scheme. As discussed in 
Appendix II the neutral chemistry will be approximated by the interplay of 
Na and NaO. From the results of Section II. 2 the concept of sodium ion 
oxides playing an important role in the sodium chemistry will be rejected. 
Rather a cluster ion scheme will be used on the basis of discussions in 
Sections 7. 5. 1-7, 4. 4. The clustering of 0^ and H 2 O with Na^ may be neglected 
since the Na^*N 2 route is much more efficient in producing water cluster ions. 
Photodestruction of Na *N 2 and Na *C02 will not be considered since the colli- 
sional breakup of these clusters is probably of greater importance. The 
cluster Na^*H20 will be assumed to be too stable to photodestruct . Loss of 
Na in the form of Na^»H20 will be taken as having a rate coefficient equal 

to = ky/10. The details of the sink need not be specified for this calcu- 
lation and are therefore not included. Tlie final reaction scheme is shown 
in Figure 7.8 and the list of equations and rates are given in Table 7.3. 

7 . 4.6 Relating experimental Na and Na profiles as a test of the 
cluster ion reaction scheme. Although several measurements of neutral sodium 
by both ground-based and rocket-borne photometry techniques have been made 
over the past few decades, near simultaneous measurements of Na and Na have 
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SINK 


Figure 7.8 The adopted chemical scheme for 
the meteor ablation-cluster ion 
theory o£ sodium, A possible 
chemical structure for the sink 
is sho\m in Figure 7.7. Specific 
reactions are listed in Table 7.3. 







Table 7.3 


Reactions and rate coefficients 
ablation-cluster ion theory 


for the meteor 
of sodium 


Reaction 

Reference Reaction 

Number 

Clustering 

Rl Na* + N 2 + X ^ Na*-N2 + X 


R2 

Na* + COj + X ? Na**C02 + X 


Switching 

H3 

Na**N, + CO 2 + Na‘^*C02 + Nj 

M 

Na^'COj + H 2 O + Na‘^*H20 + CO 2 


Recombination 

ns 

Na^'CO^ + e Na + CO 2 

R6 

Na^-ll20 + e’ -V Na + 


Neutral Reactions 

H7 

NaO + 0 Nil 4 O 2 

m 

Wa + 0^ -► NaO + O 2 


Ionization and Charge Transfer 

I 

Na + hv Na"^ + e** 

m 

Na + 0^^ ^ Na^ O 2 

lUO 

Na + NO* -► Ma* + NO 


Rate Coefficient 


fej = 2 X 

fc_l = 9 X 
fc_2 = 6 X 

= I X 
= iX 

= 1 X 

= 1 X 

Q 



«ll ^ -1 

k = 3a9 X 10 cm’^ s 

/ ^ r ,,.“12 3 -1 

= 6.S '< 10 cm s 

o 


J = 1 .7 X 10"^ s"^ 

= 6 X 10 "'*^ cni^ s~^ 

“ _T1 -I .1 

fejjj = a X 10 cm s 
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been performed only once. Goldberg et at. [1977] observed over Esrange, 
Sweden two peaks in the Na^ profile occurring at 90 and 95 kmj the Na profile 
also exhibited a double peaked structure, with the peaks occurring at the 
same altitudes. The [Na]/[Na*] ratio was near 2 at these peaks. While these 
data indicate that [Na] « [Na ] much data are lacking in the report, parti- 
cularly on the behavior of the cluster ions. 

iVhile simultaneous measurements which give us enough information have 
not been performed, we are fortunate in having two independent measurements 
which nearly satisfy the criterion. During the Geminids meteor shower on 
the night of December 13-14, 1971, Hake et at, [1972] measured the concen- 
tration profile of atomic sodittm over Menlo Park, California [37°27*N, 
122°10*W). The following day, December 14, 1971, at 12:11 UT Zbinden et al. 
[1975] launched a rocket -borne mass spectrometer at Sardinia (39**36'N, 

9®26'E) and obtained Na"*" and Na’*’*H20 profiles in the 80-110 km region. 
Selected data of Zbinden et at. [1975] are shown in Figure 7.9. In a later 
publication, Hidalgo [1977] derived a temperature profile based on these 
measurements by employing the temperature dependence of a N 0 ^*N 2 reaction 
scheme. The derived temperature, T = 22S“K between 80 and 86 km, is a 
necessaiy input into the cluster ion scheme. 

Major constituent concentrations are taken from CIRA [1972]. Although 
this is a reference atmosphere, problems should not occur since N 2 is not 
believed to be highly variable. The CO 2 and H 2 O concentrations were based 
on constant mixing ratios of 300 ppm and 1 ppm respectively. Values for 
NO and O 2 were taken from Zbinden et at. [1975] ; but since they did not 
simultaneously measure [e"] , the electron density profile for quiet sun con- 
ditions is taken from Medhtty et at. [1972] . 
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Chemical equilibrium may be assumed for these calculations since the 
cluster ion reactions have small time constants in comparison to the time 
scales associated with eddy and molecular diffusion. Using the Na profile 
of Zbinden et at. [1975], Na, Na'^*N 2 , Na'^«C02, and Na'''*H20 values are calcu- 
lated every kilometer. The results of these calculations are shown in Figure 
7.10. 

The agreement is very good. First the peak concentration of sodium is 
predicted properly, and the peak occurs at the proper altitude. It is not 
expected that the shapes of the two sodium profiles should agree, since 
dynamical effects can undoubtedly perturb the layers. However since lidar 
measurements indicate no diurnal variation in the sodium layer [Gibson and 
Sandfovdj 1972], the profile of Hake et at* [1972] should be a good estimate 
to the general features of the sodium layer that existed during the Zbinden 
et a%, launch. 

•f 

Furthermore the cluster ion scheme accurately predicts the Na *K20 pro- 
file near 80 km obtained in the mass spectrometer measurements. Since sig- 
nificant amounts of Na'''*H20 are also predicted near 90 km, it may be improper 
to identify this ion as as Zbinden et at. have done. The cluster ion 

Na^*C02 is calculated to remain less than 0.3 cm“^ except in the altitude 

-3 

range 89-98 km where values are close to 1 cm . Only at the peak of the 

+ 3 

Na layer, is Na "CO 2 estimated to be 2-3 cm . It is therefore unlikely 

that Zbinden et at. [1975] detected Na^»C02 and this is reflected in their 

+ -2 -3 

data. The cluster ion Na *N 2 is calculated to remain below 10 cm at all 

altitudes and obviously was not detected. Finally on the basis of Figure 

II. 4 which predicts that [NaO^] 10“^[Na‘^] at 85 km, it is very unlikely 

+ 4 * 39 *** 

that mass 39 was NaO but rathei* K , 


ALTITUDE (km) 



Figure 7.10. Profiles obtained from the cluster ion 
scheme using the Na^ profile of Zbinden 

et at. [1975] as input. Comparison of ! 

the calculated Na profile to the experi- " 

mental profile measured by Bake, et at. ^ 

[1972] shows good agreement, | ^ 

• I 
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In sununary the simple chemical equilibrium calculation for the cluster 
ion scheme accurately predicts the concentrations of both neutral and ionic 
foxms of sodium. The final test of the cluster ion scheme will concern 
itself with the ability of the scheme to explain the more common features 
of the sodium layer, 

7.5 Model CaZaulations 

7 . 5.1 The effect of variations in second and third body concentrations. 
Assuming chemical equilibrium, the expression for [Na] becomes 

fe_[0][Na0] + {?:.[Na'^*CO,] + ?ift[Na'^*H„0] } [e'] 

[Na] = -d 2 2 J (7.30) 

k^iO^] * I + kg[Q^ ] + fejQ[N0 ] 

It is apparent that [Na] is proportional to the atomic oxygen and electron- 
density concentrations. Thus the sharp decrease on the bottomside of the 
layer is related to the electron density and atomic oxygen ledges which 
exist in the same region. 

An examination of the daytime time constants indicates that the neutral 
Na system comes into equilibrium slightly faster than the neutral-ion system. 
Yet the latter is still of a shorter time scale than the process of eddy 
diffusion. Thus the chemical scheme behaves as though sodium acted as a 
reservoir. It is in constant equilibrium with NaO, and so the supply of 
sodium by the recombination of the cluster ions must be balanced by the loss 
of sodiiun by photo ionization and charge exchange. 

All of the reactions depend on assumed profiles; but since the important 
cluster ion processes involve three body reactions, and the ambient neutral 
atmosphere is not expected to vary greatly with time, the chemical scheme 
should be applicable to all seasons. 
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However the effect of variations in [C02]» [H2O], and [«"] should be 
investigated. The case for CO2 will be considered first. Assuming chemical 
equilibrium and using the notation of Table 7.3. 


*:,[Na'"][N2][X] 

*^2^ = r^[x] Tl:3lC02T 


(7.51) 


The switching reaction is so fast, that ^_i(X]. Hence [Na^*N2] 

varies inversely as [CO2] . The concentration of Na*»C02 may be bitten as 


[Na *€02] = 


fe2[Na’"][C02HX] + fc3[Na'"*N2] [CO2] 


fc_2(X] + k^[e~} * fe4[H20] 


(7,32) 


Since [Na *002] is more strongly supplied by the s;tfitching of CO2 id.th 

*}' ■}■ 

Na «N2 than it is by the direct clustering of Na with CO2, 


[Na**C02] « /c3[Na**N2] [CO 2 ] 


(7.33) 


which implies ivith the aid of equation (7.31), that [Na*»C02] is independent 
of CO2 concentrations on the order of 300 ppm. 

The case for H2O is a bit different. Mixing ratios on the order of 1 ppm 
are commonly assumed in cluster ion calculations {Reid^ 1977; Thomas^ 1975], 
Denoting the rate coefficient for the clustering of Na^ with H2O as ?c, one 
obtains 

fe [Na'"*CO ][H 0] + fe[Na*][H 0][X] 

[Na ‘H O] = = = 1 (7.34) 

fegM + ][X] 


It is assumed in equation (7.34) that the switching of H2O with Na^*N2 pro- 
duces less Na «H20 than the indirect route involving the switching of CO2 
with Na^*N2* This is true since [CO2] » [H2O] . The main question concerns 
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the numerator of equation (7.34) and whether 
^ [Na^.CO,] 

.J L. » 1 

fc[Na^IfX] 


C7.35) 


Referring to Table 7.1, it is apparent that ~ 10^ ^[X] where a temperature 
of 200°K is assumed and [X] = 7 x 10“^^ citi”^ at 90 km. Since all calculations 
show that [Na^*C02] - lO"^ [Na*] , equation (7.35) holds and the direct clus- 
tering of H 2 O \tfith Na should not be important. Below the electron-density 
ledge (-85 km), fegLe"! « k^[H20] and so 


[Na^'.COj] . 


a ? 85 km 


(7.36) 


Hence the production of [Na^*H20] should be independent of IH 2 O] . Above 95 
km, kg[e~] becomes comparable to k^[H20]. However, calculations indicate 
that [Na^*H 20 ] still remains relatively independent of [H 2 O] for water mixing 
ratios of up to 10 ppm. 

The scheme will be affected slightly by changes in [e~] . In e;q>ression 
(7.32), the denominator is controlled by the relative magnitudes of ^^[6"*] 
and k^[H20]. The altitude at which [Na^*C02] is affected by the electron 
density is determined by the condition 



Also from (7.34) [e~] becomes important when 


[f] - ^11 
Tn^W" ^ 


-25 

3 X 10 


(7.38) 


The altitudes corresponding tc (7.37) and (7.38) are near 95 km and 85 km 
respectively. Above these altitudes [e~] variation becomes less important 
for producing Na because increases in [e~] create decreases in Na^*C02 and 
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Na 'H20. However below 85 km, the production of Na is increasingly dependent 
on [s~l since [Na^*C 02 ] and [Na^*H 20 ] are controlled more by switching and 
three body leactions respectively. Hence one would expect the bottomside of 
the Na layer to be moderately affected by the electron-density profile, par- 
ticularly near the ledge. 

In fact a knowledge of the value for [e“] near the ledge may be an im- 
portant parameter for determining the value of Since the Zbinden et al. 

launch was performed under anomolous conditions, [e~] may have been enhanced 
in the 80-85 km region by as much as a factor of 3. A value of = kj^/2 
is then found to be more appropriate than the value used in the previous 
calculations , This higher value would also yield a more rapid 

bottomside falloff in [Na] under normal conditions. 

7.5.2 Seasonal variation. It is well known that [Na] increases during 
the winter and that the peak of the layer may lower slightly (Pigure 1.1). 

The cluster ion reactions are temperature sensitive, but unfortunately in 
the wrong directions. That is, as temperature increases, fewer cluster ions 
are produced which reduces the rate of dissociative recombination. Hence 
the temperature dependence implies that [Na] should decrease in the winter. 

If it is assumed that Na has a Gaussian profile given by 

r 100 exp [-((s - 90) /3)^] cm~^ for [Na^] > 1 cm"^ 

[Na ] (7.39) 

L 1 cm otherwise 

then this increase in summer is evident in Figure 7,11. Notice that the 
Gaussian profile is not an accurate description, since large quantities of 
Na •H 2 O are predicted. A realistic Na prefile would fall off much more 
rapidly on the bottomside. Yet the technique still has merit. 

The assumption used in deriving Figure 7.11 is that the Na^ layer does 
not move. Three body reactions are so important to this scheme, that the 


ALTITUDE (km) 
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equilibrium value of [Na]/[Na*] will be highly altitude dependent. Let us 
assume that Na’^ takes on the same Gaussian shape but that it is centered 
about 95 km in the summertime (T - 150^ K) and compare this to a Na* layer 
centered about 90 km in the winter (T - 225“ . The two resulting profiles 
are sho^vn in Figure 7.12 and indicate a lowering of the peak of the layer by 
about 5 km and a column increase in sodium by a factor of 5. The factor is 
more than that observed by Megie and Blamont [1977] and hence the 5 km sepa- 
ration may be reduced to 3-4 km which is in good agreement with the profiles 
shown in Figure 1.1. The principle then is, if the Na profile does not 
change its amplitude or shape with season, the peak of the Na layer deter- 
mines the peak of the Na layer; and the altitude at which the peak occurs 
strongly affects the equilibrium value of [Na] . It is interesting to note 
that the [Na]/[Na^] ratio varies from 30 to 90 for this example (stunmer to 
winter] and hence this fact may explain the wide range of estimates for this 
ratio [Navd.S'ij. 1968]. 

Finally a seasonal variation in the ambient background density would 
produce the same effect as that produced by the movement of the Na^ layer. 
From CIRA [1972] the atmospheric concentration at 90 km is about 2 times 
that at 95 km. Hence an increase in number density of a factor of 2-3 during 
winter is required. Model calculations show this is plausible (Tables II. 1 
and 1 1. 2). For this assumption seasonal changes in the peak of the Na layer 
would not be so pronounced unless dynamical effects became more important 
during winter. 

7.5.3 Diiamal variation. Motions play an increasing role at nighttime. 
Since photoionization of Na is removed from the scheme, and charge exchange 
is severely reduced because of the decrease in [02^] and [NO^] , a path no 


ALTITUDE {km) 


SUMMER 


WIMTER 



CONCENTRATION {crrT^) 


Figure 7.12 The effect of lowering the Na layer from 
95 km in the summer to 90 km in the winter 
yields an increase in the total column 
content of Na. A similar effect would be 
caused by an increase in the ambient neutral 
density during winter, although the altitude 
difference would be much reduced. 
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longer exists between Na and Na . Rather Nu should be depleted by the clus- 
ter ion reactions and either be converted to Na or be removed by the sink in 
the 80-85 km region. However the [NaJ/[Na ] ratio is high enough that the 
column abundance of Na* should cause a negligible increase in the column 
abundance of Na. A diurnal variation should not be observed. The layer at 
this time may be an effective tracer for atmospheric motions and could be 
used to monitor wave motion having spatial scales less than the width of the 
sodium layer (10-20 km) . The lack of diurnal variation is supported by lidar 
measurements [Gibson and Sandford^ 1972] as well as by spectroscopic measure- 
ments of sodium dayglow [Bumett et at, 1972, 1975], 

7.6 The Interpretation of Some Experimental Results using the Meteor 

Ablation-Cluster Ion Theory 

7.6.1 Variation of th^ topside scale height. Although the cluster ion 
scheme can explain a seasonal variation in the sodium column content and a 
lack of diurnal variation in the layer, one of the questions to be explored 
by this work has been the validity of the corkscrew mechanism. A major 
concept is that the corkscrew mechanism is responsible for the continual 
supply of Na to the sodium layer. The experimental observations of enhance- 
ments in the sodium topside which move downward \ri.th a velocity of about 
50-100 cm s”^ lend support to the presence of the corkscrew mechanism. The 
chemistry developed in this chapter shows that chemical time constants are 
still rapid at 95 km for the conversion of Na* to Na and hence movement of 
the Na layer will be mirrored in the Na layer. Results from other groups 
studying the sodium problem also lend support to the corkscrew mechanism. 

A characteristic feature that Megie and Bl-amcnt [1977] observed is the 
variation of the topside gradient with time. At night the layer is under 
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dynamical control since the ion chemistry gradually disappears as Na is 
depleted. Under these conditions one would expect the layer to be broadened 
under the influence of eddy diffusion. In 90% of the cases, Megie and 
Blamont found that a proper choice for the eddy diffusion coefficient could 
explain this broadening of the layer and in fact they obtained average values 
for K using this technique. 

The puzzling fact is that in 10% of the data, the topside scale height 
decreases with time particularly at the beginning of the night, and they 
postulate that this results from a source of sodium near 91 km. It may well 
be that this source is the dumping of ions into the region by the corkscrew 
mechanism. Near sunset, the 90-120 km region still has a supply of long- 
lived Na^ ions remaining from the daytime ionization of Na, It is expected 
that these ions remain layered and are sitfept out of the region during the 
night. Evidence of electron-density gradients have been seen in the 90 km 
region during the nighttime monitoring of the layer (Section 6.5) and may 
be indicative of metallic layers of ions. The maximum amount of ionization 
available naturally occurs near the beginning of the night, however higher 
altitude Na layers may enter the Na layer ttfell into the evening according 
to tliis theory. 

The time scale for the conversion of Na^ to Na is limited by the time 

constant for the formation of Na At 92 km this time constant is roughly 

20 min. A layer of Na capable of maintaining a concentration of 50 cm 

having a thickness of 3 km will then provide a source having an integrated 

4 -2 -1 

strength of 1.2 x 10 cm s , Interesting enough this compares favorably 

with a lower limit value on the source of 2 x 10^ cm"^ s~^ assuming 
6 2 

if = 2 X 10 cm s""^ [Megie and Blamont^ 1977] . During enhancements in the 
meteoric input, it is likely that the integrated source could be greater. 
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Eccnson and Donaldson [1967] explained the influence of the photoioniza- 
tion rate on the topside scale height. Their method does rot assume a source 
for sodium. If a constant source in the same altitude region is postulated 
then the effect works against photoionization and tends to increase the 
scale height. Since layers of ions driven by the tidal oscillations are 
moving through the sodium layer, diurnal variations in the topside layer 
should occur. 

7.6.2 The ratio. The cluster ion scheme may have an impact on 

the seasonal variation of the [Na]/[K] ratio as observed by Megie et at. 
[1978], K is not as reactive as Na and hence should not show as strong of a 
temperature dependence in the cluster ion-chemistry scheme. Furthermore ion 
separation is believed to le produced by the wind shear mechanism at higher 
altitudes (- 105 km] [Chimonas, 1969]. The effect is presumably lost near 

90 km but may still occur. Since the seasonal variation is tied to the tem- 
perature effects on reactions and the position of the Na* layer (or alter- 
natively on the neutral atmospheric concentration] , these temperature and 
layer effects need to be considered for the case of K before it is assumed 
that chemistry cannot explain the seasonal changes in the [Na]/[K] ratio. 

7.6.3 Removal of sodium during the day. Finally Megie and Blamont 
[1977] notice that after meteoric activity, the enhancement in the sodium 
layers is removed sometimes within 12 hrs. It should be eirphasized that 
the sink for sodium varies directly as [Na*] and is therefore continuously 
adjustable. Since [Na*«H20] varies directly as [Na*] , an enhancement in 
[Na ] produces enhancements in [Na *H20] which causes a greater loss of 


sodium. 
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7.6,4 Uct meteorie input. In this theory the source is assumed to be 
Na^ resulting from the photoionization of ablated material containing sodium. 
This is based on the assumption that meteoric material is ablated in neutral 
form. However, Lebedinets ccnd Shushkova [1970] calculated that near 110 Jon, 
meteor ionization is important. The effect of change transfer oi’ photo- 
ionization in the ablation process may also increase the input of Na^. If 
these facts are true then a source for Na would result at night, indicating 
the corkscrew mechanism should supply the Na layer at night. These ideas are 
not inconsistent with observations or theory. 

7 . 7 Summary 

In this chapter a simple ion chemistry for sodium has been developed 
under the assumption of daytime chemical equilibrium. The chemistry includes 
switching and recombination of cluster ions but does not consider their 
photodestruction. The main cluster ion is found to be Na ‘^ 2 ^' chemis- 

try is then used in conjunction with a meteor ablation source to explain the 
dynamical features of the sodium layer. 

A seasonal variation may be produced by the model and lack of a diurnal 
variation is evident. Current problems in the sodium study are discussed in 
the light of this theory. 
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8. SUMMARY AND CONCLUSIONS 

8.1 Surmcay 

A laser radar system has been built and operated to monitor the night- 
time sodium layer over Urbana. The signal-to-noise ratio has been increased 
sufficiently to obtain accurate relative concentration profiles of the sodium 
layer in 10 minutes. Thus the time history of the sodium layer may be moni- 
tored throughout the night. 

In general the layer is observed to be highly variable during the winter 

time. Wavelike structures have been found which appear to have downward 

-1 - 1 

phase velocities of from 50 cm s to 100 cm s and vertical wavelengths 
from 3 to IS km. The altitude of the peak of the layer is often found to be 
lower in the winter than in the spring and has appeared as low as 83 km. The 
bottomside usually falls off more rapidly than the topside. 

Simultaneously with the collection of data, a theory for the sodium 
layer was developed. The main problem areas about the aeronoray of sodium are 
the following: 

1. Origin. Suggested sources have been meteor ablation or thin 
aerosol layers, either of extraterrestrial or terrestrial origin. 

2. Boundary gradients. The sodium layer falls off quite rapidly on 
both the topside and bottomside. 

3. Seasonal variation. There is an enhancement in winter at midlati- 
tudes and a possible lowering of the peak of the layer at this time. 

4. Diurnal variation. Lidar measurements and spectroscopic dayglow 
measurements indicate little diurnal variation. 

5. Nootumal variation. The sodium layer seem - to show little sys- 
tematic variation during the night. 

Sporadic effects. The layer shows considerable short-term variability. 




6 . 


197 


8.2 Conolusions 

To explain all the aboveraentioned properties a meteor ablation-cluster 
ion theory was developed. This theory postulates that the source of sodium 
is meteor ablation in the 90-120 km region. The deposited sodium is then 
rapidly ionized in the daytime by photo ionization or charge exchange with 
0^ and NO . It is shotm that the lifetime of Na is long enough for Na to 
be layered by an interaction of the tidal winds with the geomagnetic field 
lines (corkscrew mechanism). The layers collect near the nodes of the wind 
profile and are then transported downward with the characteristic phase velo- 
city of the wind profile (-30-100 cm s'^). Below 100 km the ions no longer 
can maintain the same down\tfard velocity of the wind profile as the neutral - 
ion collisional frequency increases. They gradually are deposited in a re- 
gion between 90 and 95 km. 

A cluster ion chemistry was developed to explain a rapid conversion of 
Na^ into Na. This scheme had several advantages over the earlier metal oxide 
ion chemistry. The major sodium constituents then become Na and Na at higher 
altitudes with NaO and Na^»H.,0 playing an increasing role near the bottom of 
the layer. A sink is provided in the form of higher order cluster ions which 
possibly attach themselves to aerosols and thereby undergo sedimentation to 
the lower atmosphere. 

The following explanations are offered to the problem areas; 

1. The origin of sodium is taken to be meteor ablation based on the 
correlation of sodium column increases with meteoric activity. Also metallic 
ion abundances have relative ratios close to what is believed to exist for 
the neutral metals in meteoric material. 
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2. The topside gradient is controlled by photoionization, eddy diffu- 
sion, and the dumping of ions by the corkscrew mechanism. The bottoraside of 
the layer is chemically controlled and is related to the electron density 
and atomic oxygen ledges known to exist in the region. The bottomside fall- 
off of Na also depends on the atmospheric conditions which favor the formation 
of higher order cluster ions below 85 km. 

3. The seasonal variation in olumn content is a result of the sensi- 
tivity of the cluster ion reactions to atmospheric density. It may be re- 
la ted to a repositioning of the Na layer with season or an actual change in 
the atmospheric density with season. 

4. The lack of diurnal variation follows from the large [Na]/[Na ] 
ratio. Na is gradually depleted through the night. 

5. The layer shows no systematic variation during the night because 
the loss and production mechanism of the sodium column content are turned 
off. 

6. The short-term ’/ariability is undoubtedly related to the turbulence 
known to exist within the region. The influence of gravity waves and hori- 
zontal transport are undoubtedly important. 

8.3 Suggestions for Future Work 

8.3.1 Eospe/Hmental aspects. Several correlative measurements are sug- 
gested by this work and a brief discussion of each will follow. In addition, 
some obvious improvements in the lidar system are proposed. It has become 
apparent that lidar measurements taken alone are insufficient to solve the 
aeronomy of the sodium layer, and the following topics present the core ex- 
perimental effort which should be undertaken. 

1. Sirmttaneous measurements of electron density. The theory predicts 
that the bottomside of the sodium layer is tied to the electron-density ledge. 
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Lacking a method for measuring electron density at night * we probed the 
atmosphere with a 2.66 MHz signal while simultaneously monitoring the sodium 
layer. Reflections of the radio signal indicated electron density gradients 
in the same region as the sodium layer. This then suggests that measurements 
of the electron density be made with an ionosonde or incoherent- scatter radar, 
simultaneously with the lidar experiments, for precise monitoring of the 
ionization. 

2. Simultaneous measurements of neutvat winds. A meteor radar system 
exists at the Aeronomy Laboratory Field Station and can monitor winds in the 
sodixnn region. It would be important to find any correlation between the 
waves identified by meteor radar and waves observed by laser radar. Further- 
more the meteor radar system could monitor the east-west winds which most 
strongly affect the vertical ion transport below 100 km. By this joint experi- 
ment it may be possible to determine if the wave motion in the sodium layer 
is related to wave motion of the ambient atmosphere or rather to the cork- 
screw mechanism. If the former case were true, sodium could be used as a 
tracer to observe atmospheric wave motion, 

A feasibility study conducted at the Field Station revealed inter- 
ference between the systems. Also the meteor radar unit was observed to 
increase the ambient temperature of the room housing the laser to the point 
at which thermal stability of the tuning optics became a problem. If the 
systems are operated alternately throughout the night to avoid the inter- 
ference, the data may suffer because each system depends on the integration 
of returns for improved signal- to-noise. Yet the experiment can be done 
with existing equipment and should be attempted in the near future to see if 
the above objections can be removed. 
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3. Mobile Hdar applioatiom . There are other experiments which re- 
quire lidar mobility. Simultaneous lidar monitoring with a rocltetbome posi- 
tive ion mass spectrometer measurement could be accomplished over Wallops 
Island, Virginia. This would help to determine more exactly the relation 
between Na and Na*, and thus the cluster ion scheme could be refined. Fur- 
thermore a mobile lidar could be taken to the equator where the sodium layer 
is predicted to be rather broad. Neither of the measurements have ever been 
done. In any event work should progress toward reducing the volume of the 
lidar system in preparation for the day when a mobile lidar is feasible. 
Measurements from an airplane might also be considered. The weather over 
Urbana many times of the year is not conducive to frequent operation; an air - 
plane would give us the opportunity to fly above the clouds and avoid much 
of the attenuation due to molecular and aerosol scatter. 

4. Steerable lidar. Although this work assumes that the movement of 
sodium is in the downward direction, horizontal transport may be present. A 
steerable lidar system could measure the effect of horizontal transport; how- 
ever, it should be noted that the experiment has been performed over Winkfield 
(51.4‘’N, 0.7*’W) [Thomas et al. ^ 1977]. Hence the experiment should be under- 
taken only if substantial improvements over success of Thomas et al. can be 
guaranteed. It should be noted that safety regulations and common sense dic- 
tate that the lidar should be turned off when planes are in the area if the 
zenith angle of the beam is great enough, A radar unit may then have to 
accompany the steerable lidar for detection of local air traffic. 

5. Measurement of potassium. Measurements of other constituents, 
especially potassium may be conducted; but the measurements involve much lower 
signal levels. This study indicates that more work needs to be done in the 
sodium area, and a measurement of K may only dilute our ability to study the 
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Na layer. Megie et al» [1978] have conducted simultaneous K and Na measure** 
raents; and although their statistics could be improved it will be hard to 
better their results in the near future. 

6. Daytim masueernents . Daytime lidar measurements of Na have not 
been repeated. These measurements made a strong intact on the sodium theory 
because they showed no diurnal variation in the layer [Gibson and Sccndfordj 
1972]. Tae measurements should be repeated. Careful monitoring of the twi- 
light sodium layer could also be accomplished aiid would provide information 
on any rapid structural changes in the layer during this time. A smaller 
telescope with a narrow band filter would be required and increased laser 
output desirable. 

7. Laser improvements, A general upgrading of the lidar is required. 
The signal-to-noise can best be increased by enlarging the aperature of the 
receiving telescope; this plan is being implemented. In addition laser im- 
provement is necessary. Wavelength emission stability and spectral linewidth 
monitoring improvements are the two areas of prime concern. At this time we 
cannot compare returns from one night to the next with any accuracy because 
of the iincertainty in the spectral linewidth measurement; hence our approach 
has been to observe relative changes in the layer during a night. Laser 
stability can be accomplished bes*: by a two etalon tunfng in which the tuning 
elements are temperature controlled. A Lyot filter has been used with success 
by some groups as a replacement for the crude tuning etalon and should be 
considered. A high resolution spectrometer is sorely needed to measure the 
spectral width of the laser line emission. If the linewidth cannot be moni- 
tored accurately, then absolute measurements of sodium cannot be made, A 
high quality spectrometer is a necessity for excellent data. Ideally line- 
width measurements should be recorded with the data. Some thought should 
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also be given to the portability of such a spectrometer. Finally it has often 
been proposed that an amplifier be added to the current oscillator. Before 
this is implemented, a transverse flow dye laser should be considered. These 
lasers can provide high average power as well as high peak powers [Jetftwa, 

1978] and may be commercially available in the near future. The spot geometry 
is not circular, but this should cause no problem. Also any second generation 
laser must reduce the duration of the pulse to less than 500 ns for better 
range resolution. 

8. Measurement of ohemioaZ rate aoeffioients. Finally one area of 
experimental research is directed toward the scientific community in general. 
Tl\ere is an urgent need for the laboratory chemists to measure the rate coef- 
ficients for many of the sodium reactions at mesospheric temperatures. It is 
crucial to determine if the important neutral sodium reactions proceed by the 
electron jump mechanism. Also the estimates for the rate coefficients of the 
cluster ions are accurate to an order of magnitude at 300”K. The extrapola- 
tion of these rates to lower temperatures therefore has even a greater chance 
for inaccuracy. 

8.3.2 Tkeoretioal investigations. The foundational work for the 
meteor ablation-cluster ion theory has been performed and several areas have 
now been opened for future research. 

1. The effect of higher order cluster ions and the eleotron density 
ledge near 85 hn. The cluster ion chemistry should be expanded to model more 
accurately the sink. In particular the effect of higher order clustering 

.J* 

reactions on the Na »H20 topside needs to be investigated. It would be in- 

' 4 ' ^ 

teresting to determine if the Na clustering reactions can be tied to the NO 
clustering reactions. 
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This work may lead to a further understanding of the electron- 
density ledge near 85 km. This region is particularly interesting because of 
strong dynamical changes in the region and temperature gradients. 

2. Seasonal variation of [^Va]/[/^]. The [Na]/[K] ratio variation with 
season needs to be explored. Is this variation caused by changes in the source 
for Na or by chemical differences between the elements? The cluster ion reac- 
tion scheme may invalidate the assumption that K behaves in the same manner 
that Na does. 

3. Chemistry of meiattio species, ^ general study of metallic species 
in general would be a significant step forward. Although Staider [1969, 1970] 
provides an excellent review, the work needs to be updated in the light of 
cluster ion reactions. It may be possible that some metals follow an oxide- 
ion chemistry while others obey the cluster ion reactions. 

4. Atmospheric simulations. Finally atmospheric simulations which 
emphasize dynamical effects on the sodium layer need to be explored. It is 
important to be able to eventually model the sodium layer to see if experi- 
mental observations can be duplicated. The chemistry for this model though 
may have to be fairly complete. Sodium probably is tied to other clusters 
in the region and a complete neutral chemistry may also be required. Atmos- 
pheric motions would include the wind shear mechanism, gravity waves, hori- 
zontal transport, and changes in eddy diffusion. 
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APPENDIX I 

THE INITIAL LIDAR RECEIVER 

To obtain the first sodium returns over Urbana as quickly as possible, 
the first lidar receiver used the PDP-15 computer facility at the Aeronoray 
Laboratory Field Station. This computer is used with several ground-based 
experiments for data collection and its inclusion into the lidar system was 
accomplished by minimizing hardware whenever software could accomplish the 
same results. 

Since we required a high data rate (100 kHz) interface with direct 
memory access, the custom-built analog- to- digital (A/D) converter interface 
was chosen as the point of entry for our data. Existing data collection 
software was used by substituting 8 bits of lidar data for the output buffer 
of the A/D conveirter. Softitfard modification provided suitable masking of the 
unused interface bits; and when the BCD counter was used, the software con- 
verted the BCD to binary. The block diagram of this system is shoim in Figure 
I.l. The flow of information is essentially the same as that described in 
Chapter 5. Schematic diagrams of important circuits are indicated in Figures 
1.2 - 1.4. 

Because of the versatility of the PDP-IS and the ample storage capa- 
bilities, the computer enabled us to define problem areas and investigate 
data collection routines. A serious problem for this system resulted from 
combining the asynchronous BCD upcounter (University of Illinois, Elementary 
Particle Physics Research, 10 MC Scaler S102) with the sjuichronous computer 
system; this arrangement resulted in the occasional sampling of the upcounter 
while it was shifting thus yielding errors. When adjacent bins were sub- 
tracted to obtain the counts per bin, these sampling e3crors created errors 
in two adjacent bins. The ripple through BCD counter also provided ample 
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Figure 1.2 Line driver and line receiver circuitry for the 
receiving system which employed the BCD upcounter 
[Scaler S102). The actual circuit consists of 4 
copies each of circuits (a^ and (b) . 
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opportunity for sampling errors to occur since it remained in a shifting 
state for an extended period of time. This counter was later replaced by 
two cascaded 4-bit 32-MHz counters which required only 25 ns for all outputs 
to shift (Figure 1,3). These counters reduced the sampling errors substan- 
tially but did not solve the problem. The remaining errors were removed by 
software, but data were lost in the process. 

The eventual need for a mobile system, coupled with limited manpower, 
suggested that a new system revolving around the HP 2114A computer should 
be built. It was judged more efficient to incorporate the improvements into 
the new system rather than to continue to modify the old system. Furthermore, 
the range bin size for the PDP-15 system was limited to 1.5 km which is stand- 
ard for all experiments using that computer. Thus the desire for improved 
range resolution was also a determining factor in the construction of a new 
receiver system. This new system is now in use and is described in Chapter 5. 
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APPENDIX II 

SODIUM NEUTRAL AND METAL OXIDE ION CHEMISTRY 
II.l ^eutpdl Sodiim Chemistry 

The neutral sodium chemistry has been reviewed by many authors [Brom, 
1973 j KviftSf 1973; Hunten^ 1967]. In general theories have suggested that 
sodium reacts with various forms of oxygen, nitrogen, and hydrogen. Although 
most of the reaction rates have not been measured for the sodium chemistry, 
estimates may be made on the basis of laboratory experience Involving reac- 
tions which are closely similar to those of sodium. 

A map of the neutral species chemistry is shown in Figure (11.1). The 
map is useful for indicating the multitude of interactions thought to exist; 
although no study to date has incorporated all of these reactions. In fact 
most of the reactions are of secondary importance, given the current esti- 
mates for the rate constants. This is clearly indicated if we look at the 
time constants for these reactions. The assumed profiles for the second and 
third bodies of the reactions are taken from Shimasaki and Laird [1972] and 
axe shown in Tables II.l and II. 2. The specific reactions and adopted rate 
constants are listed in Tables II. 3 and II. 4; these tables also include the 
calculated time constant. Photo ionization of NaO and Na02 is neglected be- 
cause of the lack of knowledge about these reactions. The fast bimolecular 
reaction of Srivastava and Skukta [1970] which produces the excited species 
NaO* from Na is also neglected since it is likely that NaO* will dissociate, 
returning Na rather than producing NaO, 

It is apparent from these calculations that the neutral chemistry is 
dominated by the interchange of sodium with sodium monoxide via reactions 
involving ozone and atomic oxygen. This conclusion is in agreement with the 
numerical calculations of Megie [1977] and for the adopted reaction scheme 






Table II. 1 Adopted profiles for major constituents during siaraner [Shimazaki- and Laivd^ 1972]. 

To indicate diurnal changes, the top entry, if present, represents the daytime value 

(12 hrs., X = 45°] while the lower entry is the nighttime value (0 hrs., x = 135°). 

Omissions imply very low concentrations. The entry 4.1(10) denotes a concentration 

£ jt 1 11^10 -3 

of 4.1 X 10 cm . 

-3 

Concentration (cm ) 


Constituent Altitude (km) 



70 

80 

90 

100 

110 

120 

0(^P) 

4.1(10) 

1.0(7) 

6.0(10) 

1.7(10) 

6.0(11) 

6.0(11) 

2.8(11) 

2.8(11) 

1.0(11) 

1.0(11) 

6.0(10) 

6.0(10) 


4.6(14) 

1.0(14) 

1.0(13) 

1.0(12) 

1.0(11) 

1.9(10) 

°3 

4.6(9) 

1.3(10) 

3.6(8) 

3.6(8) 

1.0(8) 

7.7(8) 

1.0(5) 

7.7(6) 

1.0(5) 


OH 

1,3(6) 

7.7(5) 

7.7(5) 

2.8(6) 

3.6(3) 

3.6(4) 




H 

1.5(7) 

1.7(8) 

1.3(8) 

4.6(7) 

4.6(7) 

1.0(7) 

1.0(7) 

3.6(6) 

3.6(6) 

2.8(6) 

2.8(6) 

«2 

2.2(15) 

4.6(14) 

6.8(13) 

7.7(12) 

1.0(12) 

2.8(11) 

X = N„ + 0 

2.7(15) 

5.6(14) 

7.8(13) 

8.7(12) 

1.1(12) 

3.0(11) 


Table II. 2 Adopted profiles for major constituents during winter [after Shimazaki cmd Lait>d, 1972], 
The entry format is the same as in Table II. 1. 


Concentration (cm ) 


Constituent Altitude (km) 



70 

80 

90 

100 

110 

120 

0(^P) 

1.5 CIO) 
1,0(5) 

7.7(9) 

1.7(9) 

4.6(10) 

4.6(10) 

3.6(10) 

3.6(10) 

1.1(10) 

i.l(lO) 

8,8(9) 

8.8(9) 

°2 

3.2C14) 

6,8(13) 

1.0(13) 

2.2(12) 

3.6(11) 

1.0(11) 


6.0(8) 

1.0(10) 

2.2(7) 

7.7(7) 

6.0(6) 

1.0(8) 

1.0(5) 

3.6(6) 



OH 

2.2(5) 

1.7(5) 

1.7(5) 

1.0(6) 

1.7(3) 

2.2(4) 




H 

4.6(6) 

6.0(7) 

4.6(7) 

2.5(7) 

2.5(7) 

1.0(7) 

1.0(7) 

2.8(6) 

2,8(6) 

1.7(6) 

1.7(6) 

^2 

1.7(15) 

3.2(14) 

6.8(13) 

1.0(13) 

2.8(12) 

6.8(11) 

X + N„ + 0„ 

2,0(15) 

3.9(14) 

7.8(13) 

1.2(13) 

3.2(12) 

7.8(11) 
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Table II. 3 Summer time constants for neutral sodium reactions. 
The entry format is the same as in Table II. 1. 


Loss of Na 



70 

80 

90 

100 

110 

120 

Na ♦ Oj + HaO + 0^ 

6.5(.12) 

[M] 

3.3(1) 

1.2(1) 

4.3(2) 

4.3(2) 

1.5(3) 

2.0(2) 

1.5(6) 

2.0(4) 

1.5(6) 


Na * 0 + X NaO + X 

7.0C-33) 

[K] 

1.3(6) 

S.3(9) 

4.3(6) 

1.5(7) 

3.0(6) 

3.0(6) 

5.9(7) 

5.9(7) 

1.3(9) 

1.3(9) 

7.9(9) 

7.9(9) 

Na + Oj + X + NaO^ + X 

1.0(-32) 

[B] 

8.4(1) 

1.8(3) 

1.3(5) 

1.2(7) 

9.1(8) 

1.8(10) 

Na + H + X + NaH + X 

4.0(-32) 

[K] 

6.2(8) 

2.6(8) 

3.4(8) 

7.0(9) 

7.0(9) 

2.9(11) 

2.9(11) 

6.3(12) 

6.3(12) 

3.0(13) 

3.0(13) 

Production of Na 









NaO + 0 ■* Na + O 2 

4.0(-ll] 

[KJ 

6.U-1) 

2.5(3) 

4.2(-l) 

4.2{-l) 

4.2(-2) 

4.2(-2) 

8.9(-2) 

8,9(-2) 

2.5(-l) 

2.5(-l) 

4.2(-l) 

4.2(-l) 

NaO + 0^ Na + 20rj 

2.0(-12D 

[M] 

1.1(2) 

3.9(1) 

1.4(3) 

1.4(3) 

5.0(3) 

6.5(2) 

5.0(6) 

6.5(4) 

5.0(6) 


NaH + 11 Na + II 2 

4.01-17) 

[Ml 

1.7(9) 

1.5(8) 

1.9(8) 

5.4(8) 

5.4(8) 

2.5(9) 

2.5(9) 

6.9(9) 

6.9(9) 

8.9(9) 

8.9(9) 

NaK + OH Na t H^O 

2.0C-1S) 

[K1 

3.9(8) 

6.5(8) 

6.5(8) 

1.8(8) 

1.4(11) 

1.4(10) 




NaH + 0 + Na + OH 

2.0(-lS) 

[MJ 

1.2(4) 

5.0(7) 

8.3(3) 

2.9(4) 

8.3(2) 

8.3(2) 

1.8(3) 

1.8(3) 

5.0(3) 

5.0(3) 

8.3(3) 

8.3(3) 

Other 









Na02 + 0 NaO t 0,, 

3,0(-ll) 

[Ml 

8.0(-l) 

3.3(3) 

S.6(-l) 

2.0(0) 

S.6(-2) 

S,6(-2) 

1.2(-I) 

1.2(-1) 

3.3(-l) 

3.3(-l) 

S.6(-l) 

S.6(-l) 

NaOj + O3 * 1 - NaO + 2O2 

1.0(-14) 

[B] 

2.2(4) 

7.7(3) 

2.8(5) 

2.8(5) 

1.0(6) 

1.3(5) 

1.0(9) 

1.3(7) 

1.0(9) 


Na02 + H - NaH + O2 

3.01-12) 

[K] 

2.2(4) 

2.0(3) 

2.6(3) 

7.3(3) 

7.3(3) 

3.3(4) 

3.3(4) 

9.3(4) 

9.3(4) 

1.2(5) 

1.2(5) 


3 -1 » 6 -1 j . 

*Two and three body reactions have units cm s and cm s respectively, 
[M] = {Megie and 19771, [K] = [Kuifte, 1973], [B] = [Bfoun, 1973]. 


PAGE IS 
©I QUALITY 


215 


0RIG1NAI3 PAGES 13 
OP POOR QUAXOT 


Table II. 4 


Winter time constants for neutral sodium reactions. 
The entry format is the same as in Table II.l. 


Reaction 


Rate Constant* 


Time Constants (s] with altitude (km) 


Loss ot Na 



70 

SO 

90 

100 

110 

120 

Nu *■ O3 NaO + 0^ 

6, SC-12) 

[M] 

2.6(2) 

1.5(1) 

7.0(3) 

2.0(3) 

2.6(4} 

1.5(3) 

1.5(6) 

4.3(4) 



Na + 0 1- X NaO + X 

7.0C-33) 

[K] 

4.8(6) 

7.1(11) 

4.8(7) 

2.2(8) 

4.0(7) 

4.0(7) 

5.3(8) 

3.3(8) 

4.1(9) 

4.1(9) 

2.0(10) 

2.0(10) 

Na + O2 + X Na0.j + X 

l.OC-32) 

[B] 

1.6(2) 

3.8(3) 

1.3(5) 

3.8(6) 

8.7(7) 

1.3(9) 

Na + li + X Nail ^ X 

4.0C-32) 

tK] 

2.7(9) 

1.1(9) 

1.4(9) 

1.3(10) 

1.3(10) 

2.1(11) 

2.1(11) 

2,8(12) 

2.8(12) 

1.9(13) 

1.9(13) 

Production of Na 









NaO + 0 ->■ K;a f 0^ 

4.0(-n) 

fKj 

1.7(0) 

2.S(S) 

3.3(0) 

1.5(1) 

5.4(-l) 

S.4(-l) 

6.9(-l) 

6.9(-l) 

2.3(0} 

2.3(0) 

2.8(0) 

2.8(0) 

NaO + 0^ Wa + 2O2 

2.0C-12) 

[M] 

8.3(2) 

5.0(1) 

2.3(4) 

6.5(3) 

8.3(4) 

S.0{3> 

5.0(6) 

1.4(5) 



NaH + II -»• Na 

4.0C-17) 

[M] 

5.4(9) 

4.2(8) 

5.4(8) 

1.0(9) 

1.0(9) 

2.5(9) 

2.5(9) 

8.9(9} 

8.9(9) 

1.5(10) 

1.5(10} 

Nall + OH -► Na + II2O 

2.0C-15) 

[K] 

2.3(9) 

2.9(9) 

2.9(9) 

5.0(8) 

2.9(11) 

2.3(10) 




Nall + 0 Na + OH 

2.0C-J5) 

IM] 

3.3(4) 

5.0(9) 

6.5(4) 

2.9(5) 

1.1(4) 

1.1(4) 

1.4(4) 

1.4(4) 

4.6(4) 

4.6(4) 

5.7(4) 

5.7(4) 

Other 









NaO^ + 0 NaO + O2 

3.0C-11) 

IM] 

2.2(0) 

3.3(5) 

4.3(0) 

2.0(0) 

7.3(-l) 

7.3(-l) 

9.3(-l) 

9.3(-n 

3.0(0) 

3.0(0) 

3.8(0) 

3.8(0) 

NaO^ + 0^ NaO + 20^ 

l.OC-14) 

IB] 

1.7(5) 

1.0(4) 

4.6(6) 

1.3(6) 

1.7(7) 

1.0(6) 

1.0(9) 

2.8(7) 



NaO^ + H -I- Nall + 0^ 

3,0{- 12) 

IKJ 

7.. 3(4) 

5.6(3) 

7.3(3) 

1.3(4) 

1.3(4) 

3.3(4) 

3.3(4) 

1.2(5) 

1.2(5) 

2.0(5) 

2.0(5) 

3 •'1 6 ■'1 

*Two and three body reactions have units cm s’ and cm s' 
IM] = [Hegie mvl Blainont, J977], [K] = [Kvifte, 1973], [B] - 

respectively, 
[Braum, 1573], 
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described in Section 7.4.5 only these two reactions will be used to describe 
the neutral chemistry. 

Some question exists as to whether these two reactions proceed by the 
electron jump mechanism, in which case the rate constants are about two orders 
of magnitude larger [Koth and EZgin^ 1976]. The ablation model will still be 
useful in this case, if much of the deposited Na is converted to Na before 
it reaches the layer. This idea is pursued in Section 7.6. In the absence 
of experimental data, we will adopt those values obtained from kinetic theory 
analogies with hydrogen reactions. 

1 1. 2 Metal Ooside Ion Reactions 

11. 2.1 Introduction. Until recently the ionic reactions of the sodium 
chemical scheme have concentrated on the relation between Na , NaO , and Na02 . 
However most of the reaction rates have never been measured and some of the 
reactions may be speculative. For this reason it is important to evaluate 

the ideas behind this ion chemistry; the foundation for this work was pre- 
sented by Suiider [1969, 1970]. 

11. 2.2 Ionization processes. The primary source of Na^ is photoioni- 
zation of Na during the daytime, 

Na + ?iv ->■ Na^ + e" J=1.7xl0^s^ (II*1!) 

This value of the photoionization rate is quite typical for most studies of 
sodium which consider the ion chemistry. Snider [1969] obtained the value 
from the integrated product of the sodium absorption cross section and the 
unattenuated solar flux. The value is an upper limit, since 100% ionization 
efficiency is assumed. The calculation is of course sensitive to the reso- 
lution in the cross section and flux data and also the size of the integra- 
tion step. Photoionization of Na by HLya is negligible since unity optical 
depth of the solar ionizing flux for Na occurs near 40 Ion. 
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Charge exchange with NO'*’ and 0,* is of secondary importance in the day- 
time but dominates at night over photoionization, 

** 1 

( 11 , 2 ) 


Na + Na* + 0- 


, ^ ,--10 3 

= 6 X 10 cm s 


Na + no’*' ■+ No'*' + NO 


-11 3 -1 

= 8 X 10 cm 3 


( 11.31 

Ratos of these reactions are taken from Fcm\xghcv ct at* [I960]. Since the 
olQctron affinity for 0^'' and NO^ is about 12 ev and 9 ov I'espectively, while 
the ionization potential of Na is about 5 ev, those reactions are exothermic. 
The excess energy released by the reaction may appear as vibrational excita- 
tion in the diatomic product molecule or may even dissociate the molecule; 
it is possible that the excess energy may instead excite the sodium ion. 
Sodium oxide ions may be reduced by atomic oxygen by the process 

(II. 4) 


NaO^ + 0 ■> Na^ + 0- 


1 ,--10 3 -1 

ft e 10 cm s 


The analogous reactions involving MgO and SiO have been observed in the 
laboratory 1972] and are found to be quite fast [Ferguson and 

Fchsenfeldt 19oS]. Therefore it is likely that process (11.4) is similarly 
fast. If this is so, then processes involving NaO formation must bo considered 
as potential sources for Nn^. 

Since the ionization potential for NaO is thought to bo about 7.6 ev 
[Si^ider^ 1969] it is unlikely that photoionization of NaO is a major source 
of NaO^. However, the charge cxchajtgo px-ocesses 

NaO + 0^'' -X NaO'*’ + 0, (II. 5) 

NaO + no'*' NaO"^ + NO (II. 6) 

are exothermic. Yet as discussed in Section II, 1, the [0]/[0,] ratio keeps 
[NaO] small above 85 km. Since NO'*' and 0^^ begin to increase above 8S km. 
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if we assume the same rate for processes (II. 5) and (11,6) as we did for 
(II. 2} and (II. 3) respectively, it is unlikely that charge exchange processes 

■ 4 * 

with NaO are a significant source of Na . 

Reaction (II. 2) may have the branch 

Na + 02 "^ -V NaO'^ + 0 (II. 7) 

with the rate constant of reaction (II. 2) being about 5 times that of (II. 7) 

•f' 

[Ferguson and Fehsenfeld^ 1968]. However since NaO is returned quickly to 
Na^, reaction (II. 7) is probably not that significant. The analogous reac- 
tion involving NO^ is probably endothermic. 

Sources of Na^ resulting from Na02 will not be considered because infor- 
mation is lacking on the bond energy and ionization potential for NaO.,. How- 
ever if NaO is not thought to be an important source of Na , as discussed 
above, Na 02 is most likely of lesser significance since [Na 02 ] is thought to 
be much less than [NaO] [Broun, 1973], 

Photodissociation of the sodium oxide ions may also be a source of Na , 

NaO^ + Tiv ->■ Na^ + 0 (II. 8) 

NaO^^ + hv Na^ + 0., (II. 9) 

but experimental evidence does not exist for these reactions; these reactions 
may be invoked to create a diurnal variation in Na* [Blamont and Donahue, 
1964]. 

II. 2. 3 Loss processes of Na"^. Sodium ions may be deionized by the 
reactions 


Na* 

+ e~ -»• Na + /iv 

(11.10) 

Na* 

+ A~ Na + A. 

(iiai) 


Direct recombination \nth electrons is so unlikely that it can be ignored; a 
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1 2 3 •* 1 

typical rate constant is given as = 3 x 10** cm s [Bates and Dalgax>no, 

_7 

1962]. Mutual neutralization is much faster with a rate constant fc = 10 
3 - 1 

cm s [Brown, 1973], but concentrations of negative ions are not well 

kno\m. Certainly above 85 km. A" should fall off rapidly; [A"] during the 

3 -3 

day may be as high as 10 cm [Brown, 1973] . This reaction was also con- 
sidered to explain diurnal variation in ionization [Hunten^ 1967]. 

Indirect processes involving various forms of oxygen have also been 
suggested as being important in the deionization of Na Two three body 
reactions have been proposed vvhich form sodium oxide ions, 

Na'*' + 0 + X + NaO'*’ + X (11.12) 

Na'*' + O 2 + X ^ Na02'*' + X (11.13) 

These reactions are fairly slow at 90 km having time constants of about 10 
hours; rate constants are given for reactions (11.12) and (II. 13) as 
k = 10 ”^*^ cm^ 5 "^ and fe = 5 x lO”^^ cm^ s~^ respectively by Brown [1973] 
based on the experimental evidence of Ferguson and Fehsenfeld [1968], It 
should be emphasized that Na 02 may be unstable and therefore would break up 
easily in collisions with other atmospheric molecules; furthermore the reaction 

Na02‘^ + 0 -»■ NaO'*' + O 2 (11.14) 

is possible exothermic since the NaO -0 bond is weaker than the 0-0 bond 
[Swider^ 1969], Brown estimates the rate of this reaction ss k ~ 10~^® 
cm s” which would imply a substantial proportion of Na 02 is convei’ted to 
NaO^. Anderson and Barth [1971] have discussed the possibility of an analo- 
gous reaction involving Mg and suggest that its rate is = 10 cm^ s 
Since NaO returns to Na rapidly, as mentioned above, these loss processes 
for Na^ may not be that significant. 
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T\(fo body reactions have also been proposed as deionizing Na^ 


Na"^ 

+ Oj + NaO’*’ + O 2 

(11.15) 

Na^ 

+ 0 ->■ NaO^ + ftv 

(11.16) 

Na"^ 

+ 0^ ->■ Na02* + ho 

(11.17) 


The reaction involving ozone is possibly endothermic. Similar reactions 

4 " 4 * + 

involving Mg , Ca , Fe were all observed to be fast but the reactions in- 

+ 4 - -11 

volving Na and K produced no products, implying a rate constant /c < 10 

[Ferguson and Fehsenfeld^ 1968] . Unfortunately the thermochemistry does not 

illuminate the problem; for NaO^, the bond energy is Z?Q(Na^-0) = 1.25 ± 0,7 ev 

[Hilderd>rand and Murad^ 1970]. The reaction is exothermic if ZJ^CNaO^) > 

Z?qC 02 - 0 ) = 1 ev and hence the thermo ciiemis try is ambiguous. In any event the 

return of NaO to Na is so fast that the deionization of Na by 0^ is 

unimportant. 

Sijyider [1968] suggests that (11,16) and (11.17) are important near 100 

km with three body reactions dominating below 90 km. An appropriate value 

-17 3 -1 

for the rate constant is 10 cm s . Although this reaction may be of 
value in the neutralization of sporadic E layers it is debatable whether the 
reaction can be appropriately applied to the main Na layer where high atomic 

4 * 

oxygen concentrations will dissociate NaO . 

1 1. 2. 4 D-'ssoaiative reconibination of metal oxide ions. Dissociative 
recombination by the processes 


NaO^ + e~ ■+ Na + 0 

(11.18) 

^ * °2 

(I I. 19a) 

NaO + e 


iNaO + 0 

(II. 19b) 
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- 7 -6 

are believed to be fast having rate constants on the order of 10 to 10 
3 -1 

cm s . Based on the discussion in Section 7,4.1, the product NaO of 
reaction (II, 19b) is rapidly returned to Na at the altitudes of interest. 

The principal weakness of this ion scheme centers on the lack of know- 
ledge for rates of reactions and the existence of certain r(5actions. Fur- 
thermore the rapid return of metal oxide ions to Na via reaction (11,4) 
places a severe constraint on the theory. Tlie severity of the problem is 
indicated below. 

11,2.5 The general ateadystate solution based on reactions discussed 
in Sections II. 2 and II. 2, Figure II. 2 schematically displays the chemical 
processes involving sodium compounds which might be a priori of importance 
in the 80-120 km region. One is struck immediately by the symmetry of the 
scheme. In addition, an examination of time constants at 90 km reveals that 
the rapid return of NaO to Na and NaO’ to Na^ produces essentially two closed 
loops, one neutral ant ion ionic, with the dominant daytime connection be- 
tween these loops being photoionization of atomic sodium and dissociative 
recombination of the ionic species with electrons. 

The rates of these reactions are given in Table II. 5. Since reactions 
involving negative ions were assumed to be less important and cluster ions 
were not initially included, these ions were eliminated from Brown's scheme. 
Assumed daytime profiles are shorn in Figure II. 5. Other constituent values 
were taken from CJflA [1965] to parallel Brovm's work. 

Calculations were performed by setting the complete time expression for 
each constituent to zero. Molecular and eddy diffusion were not taken into 
consideration. The resulting set of simultaneous equations were solved for 
relative concentrations by setting [Na ] = 1 and then all constituents were 

A 

scaled appropriately so that [Na] = 1000 cm at 90 km. 
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s Table II. 5 

Rate constants for sodium reactions. 

Reaction Rate Constant 


Reaction 

Na + 0 / -* Na"" + 0, 


* 10'^° cm^ s'^ 


2, 

Na + 

NO+ Na* + NO 


8 

it 

10-" 

3. 

Na + 

P* Na* + P 


2. 

>4 

* 10"® 

4, 

Na + 

hy -t- Na"*” + e' 


J 

= 

1.7 X 10 

5. 

Na + 

0^ + X H. NaO^ + X 


1 

X 

10'^^ cm' 

7, 

Na + 

0^ -♦ NaO + 0^ 


1 

X 

io-13 

9. 

Na 

0 + X -V NaO + X 


1 

X 

10'33 

10. 

Na + 

0^* ^ KaO* + 0 


1 

X 

10-^® 

11, 

Na"^ 

+ e -*■ Na 


1 

X 

10-12 

12, 

Na**^ 

+ A“ Na + A 


1 

X 

10'^ 

13, 

Na* 

+ ©2 + X ^ Na02* + 

X 

5 

X 

10-52 

14, 

Na* 

* N, * X NaNj^ + 

X 

1 

X 

10-51 

IS, 

Na* 

+ 0 + X -► NaO* + X 


1 

X 

10-5® 

16. 

Na* 

+ Oj NaO* + O 2 


1 

X 

10-12 



4- + 




-12 

17, 

NaO^ 

+ X -*• Na + Ot + 

X 

1 

X 

10 

18. 

NaN-^ 

* + X ^ Na* + N 2 + 

X 

1 

X 

10-15 

19. 

NaO^ 

* + 0 ^ NaO* + O 2 


1 

X 

10-1® 

20, 

NaN2 

+ 0"*'Na0 ***^2 


1 

X 

10-1® 

21. 

NaN2 

+ 

+ e Na + N 2 


■ 2 

X 

10-" 

22, 

NaO^ 

* + e"^ Na + O 2 


2 

X 

10-® 

23, 

NaO* 

+ 0 ^ Na* + O 2 


1 

X 

10-1® 

24. 

NaO* 

+ e ■+ Na + 0 


4 

X 

10'^ 

26, 

NaO* 

+ A NaO + A 


I 

X 

10-2 

27. 

NaO^ 

+ 0 NaO + O 2 


1 

X 

10-11 

29. 

NaO 

+ 0 Na + O 2 


1 

X 

IQ-ll 

30, 

NaO 

+ ©2 NaO + O 2 


5 

X 

10-1® 

31, 

NaO 

+ NO* -* NtiO* + NO 


5 

X 

10-1® 

33. 

Na02 

+ Oj + NaO + 2 O 2 


1 

X 

10-1^ 

36. 

NaO 

+ 0^ Na + 2 O 2 


1 

X 

10-1® 


Note: Branching ratio between reactions 1 and 10 has been taken as 20%, 

Reactions and rate constants are taken from [1973], 
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Results o£ these calculations are shown in Figure II. 4. In interpreting 
Figure II. 4 it should be noted that profiles for these constituents are not 
being plotted. Rather the curves at each altitude give the predicted rela- 
tive proportions of the respective constituents. It is obvious though that 
Na^ should dominate at all altitudes, and this is not observed experimentally. 
Furthermore it was fovind that inclusion of the negative ions and cluster ions 
aided the falloff of [Na^]/[Na^] below 90 km. While NaO^, Na02^, and NaN 2 *^ 
were small in number density, the reactions returning these constituents to 
Na^ were sufficiently fast that they had to be included in the scheme. Finally 
it should be noted that the restriction [Na ] 1 creates the problem that a 

small change in Na^ can make a great change in the other minor constituents. 

It was obvious that the very completeness of Broom’s scheme with all of 
its uncertainties made it quite difficult to obtain agreement with experi- 
mental results. 
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APPENDIX III 
CLUSTER ION REACTIONS 

111.1 Expe:nmentat Studies 

Over the past ten years mass spectrometer measurements of ion-molecule 
reactions have been extended to higher pressures and to a wider range of tem- 
peratures, As a result it has been shown that third order kinetics play an 
important role in ion-molecule reactions. In many instances the clustering 
reactions of the form 

m'*' + A + X ^ M’^-A + X (III.l) 

where M* is an ion, and A and X neutral species, have been studied for the 
case of alkali ions. These reactions have the advantage of employing a less 
complex ion which can be copiously produced. Minimum values for rate constants 
or equilibrium constants for atmospheric molecules clustering to ions can be 
obtained by studying the clustering of alkali ions. As a result laboratory 
chemists can understand the clustering process more fully; and ionospheric 
chemists, interested in the metallic species, benefit directly. Therefore 
unlike the metal oxide ion chemistry discussed in Section II. 2, several clus- 
tering reactions have been measured for Na or may be easily extended to the 
Na case. If clustering appears to be important for the sodium chemistry, then 
much of the chemistry no longer needs to rely on speculation, 

111. 2 The Clustering Meehanism 

The mechanism of clustering is often considered as an energy transfer 
process [Goodj 1975]. Reaction (III.l) is believed to proceed in a two step 
process, ten^orarily forming an intermediate complex which may either break 
up or be stabilized by the third body of the reaction. This process is indi- 
cated by the following two reactions 
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+ A 


m'^-a* + X m'^-a + X 

lr 

a 


(III. 2) 


(III. 33 


where 

k ~ combination rate constant 
a 

k^ - decomposition rate constant 

k^ = stabilization rate constant 

k = activation rate constant. 
a 

Usually assuming equilibrium, the overall forward and reverse 

rates of this process may be written in terms of the above rate constants 

K (III.S) 

IVhat is interesting about cluster ion reactions is that they exhibit an 
apparent negative Arrhenius energy temperature dependence. Although these 
reactions are more complicated than common three body reactions, we may apply 
the general form for a reaction rate coefficient. The forward rate coeffi- 
cient can be written as [Good, 1971] 

k^ - aiT/T^)^ (III.6) 

and the reverse rate may be written as 

k^ = o{TlT^^ expC-e/D (111.7) 

where a, c, and e are constants and 2*^^ is a reference temperature, 
usually 300®K. Reimvt [1976] develops a method for estimating the errors 
incurred in the determination of these constants based on thermodynamical 
data. Because the necessary information for alkali clusters is not available. 
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the error analysis fails. However, the basic temperature dependence may 
still be determined. Equations (III, 6) and (III. 7) indicate that the tempera- 
ture dependence of the forward rate constant is the same as that of the pre- 
exponential factor of the reverse rate coefficient. The exponent, b, is 
negative as mentioned above and for Dq(M-A) » nRTf b is related to the num- 
ber of vibrational degrees of freedom in the cluster; that is |bj increases 
as the number of degrees of freedom increase [Durden et 1969], Although 

b has been shown to be large in absolute value for N0^»N2 formation [Arnold 
and Krankowskyi 1977], it is thought that b is between -1 and -2 for reac- 
tions involving Na [Keller and Beyer^ 1971] . As a rule of thumb may 
increase by a factor of 2 from 300°K, where the reaction rates are measured, 
to 200®K, appropriate for mesospheric temperature. 

Often when information is lacking specifically about k^ and k^, the 
equilibrium constant 

K = k^k^ (I II. 8) 

is measured. This quantity is independent of the third body and strongly 
temperature dependent. pressure units then K will be de- 

noted hv K I if they are in concentration units, then K will be denoted by K , 
p a 

III. 3 Determining the Equilibrium Constant for Cluster Ion Reactions at 
Mesospheric Temperatures 

The change in the Gibb’s free energy, AG, for a reaction has long been 
recognized as a criterion for the equilibrium and spontaneity of the reaction. 
By intuition it is thought that heat of reaction, Aff, or alternatively the 
change in enthalpy, is responsible for the spontaneity of the reaction; that 
is, if a reaction is exothermic, then it will occur. For an isothermal reac- 
tion the relation 
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AG = Aff - TLS 


cm. 9) 


applies where 

T “ absolute temperature 
A5 = change in entropy. 

Hence it can be seen that both concepts apply for TAG small. For cluster 
reactions, !TA5 is relatively small, and hence either AG or Afl will predict the 
spontaneity of the reaction. 

The equilibrium constant may be related to either the change in the Gibb's 
free energy or the change in enthalpy. A common expression for the variation 


of Jfp with temperature is given by 
-AG = flrin if 


cm. 103 


However AG may be a function of temperature, and hence equation Cm. 103 is 
not useful for extrapolating from room temperature, where the reactions are 
measured, to mesospheric temperatures. However since aG is independent of 
pressure, equation Cm.lO) predicts that is likewise independent of tem- 
perature. 

The change in enthalpy may be related to the equilibrium constant by 
the van't Hoff relation. 


In 


K Crj 
P 1^ 


-Aff 

R 




cm. 113 


where R - universal gas constant. 

There is evidence that hR is relatively constant with temperature for 
most cluster ion reactions over a 100“K range of temperature [Heimevt^ 1976], 
and this assumption is certainly true for hydrates of Na as shown in Figure 
(III. 13 . Hence equation Cm*23 niay be used to determine the value of at 
mesospheric temperature. Furthermore it is important to note that -Aff > 0 
for cluster ion reactions, and so K will increase as the temperature of the 


I 


2 

l/T*IO^Dog.K*)-' 


Figure III.l Van't Hoff plots for equilibrium constants 

•^n-l.n reactions Na'''*CH20] i + H^O + H2O 
Na"*"* [H203 ji + H2O [Dsidia and Kehctnle^ 1970 ], 
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reactants is lowered. 

One final point must be made about the equilibrium constant. Normally 
Xp is expressed in standard- state units of atm However for aeronoiny the 

3 

equilibrium constant, K^, has units cm . It is important to recognize the 

difference bet^veen K and K , By definition K is given for the reaction 

pa' P 




A + B -*■ C + D 


( 111 . 12 ) 


as 


p(? 

C D 

K - r (at equilibrium) 

P ^ ^ 

B 


(III. 13) 


T tfl * 

where P. is the partial pressure of the i — gas and j represents the number 

of times the i — gas appears in the stoichiometric equation. Assuming an 

ideal gas, P. = n.RTfV = a .RT , where e. is the concentration of the gas, 

we can rewrite (111,13) as 

4 


r'’ (f 

P rff. fh 

^A ^B 


r- K (RT) 


Av 


(III. 14) 


where Av is the number of moles of products minus the number of moles of 
reactants. For cluster ion reactions, Av = -1. Hence when the equilibrium 
constant is given in units of cm"^, should have an extra factor of T; 

that is if fej. varies as T~^ then should vary as [Rei.d, 1977], The 

point is a minor one because the factor of T may be assumed into the activa- 
tion temn of the reverse rate coefficient. 


